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A.  STRUCTURAL  ANALYSIS  OF  HOMOGENEOUS- 
AXISYMMETRIC  RADOME  SHELLS  UNDER 
MECHANICAL.  THERMAL.  AND  INERTIAL  LOADING 
1.  INTRODUCTION 

Up  to  the  present  time  shell  analyses  carried  to  the  point  of  numerical 
calculation  have  been  restricted  to  membrane>type  approximations,  shallow- 
shell  type  approximations,  or  to  shell  analyses  limited  to  axisymmetrical 
loadings.  The  present  development  combines  the  techniques  of  tensor  analysis, 
differential  geometry,  and  variational  calculus  with  the  capabilities  of  large- 
scale  digital  computing  machinery  to  permit  solution  of  the  shell  problen^s 
encountered  in  radome  applications  without  restriction  on  gecmetry  or  form  of 
loading  distribution,  and  with  comj^ete  consideration  of  both  membrane  and 
flexural  effects.  The  variational  approach  pemdts  treatment  of  static  and  • 
dynamic  or  vibratory  loadings  by  related  methods  and  permits  introduction  of 
thermal  effects  (both  material  property  variatiott  and  initial  stress 
effects)  without  major  complication. 

The  present  report  describes  the  formulation  of  the  radome  problem  for 
homogeneous  radome  shell  types  with  particular  emphasis  on  ogive,  cone,  and 
hemispherical  shapes.  The  techniques  developed  are,  however,  equally 
applicable  to  other  shapes  with  completely  general  loading  distributions  and 
subject  to  arbitrary  temperature  variations.  Within  the  framework  of  the 
tensor  methods  used,  the  internal  displacements  (and  velocities  in  dynamic 
problems)  are  expressed  as  first-order  tensors,  stresses  and  strains  in  terms 
of  second-order  tensors,  and  stress- strain  relations,  in  terms  of  fourth-order 
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tensors,  all  in  curvilinear  coordinates  appropriate  to  the  radome  shell  middle 
surface.  All  elastic,  thermal,  and  in  dynamic  applications,  inertial  aspects 
of  the  problem  are  formulated  in  general  tensor  notation  as  is  the  differential 
geometry  of  the  shell  middle  surface. 

Actual  numerical  computations  are  rendered  feasible  by  utilizing  the 
above-described  intrinsic  shell  theory  in'  conjxmction  with  an  energy  or 
variational  method.  Final  computation  reduces  to  the  solution  of  an  algebraic 
problem,  non-homogeneous  in  the  case  of  response  to  static,  dynamic,  or 
thermal  loading,  and  homogeneous  in  the  case  of  free  vibration.  The  entire 
development,  for  certain  conditions,  proceeding  from  the. prescription  of 
initial  shell  geometry  and  loading  to  the  final  resvilting  deformation  and  stress 
distribution  is  to  be  programmed  for  execution  upon  the  IBM  704  digital 
computer. 

The  present  section  describes  the  analytical  phases  of  an  investigation 
vdiose  ultimate  objective  is  the  development  of  an  efficient  and  accurate 
program  for  the  structural  analysis  of  thermally  and  mechanically  loaded 
axisymmetric,  homogeneous  radome  shell  structures  by  means  of  high-speed 
digital  computing  equipment.  The  end-product  sought  is  a  program  which 
reqviires  but  a  minimum  of  input  data  describing  the  shell  geometry  and  the 
imposed  thermal  and  mechanical  loading  distributions  and  \^ich  utilizes  this 
data  to  generate  solutions  of  the  structural  problem.  The  solutions  are  to 
include  adequately  detailed  descriptions  of  the  distributions  of  deformation  and 
stress  throughout  the  shell.  Moreover,  the  program  is  to  be  sufficiently 
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general  to  permit  inclusion  of  varying  degrees  of  edge  fixity  as  well  as  the 
presence  or  absence  of  an  elastic  central  boom. 

The  analysis  is  based  on  a  general  formulation  of  the  axisymmetric  shell 
problem  in  curvilinear  coordinates  and  takes  advantage,  from  the  very  outset, 
of  the  knowledge  that  the  problem  is  to  be  solved  numerically.  This  permits 
an  approach  which  relies  upon  a  minimum  of  physical  assumptions  or 
approximations.  The  only  such  assumptions  introduced  are  the  following: 

(a)  The  shell-type  assiunption  that  the  normal  stress  component  in  the  direction 
normal  to  the  shell  middle  surface  may  be  neglected  relative  to  the  other 
induced  stresses:  (b)  The  shell-type  assumption  that  normals  to  the  middle 
surface  remain  straight  under  deformation;  (c)  The  calculation  of  the  slopes 
of  the  above  normals  under  the  assvunption  that  normals  to  the  middle  surface 
remain  normal  under  deformation;  (d)  The  assumption  that  Young's  modulus  . 
c  and  Poisson's  v  may  be  considered  constant  throughout  the  shell  and  equal 
to  their  respective  mean  values  over  the  temperature  range  experienced  by  the 
shell  for  any  given  problem.  This  last  assumption  may  be  dropped  when 
axisymmetrical  loading  is  considered,  but  results  in  an  essential  simplification 
for  asymmetric  conditions. 

The  fact  that  it  is  recognized  in  advance  that  numerical  methods  are  to  be 
used  makes  it  unnecessary  to  use  any  of'the  standard  shell  approximations  of 
expanding  in  powers  of  the  normal  coordinate,  and  neglecting  various  terms  as 
small.  Incidentally,  it  is  this  type  of  approximation  v^ich  accounts  for  the 
major  portion  of  the  discrepancies  noted  between  the  linear  shell  theories  of 
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2 

various  investigators  . 

In  view  of  the  current  interest  in  the  unified  antenna- radome  concept, 
the  analysis  presented  herein  is  made  sufficiently  general  to  encompass  the 
case  of  the  truncated  axisymmetric  radome  with  central  boom.  The  formula¬ 
tion,  however,  is  completely  general  and,  aside  from  certain  problems  which 
arise  at  the  apex  of  a  complete  shell,  suffers  from  no  essential  complications 
when  applied  to  general  shells  of  revolution  with  or  without  a  central  boom. 

The  technique  employed  herein  hinge*  on  the  energy  approach  to 
structural  analysie  and  require*  first  an  extension  of  the  energy  principle 
to  permit  inclusion  of  thermal  effects.  For  this  purpose  the  so-called 
*'free  energy"  of  Hemp^  is  utilised  and  is  specialised  from  the  general  thermo¬ 
elastic  three-dimensional  form  given  by  Hemp  to  the  form  appropriate  to  ahell 
analysis.  The  use  of  the  energy  technique  together  with  a  general  tensorial 
formulation  of  the  energy  expression  in  curvilinear  coordinates  appropriate 
to  the  shell  middle  surface  virtually  eliminates  all  possibility  of  error  due  to 
neglect  of  terms  or  inaccurate  or  incomplete  resolution  of  forces  required  by 
equilibrium  methods.  Correspondingly,  the  macroscopic  stresses  or  stress- 
resultants  of  the  usual  shell  theory  are  by-passed,  the  entire  problem  being 
formulated  in  terms  of  displacements.  Furthermore,  the  difference  approxi¬ 
mation  required  by  the  numerical  solution  technique  is  viewed  as  a  direct 
method  in  the  calculus  of  variations:  The  difference  approximation  is 
introduced  directly  into  the  energy  integral  (thereby  converting  this  integral  to 
a  sunrunation)  and  the  equilibrium  problem  reduces  to  a  minimum  problem  with 


w;dd  tb  .^9-2? 


664 


respect  to  the  mesh-point  values  of  the  displacements.  This  technique  has 
several  advantages:  First,  it  insures  that  the  numerical  problem,  in 
difference  form,  has  a  positive -definite  matrix,  thereby  guaranteeing  the 
convergence  of  iterative  solution  techniques.  Second,  it  makes  the  introduction 
of  quite  general  boundary  conditions  extremely  simple  and  ijearly  automatic, 
and  finally,  it  eliminates  any  necessity  of  performing  the  operation  of  first- 
variation  to  obtain  the  Euler  equations  and  boundary  conditions  of  the 
variational  problem,  followed  by  differencing  of  these  equations.  Rather,  the 
differencing  in  the  energy  expression  proper  provides  the  matrix  of  the 
numerical  problem  directly  and  automatically  insures  the  appropriate 
expression  of  boundary  conditions  in  difference  form  whether  these  Correspond 
t*o  fixed,  ’partially  fixed,  or  free-edge  situations. 
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2.  GEOMETRY  AND  COORDINATE  DEFINITION 

Figure  1  shows,  for  the  general  shell  of  revolution,  the  coordinate 
system  utilized  together  with  the  associated  covariant  metric  tensor  (a^j) 
components.  The  metric  tensor  elements  are  expressible  completely  in 
terms  of  meridian  slope,  p,  and  the  two  principal  radii  of  curvature  Rj  and 
Rj  .  The  orthogonality  of  the  coordinates  insures  the  diagonal  nature  of  the 
metric  tensor  and  the  simple  relations  for  the  contravariant  form  a^^.  The 
quantity  J,  which  is  equal  to  the  square  root  of  the  determinant  of  aij,  gives 
the  volume  element  for  the  y^  coordinate  system  and  will  be  needed  for 
formulation  of  the  energy  integrals.  The  coordinate  y^  gives  the  inward 
distance  locally  normal  to  the  shell  middle  surface,  y^  gives  meridional 
distance  along  the  middle  surface  from  the  apex  of  the  middle  surface 
(extended,  if  the  actual  shell  is  truncated),  and  y^  is  the  longitude  angle  in 
radians . 

Figure  2  illustrates  the  same  data  specialized  to  the  case  of  the  trun¬ 
cated  cone.  In  this  case  the  variable  radii  of  curvature  R|  and  R^  are  re¬ 
placed,  respectively,  by  y^  tan  X.  and  od  ,  while  the  generally  variable 
meridian-tangent  angle,  j3,  becomes  simply  the  constant  semi-vertex 
angle  ^  .  The  metric  elements  simplify  correspondingly. 

The  contravariant  form  of  the  displacement  vector  corresponding  to 
the  y^  coordinate  system  is  denoted  by  and  is  to  be  distinguished  from 

the  physical  displacement  vector  (u^  u^  u^  ) .  The  reason  for  introducing 
both  tensorial  and  physical  components  of  the  displacement  vector  is  that 
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the  strain-displacement  relations  are  natural  in  the  tensor  form,  whereas 
the  shell-assumption  of  u^  independent  of  y*  while  u*  and  u*  vary  linearly 
with  y*  is  expressed  in  terms  of  physical  displacements.  The  entire 
problem  will  be  cast  ultimately  in  terms  of  physical  displacements. 

The  relation  between  physical  and  tensorial  displacements  is  well- 
known; 

u^  =4^  ^Ta...  =  4; 

^  n  ’  1 

not  summed  on  i 

It  is  to  be  noted  that,  in  view  of  the  orthogonality  of  the  coordinate 
system,  there  is  no  need  to  distinguish  between  contravariant  and  covaridnt 
forms  of  the  physical  displacements  u^. 
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3. 


STRESS-STRAIN  AND  STRAIN-DISPLACEMENT  RELATIONS 


The  covariant  form  of  the  strain-displacement  law  for  any  coordinate 


system,  within  linear  theory,  is  given  simply  as 


11  LJ 


(2) 


wherein  commas  denote  covariant  derivatives  and  the  are  merely  the 
covariant  components  of  the  displacement  tensor.  The  are  the 
covariant,  linear  strain  components. 

The  stress-strain  law  for  an  isotropic  material  in  general  curvilinear 
coordinates,  including  thermal  stress  effects,  is  given  by 

O'  /T*  ^  \  i j 


t'J  =  — ^  (a^^  ajl  +  a^j  a^S 

1+  V  1-  2v  l-2v 


^kl  -  To)a^J  (3) 


wherein  t'**  is  the  contravariant  form  of  the  stress  tensor,  c  is  Young's 
Modulus,  V  is  Poisson's  Ratio,  a  is  the  linear  coefficient  of  thermal 
expansion  and  (T-Tq)  is  temperature  relative  to  the  temperature  T^ 
associated  with  the  unstressed  state.  Repeated  indices  imply  summation 
throughout  unless^otherwise  specified. 
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4.  SHELL  ASSUMPTIONS 

At  this  point,  Ithe  shell  assumptions  are  introduced.  These  may  be 
summarized  as  follows: 

(a)  0  (4) 

i.e.,  the  normal  stress  (normal  to  the  middle  surface)  is  small  of  higher 
order  relative  to  other  stress  components. 

(b)  u‘  =  vaj^j(y*.yM  +  y*  (y\y*) 

s 

(y» ,  yM  +  y*  (y»  ,  y2)  (5) 

=u|^,{y^yM 


i.e.,  the  in-plane  displacements  are  assumed  linear  in  y*  while  u*  is 
assumed  independent  of  y^. 

F  inally 

(c)  the  ^\z)  computed  on  the  assumption  that  normals 

to  the  middle- surface  remain  normal  at  the  middle  surface. 

This  process  is  clearly  discussed  by  Novozhilov  .  It  is  to  be 
emphasized  that  the  foregoing  three  assumptions,  (a)  (b)  (c),  together  with 
the  use  of  constant  mean  values  of  c  and  v,  are  the  only  ones  introduced 
into  the  analysis  until  the  difference  method  is  introduced  to  permit  numeri¬ 
cal  solution.  No  intermediate  approximations,  such  as  the  common  expansiorj 
of  contravariant  metric  elements  in  powers  of  y^  ,  are  utilized. 
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The  assumption  (4)  of  vanishing  T^^is  utilized  to  eliminate  t|jj  from 
the  stress-strain  relation  (3)  .  Hereafter,  Latin  indices  i,j,k,  I  ...  will 
range  over  1,  2,  3,  while  Greek  indices  a,  p.'y,  6  ...  will  range  over  only 
the  middle-surface  indices  I,  2.  The  elimination  of  t^jj  from  the  general 
stress-strain  law  (3)  gives  the  following  shell  stress-strain  relations  for 
the  stress  components  parallel  to  the  middle  surface: 


+_il_  a‘^‘^aV^) 


1  +  V 


1  -  V 


€  a  a“^  (T  -Tq) 

1  -  V 


This  is  the  stress-strain  law,  including  thermal  effects,  which  governs  the 
general  shell-pf- revolution  problem. 

Next  is  invoked  the  assumption  (c)  of  normals  at  the  middle  surface 
remaining  normal  for  purposes  of  computing  the  .  This  is  equivalent 
to  assuming  that 


Expanding  the  indicated  covariant  differentiation. 


this  becomes: 


(8) 
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where  the  bracketed  expressions  are  the  Christoffel  symbols  of  the  second 
kind  based  on  a^j  and  where  the  well-known  symmetry  of  these  symbols  in 
the  lower  indices  has  been  utilized.  It  is  easily  ascertained  for  the  present 
metric  that  all  Christoffel  symbols  with  more  than  one  index  equal  to  3 
must  vanish.  Hence, 


Utilizing  relations  (ll  and  Che  fact  that  a**  «  •  t,  thic  may  be  expressed 


in  terms  of  the  physical  displacements 


no  sum  on  o 


New  the  relation  between  ChrietoCfel  symbols  of  the  first  and  second  kinds 
if  given  by 


due  to  the  esniohing  of  Christoffel  symbols  srith  mors  than  one  index  equal  • 
to  3.  Invoking  oncf  again  the  diagonal  character  of  a^j.  and  the  defining 
relations^  for  the  Christoffel  symbols  of  the  first  kind; 


C 
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the  relation  (10)  may  be  put  into  its  final  form; 
a 


8u 

8y3 


=  u 


(1) 


1 


8yJ 


8u^ 

"iy 


(14) 


-l„3 


y^  =  0 


no  sum  on  a 
a  =  1,  2 


This  pair  of  equations  serves  to  eliminate  from  all  further  considerations 

the  quantities  u)^..  in  favor  of  the  middle- surface  deflections  u®  .  and 

•  •  Uj  (o) 

gradients  of  u|^j  .  With  this  relation,  the  direct  consequences  of  the 
shell-assumptions  are  concluded. 
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5.  EXPRESSION  OF  STRAINS  IN  TERMS  OF  PHYSICAL  DISPLACEMENTS 
Equations  (1)  and  (2),  together  with  the  well-known  properties  of  the 
Christoffel  symbols,  permit  the  expression  of  the  strains  in  terms  of 

physical  displacements  u^.  These  take  the  following  form,  for  the  case  of 
the  metric  tensor  corresponding  to  Figure  1; 


’In  " 


8y' 


da 


u‘ 


11 


1  9a 


2 

u3 


9y" 


8y3 


“  u3 


u3 

^22  =  -r-j-  +  —  , 

dy^  2  9y3 


(15) 


*»ya,,  8ui 


12 


8y^  4  9y^  2  9yl 


These  strain  expressions  may  then  be  given  in  terms  of  middle- surface 
values  u^^j  of  displacements  by  simply  substituting  from  (5)  and  (14)  into 
the  general  displacements  u*  of  (15).  The  substitution,  while  routine,  is 
somewhat  lengthy  and  need  not  be  reproduced  at  this  point. 
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.  ENERGY  FORMULATION 


The  energy  formulation  sought  for  the  shell  problem  under  both  thermal 
and  mechanical  loading  may  be  developed  by  extending  to  general  curvilinear 
coordinates  the  "free  energy”  considered  by  Hemp  in  Reference  (1).  Hemp 
proposes  a  thermodynamic  explanation  of  this  quantity,  but  it  suffices  to 
consider  the  free  energy  as  simply  equal  to  the  variational  integral  whose 
Euler's  equations  provide  the  proper  equilibrium  equations  of  thermo¬ 
elasticity  as  well  as  the  proper  boundary  conditions.  In  this  sense,  the 
free  energy  serves  as  the  thermoelastic  analog  of  the  strain  energy  in 
ordinary  isothermal  elasticity,  and  this  is  the  context  in  which  it  will  be 
utilized  in  the  present  analysis. 

Hemp's  free  energy,  extended  to  general  curvilinear  coordinates  y^. 


Under  the  shell  assumptions  (a),  (b),  (c)  of  the  present  analysis,  the 


expression  for  F  may  be  re-derived  with  the  following  result: 
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This  relation  in  conjunction  with  a  work  integral  W  forms  the  basis  for 
the  shell  analysis  developed  herein.  The  work  integral  for  the  shell,  assum¬ 
ing  the  mechanical  loading  to  be  due  to  a  normal  pressure  distribution, 
p(y‘  •  y*).  is  given  by 

W  dyldyi  J(y^y^,o)  p  (y^  y*)  u^  (y*  y^)  (18) 

where  p  is  considered  positive  when  directed  inward;  i.e.  ,  p  is  actually 
the  local  value  of  (outer  pressure  minus  inner  pressure). 

The  combined  problem  of  mechanical  and  thermal  loading  of  a  shell  is 
then  stated  in  variaticwial  form:  The  solution  distribution  of  deflections  u^ 
is  such  that  (F-W)  is  stationary,  i.e., 

6  (F  -  W)  «  0  (19) 

Ul«  Iftcl  Uuit  linear  theory  is  used  throughout  implies  that  super- 

ponition  of  thermally  and  mechanically  induced  stresses  is  valid.  Hence, 
thd  tbennal  response  reeolting  from  the  temperature  term  in  F  and  the 
mechanical  responae  resulting  from  w  may  be  ca,lculated  separately  and 
added.  Care  must  be  taken*  however;  to  insure  that  the  proper  mean  values 
<ai  s  and  s  are  used  in  the  preaeure  problem  if  the  pressure  loading  is  to 
occur  under  high  temperature  conditions. 
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.  ELIMINATION  OF  y‘  -DEPENDENCE  BY  FOURIER  ANALYSIS 


The  thermoelastic  response  problem  defined  by  relations  (17),  (18), 
and  (19)  is  of  three-dimensional  character  in  its  present  form.  The  purpose 
of  the  present  section  is  to  reduce  the  problem  to  a  series  of  two-dimensional 
problems  in  which  the  y*  -dependence  is  lacking.  The  following  section  will 
describe  the  technique  for  reducing  each  of  these  two-dimensional  problems 
in  turn  to  a  one-dimensional  problem  by  integration  through  the  shell  thick¬ 
ness  with  consequent  elimination  of  y^  -dependence. 

The  reduction  to  two  dimensions  is  straightforward  and  is  accomplished 
by  Fourier  resolution  of  mechanical  and  thermal  loading  distributions  as  well 
as  of  the  corresponding  response  deformation  distributions.  Thus,  the 
independent  variable  y*  is,  for  reasons  of  clarity,  denoted  by  0  and  the 
pressure  load  p(y^y^)  is  expressed  as  a  cosine  series  in  0: 

OP 


p(yl  ,  y2)  =  p(y2 , 0) 


mP(y^  )  cos  mO 


UO) 


The  use  of  a  cosine  series  is  but  a  minor  restriction  since  the  origin  of  0 
may  be  chosen  at  will,  and  the  radome  loads  to  be  expected  will  generally 
exhibit  symmetry  about  some  meridian  plane. 

Similarly,  the  temperature  loading  may  be  expanded  in  a  cosine  series 
with  respect  to  0  where,  due  to  the  validity  of  superposition,  the  origin  for 
this  0  need  not  correspond  to  that  of  the  0  used  in  the  pressure  expansion, 
provided  the  responses  are  obtained  separately  and  superimposed  correctly. 
Thus, 
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where  the  coefficient  of  expansion,  a,  which  is  in  general  temperature- 
dep<»ndent,  is  included  with  (T-Tq)  to  simplify  the  analysis.  Corresponding 
to  loading  expansions  (20)  and  (21),  not  necessarily  considered  simultaneously, 
the  displacements  u^  are  resolved  into  Fourier  components  as  follows: 


3  3 

“(o)  =  "  - 


,u^  cos  mO 


m=  0 


,u^  cos  mO 


m=  0 


,u*  sin  mO 


m  =  1 


where  =  u^(v^) 

Now,  from  equations  (5)  and  (14),  together  with  the  fact  that  all 

geometric  properties  are  0 -independent,  it  follows  that  u^'^  ^  depends 

linearly  upon  u}  .  and  —  — ;  ,  both  of  which  are  given  by  sine  series. 

VO)  9y* 

2  2 

Similarly,  j  linear  in  u^^^  and  — —  ,  both  of  which  are  given  by 

cosine  series.  Hence  and  u^  are  expressed  as  cosine  series  in  0  while 
u*  is  a  sine  series,  which  state  of  affairs  corresponds  to  the  even  charac¬ 
ter  of  the  loading  (thermal  or  mechanical)  in  0. 

An  examination  of  strain  expressions  (15)  in  the  light  of  the  above 
remark  reveals  that  and  must  be  pure  cosine  series  while 

must  be  a  pure  sine  series. 
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Now,  from  the  diagonal  character  of  the  metric  tensor  a^^  (hence,  also 
of  a  which  is  merely  the  inverse  of  a-),  it  may  be  easily  seen  that  the  free 
energy  expression  (17)  involves  quadratic  forms  in  the  strains  of  only  the 
following  types; 

(■Hu)  (ll22  )  (ili2  )  (Till  ’Izz  ) 

and  that,  moreover,  these  expressions  are  multiplied  only  by  quantities 

independent  of  y^  (or  0).  Furthermore,  the  thermal  loading  term  in  F 

involves  only  products  of  a(T-TQ)  and  rin  or  t|22  together  with  0 -independent 

quantities  (recalling  the  assumption  of  a  constant  mean  value  for  c  and  v). 

Finally,  the  work  integral  W  involves  only  products  of  pressure  and  u^  , 

together  with  0 -independent  factors.  The  result  of  all  this  is  a  complete 

decoupling  of  distinct  harmonics;  i.e.,  associated  with  the  m*’^  harmonic 

of  the  loading  (thermal  or  mechanical)  there  correspond  only  the  m^^ 

harmonics  of  u*  u^  u^  and  of  u!  ,  uf  ,  and  u,  ,  and  of  the  strains  and 

(o)  (o)  (o) 

•  tresses.  This  is  a  simple  consequence  of  the  fact  that  F  and  W  involve 
integration  with  respect  to  0  over  the  interval  (0,  Ztf)  and  the  fact  that  the 
sines  and  cosines  form  an  orthogonal  set  with  respect  to  a  constant  weighting 
function  over  the  interval  {0,Z-it).  It  follows  from  the  relations 


sin^  mO  dO  = 


cos^  mO  dO  =  It 


m 


/O 


X 


Ztt 


(23) 


cos^  mO  dO  =  Zn 


m  =  0 
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that  the  free-energy  form  ma -  be  reduced  to  a  system  of  two-dimensional 


problems  independent  of  y*  or  0  by  introducing  the  simple  notation; 


2ir  C, 


■jr 


dy*  dy*  (  ) 


(24) 


(m)  m 

Likewise,  the  work  form  W  reduces  to  a  set  of  one-dimensional  integrals: 


*(m)  =  Cm  I  ■‘V'  < 


J 


) 


(25) 


with 


=  I 

m 

_  ± 
~  1 


m  =  0 
m  >0 


(26) 


wherein  the  m^“  integrand  involves  only  the  harmonic  of  pressure, 
temperature,  and  strain. 

The  problem  corresponding  to  (19)  for  the  m^^  harmonic  is  then 
simply: 


6  (F, 


-  W, 


(m)  (m) 


)  =  0 


(27) 


The  problem,  originally  three-dimensional,  is  thus  effectively 
reduced  to  a  set  of  decoupled  two-dimensional  problems.  In  practice,  it 
will  suffice  to  consider  only  zero*’^  and  first-order  harmonics  (m  =0,  1) 
unless  the  loading  exhibits  unusually  severe  0 -dependence  associated  with 
some  sort  of  load  or  temperature  concentration. 
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8.  REDUCTION  TO  ONE-DIMENSIONAL  PROBLEM  BY  INTEGRATION 


The  work  forms  of  expression  (25)  are  already  one-dimensional, 

and  it  remains  to  put  the  free-energy  harmonics  of  expression  (24) 

into  similar  form.  This  is  accomplished  simply  by  integrating  with  respect 
to  y*  through  the  shell  thickness  from  -h/2  to  +h/2.  The  integration  oper¬ 
ation  involves  no  essential  difficulty  for  the  portions  of  quadratic  in 

the  strains,  since  all  y* -dependence  is  completely  defined.  The  tempera¬ 
ture  contribution  to  is  treated  simply  by  defining  several  required 

weighted  integrals  of  oi(T-Tq)  through  the  shell  thickness.  In  view  of  the 
fact  that  the  distribution  of  a(T-TQ)  is  input  data  for  the  thermal  stress 
problem,  any  such  weighted  integrals  will  be  known  as  required.  It  should 
be  noted  again  that  no  approximations  or  expansions  in  powers  of  y*  are 
necessary  at  this  stage  of  the  reduction  in  view  of  the  fact  that  the  problem 
is  eventually  to  be  solved  numerically.  The  presence  of  variable  coefficient 
of  logarithmic  character,  for  example,  causes  no  difficulty  in  the  numerical 
solution. 


The  type  of  result  which  is  to  be  expected  from  this  operation  may 

be  indicated  as  follows:  Equation  (14)  gives  u|^j^  as  a  linear  combination 

of  u?  .  and  involving  no  y* -dependence.  Thus,  making  use  of  the 

(o)  dy'^ 

Fourier  resolution  of  the  preceding  section: 
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Now,  from  the  free-energy  expression  (17)  and  its  two-dimensional 


Fourier  resolution  (24),  it  can  be  seen  that  the  only  y*  -dependence  in 
^(m)  through  terms  of  the  following  types: 


J  Tiii)^  J  {a“  Ti22  y 


J  (^‘  Till  a“  TI22  ) 


J  a‘‘  a“  (1112)^ 


q(T-To)  Ja“  Till 


a(T-To)  J  a“  t], 


Furthermore,  examination  of  strain  expressions  (15)  shows  that 


a"  Til,  =  '/in' 


\dY  9  y^  J 


TI22  =  ^ 


a**  tii2  ■ 


{  ^  9  ^ 

\  9y*  9y*  J 

t'a^^  f  ^1 -  9u*  u*  9 


a  ^yi^7  . —  au*' 

- +  ^22  — 

9y^  ay» 


Moreover,  J  is  given  simply  by  '^slu  ^zz 

Therefore  the  terms  (29)  which  are  quadratic  in  the  strains  are  all  of 
the  form  _ 


N/a« 


fall  a“ 


multiplied  by  polynomials  in  y^ 
of  up  to  4th  order 


The  temperature  terms  of  (29)  involve  only  the  temperature,  or  a(T-To) 


multiplied  by  polynomials  in  y*  of  second  order  or  lower. 
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Now  terms  of  the  types  occurring  in  (31)  are  all  easily  integrated 
with  respect  to  y^,  resulting  in  nothing  more  formidable  than  logarithms 
whose  arguments  may  depend  upon  y^.  The  temperature  terms  may  be 


handled  completely  by  defining  the  following  three  weighted  averages  of 
the  general  m^^  harmonic  of  uIT-Tq)  (recalling  expression  (21)  ). 


These  weighted  averages  are  all  known  in  terms  of  prescribed  data 

th 

for  the  thermal  stress  problem.  Thus  the  problem  (27)  for  the  m 
harmonic  pf  the  displacement  distribution  and,  correspondingly,  of  the 
stress  distribution,  is  reduced  to  a  one-dimensional  problem  with  y^  , 
the  middle- surface  arclength  along  the  meridian,  Serving  in  the  role  of 
independent  variable.  As  a  final  parenthetical  remark,  it  may  be  noted 
that  the  logarithmic  terms  resulting  from  the  y*  -integration  will  all  be 
of  the  form 


a  =  1  or  2 


and  hence  are  well-behaved  as  long  as  the  local  semi-thickness  —  is 
smaller  than  the  smaller  of  the  two  local  principal  radii  of  ^ 

curvature,  Rj  ,  Rj  . 
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9.  BOUNDARY  CONDITIONS 

Only  one  set  of  conditions  remains  in  order  to  complete  the  problem 
statement,  and  this  is  the  prescription  of  boundary  conditions  at  the  for¬ 
ward  and  after  edges  of  the  radome  shell.  Within  the  framework  of  the 
energy  method  it  is  possible  to  prescribe  in  a  routine  manner  very*general 
boundary  .conditions  including  the  presence,  if  necessary,  of  an  elastic 
central  boom. 

Considering  first  the  base  of  the  radome  (y^=i)  the  shell  is  assumed 
simply  supported  but  with  a  distributed  moment  spring  of  spring  constant 
kj  inch-pounds  per  radian  per  inch  of  circumference  (see  Figure  3(a)  ). 

The  presence  of  the  spring  kj  corresponds  to  partial  fixity  and  may  be 
varied  from  zero,  corresponding  to  simple  support,  to  infinity  (see  the 
end  of  the  present  section)  corresponding  to  the  fully-clamped  condition. 

The  base  boundary  conditions  are  thus  specified  as  follows: 

at  y*  =f  ;  u*  *  _u*  =  _u’  =  0  (34) 

•'  m  m  m 

And,  to  the  free-energy  component  add  the  base-spring  energy: 
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Conditions  (34)  are  specified  in  the  conventional  manner,  whereas  the 


edge  moment  condition,  instead  of  requiring  the  usual  development  of  an 
expression  for  macroscopic  or  resultant  moment  in  terms  of  displacements, 
appears  in  the  form  of  a  simple  addition  to  the  free  energy.  These 
conditions  complete  the  base  boundary  specification. 

The  nose  boundary  conditions  are  inserted  in  a  manner  similar  to 
that  illustrated  by  expressions  (35),  there  being  in  general  no  nose  con¬ 
straint  equivalent  to  the  fixed  conditions  (34) .  Figures3(b),  (c),  (d),  and 
(e)  illustrate  the  equivalent  springs  for  the  nose  boundary  conditions. 

Figure  3(b),  with  associated  spring  kj ,  represents  the  resistance  of  a 
central  boom  (and  nose  attachment  to  the  boom)  against  axial  stretching. 

The  mean  axial  displacement  at  the  middle  surface  for  the  nose  is 
given  by 

2it 

(u^  cosP  +  u*  sinp)  dy^ 

3 

« 

=ai 


qI#  cos  sin^l 


ly’  =  0 

/  =ai 


(36) 


Thus  the  energy  addition  appropriate  to  axial  displacement  of  the 


nose  is:  Add  nothing  to  Fjyj  for  fn  >0,  and  to  Fq  add, 

"n 

k. 


(qU^  cosP 


+  qU*  sinP)^ 

Jy*=0 


(37) 


y*  =  ai 
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If  there  is  no  axial  restraint  from  a  central  boom,  k|  =0,  while  nonzero 
kj  -values  give  various  degrees  of  axial  stiffness.  Similarly,  the  mean 
value  of  nose  transverse  displacement,  which  is  resisted  by  elastic  bend¬ 


ing  of  the  central  boom,  if  such  a  boom  is  present,  is  given  by 

cosj3)  cos  y*  -  u'  sin  y'J  dy* 


^  r^(u^  sinp-  u> 

'^0 

=  y  ^jU2  ainp-  jU*  cos  p  -  ^  u^ 


i 

y  =0 


y  =ai 


(38) 


y*  =  0 


y  =ai 


And  the  corresponding  energy  addition  is  (see  Figure  3(c)  )  zero  for 
m  /  1  and 


(39) 


where  k^  is  the  spring  constant  for  the  boom  in  transverse  deflection.  If 
the  boom  is  absent,  k^  =  0,  while  nonzero  k^  -values  correspond  to  varying 
degrees  of  boom  flexural  stiffness. 

The  central  boom  ,  if  present,  also  resists  a  mean  value  of  shell 
transverse  slope  at  the  nose  (Figure  3(d)  ).  This  slope  is  given  by 
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i—  -^6.  r  2 

08  p  9y^ 


jU*  COsP 


-3) 


y"=  0 


y  =  ai 


the  computation  of  which  is  facilitated  by  the  fact  (Figure  1)  that 


ap  =  _  1 


R*  • 


The  energy  addition  corresponding  to  this  is  simply 


"4  /  = 


for  m  =  1  and  zero  for  m  /  1 


Again,  in  the  absence  of  the  boom,  =  0. 

Finally,  Figure  3(e)  represents  the  nose  equivalent  of  the  base-edge 
clamping.  The  distributed  spring  kg  provides  the  equivalent  of  the  degree 
of  clamping  of  the  nose  edge  to  the  boom  or  to  any  closing  member  even  if 
the  boom  is  absent.  For  purposes  of  computing  the  contribution  to  F  from 
this  equivalent  spring  distribution,  it  is  necessary  to  have  an  expression 
for  the  y*  -component,  ,  of  physical  rotation.  This  is  given  by  the 


following  relation: 


V  =  0 


\dY^  Ri  J 

y’=  0 


The  contribution  to  F  due  to  distributed  nose  spring  kg  is  given  by  the 
integrated  square  of  the  difference  between  this  Local  rotation  and  the 
angle  (-acosy*)  corresponding  to  mean  rotation  of  the  nose  structure. 

R|  cosp]  ■  +  JIL  +  a  cos  y‘l  dy*  (43 

■  2  J  R,  J 


ly"=  0 
y^  =  af 
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This  is  easily  integrated  to  give 


(44) 


=  al 


y*  =  0 


where  a  is  given  by  (40),  and  terms  in  (44)  containing  d  contribute  only 
to  whereas  the  series  expression  contributes  to  for  all  m  =0,  1,  2,  .  .  . 

In  the  event  that  the  nose  edge  is  free  or  is  simply  supported  at  the 
connection  to  a  central  boom,  then  k5  is  simply  set  equal  to  zero,  whereas 
a  fully-clamped  nose  condition  is  achieved  by  lettinp  become  indefinitely 
large. 

The  approach  to  infinity  of  any  of  the  aforementioned  springs  ir. 
order  to  simulate  fully  fixed  conditions  may  be  handled  within  ?  numerical 
difference  scheme  by  scaling  certain  of  the  mesh-pcint  values  of  displace¬ 
ment  or  differences  thereof  in  such  a  way  that  numerical  data  remains 
bounded.  This  is  numerically  equivalent  to  inserting  the  constraints  direct¬ 
ly,  as  in  equation  (34)  for  the  base  edge  displacement.  If  this  approach  is 
considered  undesirable  due  to  other  features  of  the  numerical  solution  tech¬ 
nique,  then  the  fully  fixed  condition  may  be  placed  in  the  direct  form  of 
equation  (34).  In  other  words,  in  place  of  letting  kj  approach  infinity  and 
including  in  the  energy  expression  the  contribution  (35)  with  the  consequent 
necessity  for  scaling,  the  condition  of  full  fixity  is  equally  well  inserted 
by  setting 
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9t_u* 

m 

dy^ 


y 


i 


=1 


0 


(45) 


and  imposing  the  corresponding  constraint  on  the  numerical  difference 
problem. 
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10.  FINITE-DIFFERENCE  SOLUTION  TECHNIQUE 


The  complete  shell  analysis  is  now  formulated  as  a  series  of  uncoupled 

variational  problems  of  the  form  (27)  in  which  the  energy  expressions  ^^n) 

« 

and  consist  of  one-dimensional  integrals  over  the  shell  meridian 

plus  possible  boundary -condition  energy  additions  to  F|jj^ which  are  expressed 
as  quadratic  combinations  of  edge  displacements  and  derivatives  thereof. 
Starting  from  this  point,  a  variety  of  solution  techniques  is  available.  . 

One  might  use  the  Rayleigh-Ritz  method  of  expanding  unknown  deflections 
in  a  series  of  known,  fixed-boundary-condition  satisfying  functions  and 
perform  the  variational  operation  directly  by  minimization  with  respect 
to  the  displacement  expansion  coefficients. 

The  existence  of  edge  effects  which  exhibit  strong  gradients  in  the 
edge  regions  militates  against  this  method  as  does  the  inaccuracy  inherent 
in  utilizing  the  method  to  obtain  stresses  which  involve  derivatives  of  the 
displacement  expansion  functions.  Still  another  disadvantage  of  the 
Rayleigh-Ritz  technique  lies  in  the  tremendous  number  of  integrations 
which  must  be  performed  to  gene. ate  the  matrix  of  the  numerical  problem. 
Alternatively,  one  may  perform  the  variational  operation  analytically  to  • 
obtain  systems  of  three  ordinary  differential  equations  for  each  harmonic 
(each  value  of  m).  These  (Euler)  equations  may  then  be  solved  by  finite 
difference  techniques. 

There  are  several  disadvantages  associated  with  this  technique.  First, 
the  Euler  equationr  must  be  found  analytically,  as  must  the  natural  boundary 
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conditions  associated  with  free  or  spring-restrained  edges.  Then  both  the 
equations  and  the  boundary  conditions  must  be  placed  in  difference  form, 


an  operation  which  leaves  a  great  deal  to  personal  judgement.  Finally, 
whin  thij  ij  complete,  there  is  no  assurance  that  the  matrix  of  the  system 
of  equations  will  be  symmetric  and  positive  definite,  and  these  character- 
is  ics  are  helpful  in  guaranteeing  the  convergence  of  iterative  solution 
techniques^.  The  large-scale  character  expected  for  the  system  of  equa¬ 
tions  indicates  that  iterative  techniqr.ea  will  be  needed  for  refinement  of 
solutions  if  not  for  actual  solution  itself. 

The  approach  followed  here  is,  in  a  sense,  a  combination  of  the  above 
two  techniques.  The  continuous  range  (aJ  ,  f )  of  y*  is  replaced  by  a  dis¬ 
tribution  of  discrete  points  and  the  values  of  the  displacement  components 
m’^(o)  nV(o)  m^(o)  these  points  are  considered  as  degrees  of  freedom  for 
the  variational  problem.  The  point  spacing  need  not  be  uniform  and,  in 
fact,  in  order  to  describe  adequately  the  edge-effect  or  "boundary-layer" 
effect,  this  spacing  should  be  considerably  smaller  in  the  base  and  nose 
edge  regions,  changing  to  a  larger  value  for  the  major  internal  portion 
of  the  range.  The  decisive  point  here,  however,  is  the  fact  that  the 
difference  approximation  is  introduced  directly  into  the  energy  integrals 


{which,  in  turn,  become  sums)  and  the  boundary  energies,  thereby  convertin 

the  total  energy  expression  (E  ,  -  W,  .  )  into  the  sum  of  a  quadratic  form 

(m)  {m) 

and  a  linear  form  in  the  point  values  of  the  displacements  •  This 

serves  to  make  the  differencing  more  routine  and  less  a  matter  of  choice. 


but  more  important,  it  guarantees  the  symmetry  and  positive  definiteness 
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of  the  matrix  of  the  equation  system,  since  the  matrix  arises  from  the 
positive  definite  quadratic  term  in  the  free  energy.  Finite-difference 
approximation  of  this  term  should  not  alter  its  positive  definite  character, 
and  will  result  in  a  symmetric  system  matrix.  The  finite  difference  form 
of  the  Euler  equations  or,  what  is  equivalent,  the  equilibrium  equations 
is  then  given  simply  by  equating  to  zero  the  derivative  of  the  total  energy 
form  with  respect  cO  each  of  the  degrees  of  freedom.  Thus  the  matrix  of 
the  quadratic  energy  form  (multiplied  by  2)  becomes  the  matrix  of  the 
difference  equations,  a.id  the  vector  whose  inner  product  with  the  unknown 
displacement  gives  the  linear  form  in  the  total  energy,  becomes  the  forcing 
vector  for  the  system.  Hence  the  above-mentioned  differentiation  with 
respect  to  the  point-values  of  the  deflections  need  not  be  carried  out, 
since  the  results  are  known  beforehand,  and  the  matrix  of  the  quadratic 
energy  form  together  with  the  vector  of  the  linear  form  serves  to  describe 
completely  the  system  of  nonhomogeneous  linear  equations  which  con¬ 
stitutes  the  difference  form  of  the  equilibrium  problem.  The  finite- 
difference  equivalents  of  all  boundary  conditions  are  automatically 
incorporated  into  the  matrix  of  the  quadratic  form,  and  no  further  attention 
to  boundary  condition  specification  is  required.  This  statement  holds  true 
whether  the  boundary  conditions  correspond  to  free-edge,  partially  fixed, 
or  fully  clamped  situations,  although  in  the  latter  case  a  numerical  scaling 
transformation  will  be  necessary  if  the  full  fixity  constraint  is  not  incor¬ 


porated  as  an  absolute  constraint. 


In  order  to  illustrate  in  somewhat  greater  detail  the  type  of  difference 
replacement  which  is  introduced  into  the  energy  integral,  the  following 
exan,ples  of  typical  terms  should  serve  as  models.  Typical  of  the  quad¬ 
ratic  first  derivative  terms  which  appear  in  the  energy  expression  is  the 


following: 


wherein  w  is  the  generic  displacement  variable.  If  the  end-point  (af)  is 
numbered  n  =  U  and  the  point  i  by  n  =  N,  this  term  would  be  replaced  by 
the  difference  approximation. 


n,  n-1 


-  w„  , 
n  n-1  I 

Vn-1 


where  w^  is  the  value  of  w  at  point  n,  and  is  the  spacing  between 

points  n  and  n-1.  Form  (47)  holds  regardless  of  nonaniformities  in  point¬ 
spacing,  i.e,,  regardless  of  relative  magnitudes  of  the  spaces  n-1* 

A  term  of  the  form 


pw  dx 


which  typified  the  linear  forms  contributing  to  the  forcing  of  the  system 
would  be  replaced  by 


>  Pn  "n  ^'"'0 
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N,  N+1 


are  both  to  be  taken  equal  to  zero.  An  examina- 


where  A  and  ^ 

Uf*l  *,i*^i* 

tion  of  the  difference  forms  resulting  from  differentiation  of  these  expres¬ 
sions  with  respect  to  a  point  value  shows  that  proper  difference  fo,  mulas 
for  the  second  derivative  of  w,  even  for  the  case  of  an  irregular  mesh, 
result  automatically.  (The  second  derivative  of  w  is  the  proper  differen¬ 
tial  form  associated  with  the  variational  integral  (47)  ). 

For  terms  such  as  the  following: 


X'&) 


dx 


(50) 


which  typifies  the  highest  order  contributions  to  free  energy,  the  differ¬ 
ence  approximation  is: 

N-1 


2 


^  ^n,  n-l  ^^n,  n4l^ 


n=  1 


^n,  n-1  ^  ^n,  n-fl/  \^n,  n-1  ^n,  n+1 


(51) 


A  .i"'  i  +  A  ,w,,-(A  ,+  A  .,)  w_ 

n,  n+1  n-l  n,  n-1  n+1  '  n,  n-1  n,  n+l  n ; 


Once  again,  an  examination  of  the  difference  form  resulting  from 
differentiation  with  respect  to  some  internal  Wj^  shows  that  the  result  is 
precisely  the  centered  difference  form  (for  a  general  irregular  mesh) 
corresponding  to  the  fourth  derivative,  which  is  the  proper  Euler  expres¬ 
sion  corresponding  to  the  integral  (50).  Furthermore,  it  is  not  difficult  to 
establish  the  fact  that  differentiatio.n  with  respect  to  end  point  values  of 
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Wjj  in  (47),  and  both  end-point  and  next  interior  point  values  in  (48)  leads 
to  backward  difference  formulae,  for  the  right  end  of  the  range  (and  for¬ 
ward  difference  formulae  foi  the  left  end  of  the  range)  corresponding  to 
the  natural  bouncer^  conditions  associated  with  the  variational  expressions. 

Thus  the  replacement  of  the  integral  expression  for  “  ^(m)  ^ 

by  difference  forms  such  as  the  preceding  examples  (with  similar  expression^ 
for  mixed  types)  and  addition  of  difference  equivalents  of  the  bdundary  spring 
terms  serves  to  plrce  the  entire  problem  in  proper  difference  form  and, 
at  the  same  time,  insure  the  symmetry  and  positive  definiteness  of  the 
formulation.  Solution  for  the  response  deformations  then  consists  in 
solving  the  rystems  (one  system  for  each  harmonic,  .n,  considered)  of  equa¬ 
tions  for  the  point  values  of  the  displacements  rn'^(o)  *  systems  will, 

of  course,  involve  a  large  number  of  simultaneous  equations,  but  the  special 
character  of  the  matrix  which  can  be  achieved  by  proper  ordering  of  the 
unknowns  permits  efficient  solution  by  any  one  of  several  well-established 
numerical  techniques. 

The  solution  of  the  system  of  equations  will  be  expressed  in  terms  of 
the  displacement  components  at  the  middle  surface  at  the  discrete  number  of 
points  chosen  for  the  difference  approximation.  These  middle-surface  valueqj 
together  with  the  relations  for  the  slopes  ,  in  terms  of  these  middle- 
surface  displacements  and  gradients  thereof,  permit  calculation  of  the 
displacement  through  the  thickness  of  the  shell  at  the  specified  mesh  points. 
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Likewise,  availability  of  -he  displacements  permits  calculation  of  the 
strains  and  stresses,  with  required  gradients  replaced  by  the  appropriate 
difference  approximations.  Strains  and  stresses  may  be  obtained  by  means 
of  higher  order  difference  formulae  in  order  to  minimize  loss  in  accuracy 
due  to  the  difference  approximation  of  gradients. 


f»95 
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11.  CONCLUSIONS 

The  objective  of  the  development  reported  in  the  foregoing  sections  is 
the  construction  of  a  digital  computer  program  capable  of  generating  stress 
and  deflection  solutions  to  the  thermally  and  mechanically  loaded  axi- 
symmetrical  radome.  The  technique  developed  stems  from  a  general 
formulation  of  the  appropriate  energy  expressions  in  terms  of  deflection 
distributions  and  proceeds  to  a  difference  formulation  of  the  numerical 
problem  by  way  of  a  direct  variational  method.  The  energy  formulation  is 
sufficiently  powerful  to  permit  inclusion  of  effects  due  to  an  elastic  central 
boom  with  no  essential  additional  complexity.  Varying  degrees  of  edge 
restraint  are  also  capable  of  inclusion  without  difficulty.  Both  iterative  and 
exact  inversion  techniques  have  been  developed  for  solving  the  resulting  large 
scale  algebraic  problem  in  a  matter  of  a  few  minutes  on  a  high-speed  digital 
computer,  and  the  problem  is  organized  into  a  positive-definite,  symmetric 
form  which  insures  convergence  of  iterative  methods. 

The  programming  of  the  analysis  has  been  carried  out  for  the  cone  and 
ogive.  The  cone  program  has  been  written  for  the  closed  nose  condition  as 
well  as  the  truncated  cone  with  general  boundary  condition  capable  of  including 
a  central  boom.  The  ogive  is  treated  for  the  closed  nose  condition.  Details 
of  these  programs  are  given  in  a  subsequent  section. 

An  Appendix  B  is  included  with  the  present  report  to  indicate  an  approach 
based  on  physical  considerations,  to  the  apex  boxmdary  conditions  on  non- 
truncated  shells  of  revolution.  As  noted  in  the  appendix,  the  physical  apex 
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boundary  conditions  should  be  augmented  with  a  development  of  mathematically 
proper  apex  conditions. 

In  connection  with  the  problem  of  boundary  conditions,  however,  several 
important  points  should  be  kept  in  mind.  First,  some  of  the  physical  apex 
boundary  conditions,  or  combinations  thereof,  developed  in  Appendix  B  prove 
to  be  essential  for  boundedness  of  stresses,  rotations,  and  deflections  off  the 
shell  middle  surface,  as  well  as  of  the  free-energy  integral  itself.  Thus, 
these  conditions  must,  under  all  circximstances,  be  enforced  at  the  apex. 
Second,  the  use  of  direct  energy  methods  such  as  is  done  in  the  present 
analysis,  is  a  powerful  tool  in  the  handling  of  boundary  conditions.  Apart  from 
direct  constraints,  such  as  fixed  edges,  the  inclusion  of  the  free  energy  of 
potential  energies  due  to  all  sources  guarantees  the  satisfaction  of  all  natural 
boundary  conditions.  Thus,  the  incorporation  of  all  physical  apex  conditions, 
as  derived  in  Appendix  B,  which  relate  to  the  direct  physical  connection  of 
the  shell  with  itself  at  the  apex,  should  give  automatically  all  the  apex  boundary 
conditions,  those  not  forced  explicitly  being  implied  by  the  variational  procedur 
and  hence  appearing  naturally. 


APPENDIX  A  -  LIST  OF  SYMBOLS 


a-.j,  aij 


'm 


h 

k 

i 

P 


U'  11  ^ 


w 


yl  yJ 


z'  z* 


F 

J 

Rl.  Rz 

T 

W 

a 

a 

t3 

€ 

V 

6 

6 

X 

jiJ 

’’ij 

(j'j 

^n,  n+1 


fraction  of  meridian  at  which  truncation  occurs 
metric  tensor 

coefficient  -  1  for  m  =  0,  =1/2  for  m  >  0 
shell  thickness 
spring  constants 

meridian  length,  extended  to  apex 
pressure 

physical  deformation  components 
generic  symbol  for  u^u^u^ 
coordinates 

phyrical  rotation  components 
free  energy 

Jacobian  or  volume  element 
principal  radii  of  curvature 
temperature 
virtual  work 

coefficient  of  linear  expansion 

mean  rotation  at  shell  nose 

local  meridian-curve  slope 

Y  oung  '  s  modulus 

Poisson' s  ratio 

first  variation 

longitude  angle  =  y* 

cone  semi-angle 

stress  tensor 

strain  tensor 

physical  components  of  stress 

finite  difference  point  spacing  between  points  n  and  n+1 
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APPENDIX  B  -  PHYSICAL  BOUNDARY  CONDITIONS  AT  THE  APEX 
OF  A  NONTRUNCATED  SHELL  OF  REVOLUTION 
Following  is  a  discussion  of  the  physical  boundary  conditions  appro¬ 
priate  to  the  apex  of  a  nontruncated  shell  of  revolution.  In  a  sense,  the  uce 
of  the  term  "boundary"  is  incorrect,  since  the  closure  of  the  shell  at  the 
apex  implies  that  the  apex  is  no  longer  a  boundary  of  the  shell  but  is  an 
internal  point.  However,  with  respect  to  the  coordinate  system  utilized 
(Figures  1  and  2),  the  apex  constitutes  the  lower  limit  of  the  variable  y^ 
and  hence  a  boundary  of  the  ,  y^  domain.  The  term  "physical"  ic  used 
advisedly,  since  the  conditions  derived  herein  are  devised  on  phyfical 
grounds.  It  is  well-known  in  plate  and  shell  theory  that  the  physical  bound¬ 
ary  conditions  are  not  always  the  mathematically  appropriate  boundary  con¬ 
ditions.  The  question  of  mathematically  appropriate  boundary  conditions  at 
a  free  edge  (not  the  apex  of  a  closed  shell)  are  resolved  automatically  by  the 
use  of  a  variational  treatment;  i.e.,  the  natural  boundary  conditions  of  the 
variational  formulation  are  the  mathematically  appropriate  ones  and  there 
need  be  no  further  concern  about  edge  physical  conditions.  At  the  apex  of 
a  closed  shell,  however,  the  natural  boundary  conditions  cannot  be  expected 
to  provide  all  the  boundary  conditions  unless  the  apex  is  punctured,  and 
this  puncturing  process  admits  solutions  which  should  not  be  present  for 
the  closed  shell. 

The  following  development  of  physical  conditions  appropriate  to  the 
apex  of  a  closed  shell  is,  therefore,  only  part  of  the  story.  The  degree  to 
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which  these  conditions  coincide  with  the  mathematically  appropriate  ones 
will  require  further  investigation. 

The  types  of  shell  shapes  of  practical  interest  fall  into  two  categories; 
those  having  a  physical  singularity  at  the  apex,  such  as  cones  and  ogives, 
and  those  having  no  such  physical  singularity,  such  as  ellipsoids  and  hemis¬ 
pheres.  Within  the  coordinate  scheme  used  herein,  hov/ever,  it  must  be 
understood  that  a  singularity  is  present  no  matter  what  the  shape  of  the  shell. 
This  is  evidenced  by  the  singular  behavior  of  the  metric  element  a**  due  to 
the  presence  in  this  element  of  the  factor  l/cos^  |3  ,  for  the  blunt-nosed 
shell  and  due  to  the  fact  that  R|  — ^0  for  the  pointed  nose  shape.  Now  this 
fact  leads  to  the  preliminary  conclusion  that  certain  terms  in  the  free-energy 
expression  contain  nonintegrable  singularities.  This,  of  course,  is  an  im¬ 
possibility  since  the  free  energy  must  be  bounded.  It  will  turn  out  that  the 
physical  boundary  conditions  derived  herein  are  precisely  the  conditions 
necessary  to  keep  F  bounded.  It  should  be  emphasized  that  analysis  of  the 
apex  region  is  essentially  a  three-dimensional  problem  in  which  the  shell 
assumptions  are  invalid.  Nevertheless,  the  apex  conditions  derived  herein 
are  conditions  which  should  be  approached  in  some  sense  by  the  shell 
solution  as  the  apex  is  approached.  It  is  in  this  context  that  the  present 
appendix  is  to  be  interpreted. 

The  physical  boundary  conditions  at  the  apex  may  be  stated  as  follows: 

As  the  apex  is  approached,  y^ — >  0,  and  the  entire  range  (0,  2ir)  of  y* 
approaches  a  single  common  point;  i.e.,  as  y^  — >0,  all  values  of  y* 
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between  0  and  2ir  are  physically  coalescing  to  a  common  point.  It  follows 
directly  from  this  fact  that,  as  y*  —♦0,  both  the  physical  displacement 
vector  u  and  the  physical  rotation  vector,  which  will  be  defined  as  z  (the 
rotation  tensor  will  be  denoted  by  ,  and  its  physical  components,  by  z^), 
become  independent  of  y*  or  equivalently,  of  0.  This,  however,  is  pre¬ 
cisely  equivalent  to  the  statement,  following  Sokolnikoff^,  that  u  and  z 
become  parallel  vector  fields  with  respect  to  latitude  circles  as  y^ — ^  0. 
This  may  be  stated  mathematically  by  prescribing  the  condition  th*t  the 


intrinsic  or  absolute  derivatives  of  z  and  u  with  respect  to  0  must  vanish 
as  y^  — >0.  The  statement  then  becomes: 


as  y 


•J  d9 


(B-1) 


0  j  =  2.3 


=  1  j  =  1 


(B-2) 


Thus  (B-l)  becomes; 


as  y 


Expanding  the  first  of  these  (the  second  follows  by  analogy) 

z  .  .  r  '  1  J 


as  y 


I  j 


-ai_+ 

9y‘ 


[ij.k]  6^ 


(B-3) 


(B-4) 


VArD  TR  3 


701 


Expanding  still  further,  this  becomes: 


Replacing  the  tensor  quantities  by  the  physical  displacements  u^ 
by  means  of  (1),  the  above  relations  become; 


Now,  in  view  of  the  fact,  embodied  in  equations  (14),  that  the  dis¬ 
placement  components  u^  are  calculated  in  a  manner  implying  locally  rigid 
rotation  through  the  shell  thickness,  it  suffices  to  satisfy  relations  (B-6) 
on  the  middle  surface  y*  =  J.  Thus,  equations  (B-6),  specialized  to  the 
middle  surface,  y*  =  0,  constitute  three  of  the  apex  physical  boundary  con¬ 
ditions.  The  other  three  conditions  follow  from  the  second  of  equations  (B-3)i 
by  direct  analogy: 
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as  y 


^/a^■ 


—  +  ^ 

2  Vj?” 


\ay* 


0y3  ^ 


y  =0 


^  1^1  -  ±!1 
ft,.*  2 


n/7*^  z* 


y  =  0 


(B-7) 


8z^  i  / 

8y»  2  \^  8y’ 


z* 


V’  =  0 


Equations  (B-7)  take  on  a  greater  significance  when  the  actual  expres¬ 
sions  for  the  rotation  components  z*  are  known  in  terms  of  the  u^. 


)  8u’ 


These  relations  are  obtained  as  follows 

^  *  1  flk  .  11^  =  i  .  a  (u^  ^  ^ 

^  J  \8y*  dyy  J  \ay^  8y*  / 

2  .  1  /"a^i  aij^  1  /^8  (u*  nT^i,  )  8u’  \ 

'  ■  7i.i7  ■  T77  %■  — i7y 

.-L  Ail .  »ii\  =  i  i^L .  °t  ''■)  ^ 

^  J  \^8y'  8y‘y  J  “  3y‘  8  y'  J 

And,  since  the  relation  between  the  physical  and  tensor  components  of  z  is 


‘  =  4' 


(B-9) 


not  summed  on  i 
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it  follows  that: 


z'  = 


f 


(u*  ) 


z2  =  ^  ( -i-  ^ 

V  ay^  ay* 


z'  =  ^ 


l"  a“  ^ 


*  9 

.rr-  8(u*  ) 

^£L^>  "■  “ 


(B-IO) 


Now,  from  equation  (14) 


ay“  ^  8y3 


-  i  u“  -f  -S^.. 


0 


(B-ll) 


no  sum  on  a 


au  a 

ay“  ay* 


Nr?^u“_£iaa 


y*=  0 


(B-12) 


Hence 


'y*  =  0 


„  _  9a  9u 

-.  -  _ 

1  9y*  ay'* 


y*  =  0 


(B-13) 


Specializing  the  rotations  of  (B-10)  to  the  middle  surface,  the  y^ 
deflection-derivatives  which  occur  in  these  expressions  may  be  replaced 
using  (B-13)  to  give  the  final  relations  for  the  physical  rotation  components 


on  the  middle  surface. 
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=  0 


/du  \ 

VyV 


sf.iz 


9y 


(u‘  ) 


-^Y^  =  0 


Now,  the  first  set  {B-6)  of  boundary  conditions  has  an  alternative  interpre¬ 
tation;  If  equation  (14)  for  —  =  u*  is  examined,  it  is  seen  that  this 

OyJ  (*) 

quantity  appears  to  diverge  as  y^ - >  0.  For  rewriting  (14)  in  the  form: 

3 


(1) 


3  y' 


3 

3y 


S) 


(B-15) 


y  =0 


the  factor  implies  an  infinity  as  y^  - ^0  .  However  the  third  of 

boundary  relations  (B-6)  is  precisely  the  expression  needed  to  keep  u|^^ 
bounded  at  the  apex.  Without  going  into  the  details,  it  may  be  stated  that 
the  first  two  boundary  relations  (B-6)  are  precisely  the  requirements  needed 
for  boundedness  of  the  physical  stresses  at  the  apex  and  of  the  free-energy 
expression  as  well.  The  physical  stresses  (T  can  be  shown  to  be  related 
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to  the  tensor  stresses  T^j  by  the  expressions  : 


O'*'  = 


4-  t’ 

,11  11  JJ 


(B-16) 


in  the  present  case. 

(no  sum  on  i  or  j) 

The  six  general  apex  boundary  conditions  (B-6)  and  B-7)  are  next 
examined  for  the  individual  harmonics  (values  of  m)  under  the  Fourier 
resolution  of  Section  VII.  Considering  first  the  axisymmetric  harmonic, 
m  =  0,  it  is  seen  that 


ii  «  0 


ii-Ls  0 

8  y‘ 


(B-17) 


throughout  the  shell. 


In  consequence  of  these  relations,  it  follows  that 


=  0 


J  =  0 


(B-18) 


throughout  the  shell. 


Now  the  second  and  third  of  apex  conditions  (B-6)  vanish  identically  as  does 
the  first  of  apex  conditions  (B-7).  Thus,  for  the  axisymmetric  (zeroth) 
harmonic,  the  apex  conditions  reduce  to  two  relations; 


as  y  _ ^  0 


2  V3y  ° 


a  “va, 


y*=  0 


+  - —  u  J  - 

ayj  o  y 

-^y  =0 


(B-19) 
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The  second  relation  of  (B-19)  implies  that  as  y  - ^0 


(B-20) 


=  0 


for  both  pointed  and  blunt-nosed  shapes  since 


8^ 


ay" 


is  nonzero  at  the 


apex  in  both  cases.  The  relation  (B-20),  in  turn,  implies  as  - >  0 

,  a  (  ^  ) 


9o-^ 


ay" 


a  y" 


1: 

y  =  0 


0 


(B-21) 


Now  the  two  apex  conditions  may  be  expanded  to  give: 


as  y 


a 


J2 _  +  9„ 


ay^ 


(B-22) 


1 


Rj  cosp 


id  3 

o'i  - ^  (Ri  cos(3)  -  qU  cos  P 


ay^ 


y’=  0 


But  a  (Rj  cosp) 


a  y" 


is  equal  to  the  sine  of  the  semi-angle  p 


I  z 

y  -►o 


at  the  apex. 

Therefore,  the  second  of  relations  (B-22)  may  be  expressed  as  follows: 
as  y"  — >  0 

1 


Ri  cos  p 


sinp  -  qU^  cos  ^ 


(B-23) 


y  =  0 
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which  simply  states  that  the  radial  displacement  at  the  apex  is  zero.  The 
first  of  relations  (B-22)  states  that  the  slope-change  must  vanish  at  the 
apex.  In  other  words,  the  apex  angle  must  not  change  under  deformation, 
and  the  apex  material  must  not  open  up. 

Next,  the  implications  of  (B-6)  and(B-7)  for  the  first  harmonic 
{m  =1)  are  considered.  Relations  {B-6),  give  the  result: 


9 

9y" 


(B-24) 


0 


3 

y  = 


0 


Now 


3  V  a 


11 


9y' 


=  sin  p 


y^=  0 


and 


3 


11 


9y^ 


=  -  cos  |3 


y  =  0 


(B-?5) 
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So  that  relations  (B-24)  may  be  rewritten: 


Rj  cosp 


jui  +  lU* 


sin  p  -  j  u*  cos  p 


y  =  0 


+  j  u'  sin  p 


p  - >  0 

-l„J=  n 


(B-26) 


3  1  n 

1  U  -  1  u  cos  P 


-y"=  0 

But  any  one  of  these  three  relations  is  derivable  from  the  other  two,  so 
that  only  two  are  independent. 

Following  from  (B-7),  (B-10),  and  (22),  similar  relations  must  hold  for 
the  ,z^: 


as  y 


Rj  cos  p 


j  z*  +  1  sin  p  -  j  cos  p 


0 


j  z  ^  -  j  z*  sin  p  J  - ^  0 


(B-27) 


1  Z^  +  1  z*  cos  p 
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and  again  only  two  of  the  three  relations  are  independent. 

An  examination  of  expressions  (B-14)  for  the  at  y3  =0  indicates 
that  1  z^  and  jz*  become  indeterminate  as  y*  — ^  0.  For 


1  z 


2 

-  (  -  jU*  COSp  +  jU^) 

Rj  cos  (3 

o 

II 

>> 

- 

(B-28) 


y  =0 


1 


Rl  cosp 


y^=  0 


2  1  ®  1'^' 
(  -  jU  -  jU  sin|3)  -  - : 


9y' 


-^y^=  0 


However,  some  thought  on  this  matter  indicates  that  the  numerators  and 
denominators  of  the  indeterminate  forms  must  proceed  to  zero  at  the  same 
rate^  as  y^ - >  0,  so  that  these  forms  are  in  general  finite,  nonzero  quanti¬ 

ties,  and  two  of  expressions  (B-26)  and  two  of  expressions  (B-27)  serve 
as  the  apex  boundary  conditions  for  m  =  1. 

For  m  >  I  the  consequences  of  (B-6)  and  (B-7)  are  as  follows: 
as  y^ - ^  0 


1 


Rj  cos  fi 


{  m  j„u*  +  j„u*  sinl3  -  j.^uJ  cos^) 


y*=  0 
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“  rn  -  sin  3  ^ 

m  m 


ui  _ ^ 

-  y*=  0 


-  m  +  co8(3  - >  0 

Jv^=  0 


(B-29) 


with  analogous  expressions  holding  for  physical  rotations  z^. 


Now  if  the  second  and  third  relations  of  (B-29)  are  substituted  into 
the  first,  and  the  same  operation  is  performed  on  the  analogous  z^-equations, 
the  following  results  are  obtained: 


I  m  u^  -  (-i— )  uM  _ 
m  '  m'  m  I 

■"m*'  ->;k»  m»‘l  — j 

-1  5  _ 


(B-30) 


Hence,  as  y 


V  =  o 


(B-31) 


But  then  it  must  follow  from  the  second  and  third  relations  of  each  of  the 
_u^  conditions  and  the  j^jZ^  conditions  that,  as  y^  — ^0 


i  =  1.  2.  3 


i  =  1.  2.  3 


(B-32) 
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The  apex  physical  conditions  for  harmonics  higher  than  the  first  (for  m  >1) 
are  thus  given  by  the  condition  that  all  middle  surface  displacements  and 
rotations  approach  zero  as  the  apex  is  approached. 

The  expressions  for  the  may  be  obtained  from  (B-14) 


Once  again,  the  physical  rotations  ^z*  and  are  indeterminate  at  the 

apex. 

The  questions  which  must  be  answered  by  considerations  of  mathe¬ 
matically  proper  boundary  conditions  are  twofold.  The  first  relates  to  the 
indeterminacy  in  the  rotations.  The  second  relates  to  the  number  of  apex 
boundary  conditions  appropriate  to  each  harmonic.  It  would  appear  that 
there  should  be  three  apex  boundary  conditions  on  the  zeroth  harmonic  and 
four  for  all  higher  harmonics.  The  physical  conditions  lead  formally  to 
two  apex  conditions  on  the  zeroth  harmonic,  four  on  the  first  harmonic, 
and  six  on  higher  harmonics.  The  six  conditions,  however,  may  overlap 
due  to  indeterminacy  in  the  rotation  components. 
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The  determination  of  mathematically  appropriate  boundary  conditions 
at  a  singular  point  of  a  differential  equation,  such  as  the  apex  in  the  prefent 
problem,  is  discussed  formally  by  Friedman  .  The  mathematically  proper 
conditions  are  shown  by  Friedman  to  be  the  specification  of  the  conjunct 
(or  bilinear  concommitant)  of  the  desired  solution  with  the  initials  of  a^ 
linearly  independent  set  of  solutions  of  the  homogeneous  differential  equa¬ 
tions.  Actual  implementation  of  such  a  program  forms  a  considerable 
effort  in  itself  and  will  not  be  attempted  herein.  The  ultimate  results  of 
such  an  investigation  should,  hov/ever,  agree  generally  with  the  physical 
apex  conditions  derived  herein  and  should  provide  the  additional  informa¬ 
tion  needed  to  explain  the  counting  and  indeterminacy  questions  which  the 
^aysical  derivation  leaves  unanswered. 
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APPENDIX  C 


DYNAMIC  PROBLEMS 

In  dynarr.ic  problens  r.noh  as  vibrations  or  transient 
loading  (pressure  or  thernal)  It  is  necesaary  to  fonr.nlate  a 
kinetic  energy  integral  for  use  in  conjunction  with  the  free- 
energy  expression  (31)  or  with  the  simpler  strain  energy,  V, 
This  kinetic  energy  expression  is  easily  cast  in  the  tensor 
form  apnropriate  to  curvilinear  coordinates  and  takes  the 
folio’-ring  form  (dots  denoting  time-derivatives): 


T  =  i  p  J  dyMy^dy^  a^i  , 

=  V'  P  J  dy’dy^dy^  u^n^ 


(C-l) 


in  the  oar.^  r-f  oT^thogonal  coordinates  on]y» 

• 

v/heroin  p  is  ''boll  material  density,  ^  denotes  tensor  velocity, 
and  u  is  physical  velocity.  Tinder  vibration  at  a  sinr»le 
frequency  (o,  the  kinetic  energy  becomes 


T 


1 


(o 


J  dyMy^dyl  a^^ 

J  dy’dy^dy^  u^u^ 


(C-2) 


VADD  ra  59-22 


714 


The  .Tvailability  of  the  foregoing  energy  and  v/ork  forms 
permits  formulation  of  the  vari-^us  classes  of  problems  to  be 


encountered  in  radnme  anali cations .  The  underlying  orincinle 
in  the  formulation  is  Hamilton’s  or,  v.’hnt  is  equivalent  in 
the  direct  variational  formulation,  the  Lagrange  equations  of 
motion.  In  the  presence  of  thermal  effects,  the  nctential 
energy  is  to  be  replaced  by  the  free  energy.  Follcv.’ing  is  a 
classification  of  the  v'’riational  orinciples  to  be  used  for 
various  cases  of  interest. 

a.  Free  vibrations,  v/ithout  thermal  effects: 


/•ti 

6  (  (T-V)  dt  =  0  (C-3) 

where  b  indicates  first  variation  and  integration  vntli  resnoct 
to  time  is  performed  between  two  arbitrary  limits,  t^  ard  t,. 

b.  static  pressure  loading,  v/ithout  thermal  effects: 

6  (V-W)  =  0  (C-4) 


where  6W  =  ^^dy’dy^  J(y\y^,0)  p, 


(C-5) 


c.  Thermal  loading,  without  pressure  loading: 


6  (F)  =  0 


(C-6) 
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It  should  be  noted  that  case  (?»,)  is  actually  all-inclusive, 
provided  the  vanishing  of  T  in  static  cases  and  equivalence  of 
F  and  V  in  isothermal  cases  is  recognized. 

This  concludes  the  classification  of  the  various  problems 
which  fall  within  the  category  of  linear  analysis,  ’-/ithin  this 
limitation  of  linearity  it  is  permissible  to  superimpose  effects. 
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For  exDpipln,  the  Fhorrnl  v/hlch  '''.9kos  it;s  nopearnnce 

in  the  free-er.er^y  function,  F,  may  bo  interpreted  simply  as 
a  loading  equivalent  to  sore  prossu^e  distribution  and  vri  th 
sere  associated  boundary  conditions.  The  response  to  this 
thermal  load1n,y  may  be  calculated  separately  and  then  super¬ 
imposed  upon  the  response  to  other  loadings  such  as  pressure. 
This  linear  analys is,  ^-/ith  its  consequent  validity  of 
superposition,  breaks  drun  in  certain  cases  vhen  one  loading, 
be  it  thermal  or  otherv/lse,  induces  significant  membrane 
stresses  in  a  shell,  and  a  second  loading  then  acts  to  deform 
the  shell  in  a  locally  normal  direction.  This  set  of  circurr- 
stances  is  handled  by  an  initial-stress  type  of  analysis  and 
is  similar  to  situations  arisinf'  in  bucklinf:  problems  and 
centrifugal  stiffening  effects.  Its  treatment  requires  the 
introduction  of  nonlinear  strain-displacement  relations  and 
v/ill  not  be  handled  herein. 
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Metric  Tensor  for  General  Shell  of  Revolution 


ds*  =  dy^  +  (1  --X_)  dy*  +  (R^  -  y* )  cos^  P  dy* 
^2 

a. j  =  a^j  =  0  :  i  /  j 


3^11  ■  (^1  ■  y^)  cos^p  =  1/a” 

a,,  =  (l.-Xl)^  =  1/a^^ 

Ra 

a33  =  1  =  1/ a^^ 


J  =  '^la-  • 

'  ij 


= 


^22 


(Rl-  y')  (1 


y3 

n  )  cos  p 
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Metric  Tensor  for  Cone 


Rj  =  tan  \ 

Ri  =  ® 

p  =  \ 

aij  =  a'-*  =  0  :  i  /  j 
a22  i  1  =  l/a“ 

an  =  (y^  sin  X.  -  coB\f  -  l/a“ 

I  1/33 

a33  ~  ^  “  Va 
J  =  (y*  sinX  -  y^  cosX) 
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FIG.  4  TYPICAL  MESH  SCHEME  FOR 

CONICAL  SHELL 
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B.  TRUNCATION  ERROR  GROWTH  IN  EXPLICIT  DIFFERENCE 


SCHEMES  FOR  NUMERICAL  SOLUTION 


OF  HEAT  CONDUCTION  EQUATIONS 


SUMMARY 


An  analysis  is  made  of  the  possible  growth  or  accumulation  of  trunca¬ 
tion  error  for  explicit  difference  numerical  calculations  of  the  time -dependenl 
temperature  distribution  in  one-dimensional  or  slab  configurations  with 
constant-property  materials.  Two  limiting  cases  are  studied:  (a)  the  case 
of  homogeneous  or  zero  initial  temperature  values  and  step  heat  input  at 
one  wall  with  the  other  wall  insulated,  and  (b)  the  case  of  arbitrary  nonzero 
initial  temperature  distribution  and  both  walls  insulated.  For  both  cases, 
it  is  proven  that  as  long  as  stability  criteria  are  satisfied,  the  truncation 
error  does  not  accumulate.  Moreover,  since  any  slab  problem  involving 
constant-property  materials  with  variable  heat  flow  through  one  wall  and 
insulated  opposite  wall  may  be  formulated  as  a  superposition  (of  Duhamel 
type)  of  these  two  problems  (a)  and  (b),  it  is  apparent  that  truncation  error 


will  not  accumulate  in  general  providing  stability  criteria  are  satisfied. 
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The  interpretation  of  these  results  for  variable-property  materials 
must  be  done  somewhat  intuitively,  but  it  is  expected  that  they  hold  in  this 


case  as  well.  One  may  resort,  for  this  purpose,  to  the  argument  which 
replaces  the  variable  material  properties  by  local  or  instantaneous  values 
and  invokes  the  constant-property  result  in  the  small.  Proceeding  in  this 
way,  the  result  should  hold  in  the  large. 
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TRUNCATION  ERROR  GROWTH  IN  EXPLICIT  DIFFERENCE 
SCHEMES  FOR  NUMERICAL  SOLUTION 
OF  HEAT  CONDUCTION  EQUATIONS 


1.  INTRODUCTION 

In  the  finite -difference  computation  of  temperature  distributions 
through  radome  shells  under  the  influence  of  aerodynamic  heating,  the 
thinness  of  the  shell  together  with  the  values  of  specific  heat  and  conductivity 
exhibited  by  ceramic  materials  dictate  that  the  computations  be  performed 
using  a  rather  course  space-mesh  (dividing  the  shell  thickness  into  about 
six  equal  intervals)  relative  to  shell  thickness  and  a  time  step  chosen 
relative  to  the  space  mesh  such  that  the  stability  criterion  for  error  growth 
is  satisfied.  The  time  step  required  for  stability  may  be  so  small,  due  to 
shell  thinness,  that  for  a  typical  trajectory  of,  say,  ten  minutes  duration, 
as  many  as  20,000  time-steps  may  be  necessary  to  span  the  trajectory. 
(Incidentally,  this  requires  only  three  min’ites  of  IBM  704  time.)  The 
occurrence  of  nonlinear  phenomena,  such  as  radiative  boundary  conditions 
depending  on  the  fourth  power  of  surface  temperature,  dictate  that  the  ex¬ 
plicit  scheme  be  utilized  rather  than  an  implicit  scheme  which  would  permit 
larger  time  steps.  The  advantage  gained  by  the  use  of  fewer  time  steps 
would  be  more  than  offset  by  the  nonlinear  algebraic  problems  associated 
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with  the  implicit  scheme  in  the  presence  of  nonlinear  boundary  data  or 
nonlinearities  imposed  by  temperature -dependent  specific  heat  and/or 
conductivity. 

The  truncation  error  associated  with  the  relatively  coarse  space 
mesh,  and  the  large  number  of  time  steps  required  to  compute  a  trajec¬ 
tory  raise  the  question  of  the  meaning  of  such  a  computation.  Does  the 
accumulation  of  truncation  error  render  meaningless  such  a  calculation? 
Following  is  an  indication  that  in  at  least  two  limiting  cases,  the  trunca¬ 
tion  error  does  not  build  up  so  as  to  render  results  useless  after  a  large 
number  of  time  cycles.  The  two  cases  treated  are  (a)  the  case  of  constant 
heat  input  at  one  wall  with  the  other  wall  insulated  and  (b)  the  case  of  zero 
heat  input  or  output  at  both  walls  such  as  would  be  expected  after  radiative 
equilibrium  is  reached  during  a  flight  trajectory  calculation.  In  this  dis¬ 
cussion  only  truncation  error  is  considered  since  the  simple  difference 
scheme  and  the  small  number  of  space  intervals  utilized  raise  the  truncation- 
error  question.  Roundoff  error,  on  the  other  hand,  is  small  for  the  classical 
explicit  difference  scheme  and,  for  a  time  step  which  satisfies  the  stability 
criterion,  cannot  be  magnified  by  repeated  time  cycles.  There  still  remains 

I 

the  question  of  accumulated  roundoff  error,  but  this  is  not  treated  herein. 
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2. 


PROBLEM  A 


Consider  first  the  case  of  a  slab  at  zero  initial  temperature  subjected 
to  a  step  heat  input  QH  (t-o)  at  the  outer  wall  and  insulated  at  the  inner  wall. 
The  exact  solution  to  this  problem  is  the  following: 


c  \h 


-2  ch 
k 


_n*  ir^  ht/ch\ 


where  h  =  wall  thickness 

c  =  (specific  heat  )  multiplied  by  density 
k  =  conductivity 


(The  assumption  of  zero  initial  temperature  is  not  restrictive  since 
case  (b)  will  treat  the  problem  of  nonzero  initial  temperature  with  zero  heat 
flow  at  the  two  walls.) 

Now  as  an  indication  of  the  degree  to  which  the  exponential  terms  in 
the  above  series  (1)  decay  with  time,  it  should  be  noted  that,  in  the  difference 
method,  the  ratio  p  given  by 


where  At  is  time  step 

Ax  is  space  interval 
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is  of  the  order  of  0.5.  Thus  the  largest  exponential  term 


.-ir^kt/ch* 


may  be  estimated  after  N  time  steps  if  the  space  mesh  is  -  times 

M 


the  wall  thickness. 


-It*  k  t/ ch^ 


«  e 


N  time 
steps 


/  NkAt 

c(Ah)^ 


ir^  N  ^  _l 
M*  2 

e  (3) 


With  a  space  mesh  of  ^  times  the  thickness  of  the  wall,  M  *  6  and  after 


100  time  steps 


IT  kt/ch* 


TT^  X  IQO 
72 


.  13.7 


*  e 


(4) 


1 100  time  steps 
M  *  6 


This  term,  as  is  seen  from  (1),  is  to  be  compared  with  terms  of  the  order 
of  unity.  Thus  after  a  time  equivalent  to  100  time  steps  in  the  difference 
method,  the  exponential  terms  in  the  series  are  insignificant.  Hence  for 
times  greater  thsui  this,  the  difference  in  temperature  distributions  between 
two  times  tj  and  tj  (both  greater  than  100  At)  is  given  very  closely  by 


This  is  obviously  independent  of  x. 
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Next  the  difference  approximation  is  scrutinized.  The  difference 


scheme  appropriate  to  the  present  problem  is,  in  terms  of  the  sketch, 


2  3  4 


^  ^ 
till 

S  S  S  S 


n+l  n  ,  n  n  . 

+  2P  (Ti  -  To  ) 


n+l 

T 

m 

n 

X  T  + 
m 

..n 

m+l 

n+l 

T 

M 

n 

*  T  + 
M 

rT,n  V 

o 

T 

*  m 

X  0 

3 

H 

O 

z  ... 

where  d  =  - r-  as  noted  earlier 

c(A  x) 

This  may  be  expressed  in  matrix  form 

n+ 1  n 

T  *  A  T  +  b 
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where  A  * 


and  where  b  * 


and  where  T°  x 
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Now 


T»  »  AT  +  b 
T'  X  at‘  +  b 
T*  X  at*  +  b 


Ab  +  b  *  (A  +  I)  b 
(A*  +  A  +  I)  b 

etc. 


(11) 


Thus,  the  vector  T*'  at  the  time  step  is  simply; 


n  ,  n-l  n-2 
:  «  (A  +  A  + 


+  A  +  I  )  b 


(12) 


Likewise  the  difference  vector  between  the  temperatures  at  two  times 
nj  At  and  n2  At  is  given  by 


-  T  ' 


-  ^  (A 


n,-l  .nj-Z  .n,+l  .n,  ,  , 

^  +A*+....+A^  +A*)b 

_ _ i 


(13) 


(n^  -  n|)  terms 


Next,  it  is  necessary  to  examine  the  eigenvalues  ^  and  eigen  vectors  of 
the  matrix  A.  This  is  best  performed  by  considering  the  difference  equation 


P  +  (1-2  (3-X)  +  pq™+^ 


(14) 


m  X  l,  2 


M-1 


with  initial  and  final  conditions 


(1-2  P-X)  q°+  2pq^  .  0 


M-1  .  .  M  _ 

2p  q  +(l-2p-X)q  rQ 


(15) 


Now  let 


1  -2  P  -X 


X  -2  cos  0 


(16) 
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where  it  might  turn  out  that  d  ie  imaginary,  but  this  possibility  is  accepted 
for  the  present.  Then  the  general  solution  of  (14)  is 

q*”  *  a  cos  mO  +  b  sin  m0  (17) 

Introducing  (17)  into  initial  and  final  conditions  (15),  there  results: 

-  a  cosO  +  a  cosO  +  b  sinO  z  0 

(18) 

-  ccsO  (a  cos  M®  +  b  sin  M®  )  +  a  cos  (M-t)  ®  +  b  sin  (M-l)O  *0 

From  the  first  of  (18)  there  results: 

b  *  0  (19) 

since  9  >  0  leads  to  a  result  which  will  be  included  in  those  which  follow 
from  the  vanishing  of  b.  The  second  of  equations  (18)  then  becomes 


-  cos®  (cos  M9)  +  cos  (M-l)  9*0  (20) 

Or,  expanding  out  the  cos  (M-l)  ®  term, 

-  cosO  cos  MO  -I-  cos  MO  cosO  •¥  sin  MO  sinO  >0  (21) 

so  that  sin  MO  sinO  «  0  (22) 

The  proper  eigen  vectors  are  then  given  by 
r,  nir 

®n  *  -TT  "  *  °  I,  2,  3 - M  (23) 


where  it  can  be  shown  that  n  >  M  leads  to  eigenvalue  results  which  overlap 
those  of  (23);  e.  g. 

V  »  I,  2.  3 . 
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Now  it  it  easily  teen  from  (26)  that  for  m  >  0  the  eigenvalues 
are  lest  than  unity 

<  I  for  m  >  0  (32) 

Moreover,  for  |3  ^  y  which  it  necessitated  by  the  stability  criterion, 


>  -  I 

When  ^  <  y  the  inequality  holds  in  (33). 

Thus, 

for  p  <i 

I  <  I  for  m  >  0 

For  example,  when  M  *  6,  p  *  0.4 

i 

Xj  «  I  +0.8  (cos  -J-  -  1)  *  .894 

Xj  .  I  +  0.8  (cos  -  1)  *  0.60 

6 

Xj  X  I  +  0.8  (cos  in  -  i)  .  0.200 

6 

X  I  +  0.8  (cos  ^  -  1)  *  -  .200 
Xg  X  I  +  0.8  (cos  -  1)  *  -.492 
X^  X  I  +  0.8(cos«-l)  *  -0.6 


(33) 


(34) 


(35) 


Now,  consider  the  quantity 

a"  b 

which  is  typical  of  terms  occurring  in  (13). 
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A"b 


A  q 


m*  0 


100  1 

Now  (.894)  is  less  than  -  so  that  for  the  case  cited  in  (35)  , 

30. 000 


when  n  >  100,  recalling  that  Xq  *  1, 


A^b  i  (Xq)"  b^q°  «  b^q* 


t  AxQ 
Mk 


by  reference  to  (31). 


Finally,  referring  back  to  (13),  for  n^ ,  n^  both  greater  than  100, 


-  T 


+  a"'-‘  t . + 


*  (nj  -  nj) 


>AxQ 


A  X  A 

A  +  A 
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but 


P  AxQ 
M  k 


kAt  Ax  0 
c(Ax)^  Mk 


Q  At 


h  c 


(39) 


>  o  that,  for  Hj  ,  n|  both  >  LOO 

T*"*  -t"'  i  {tz  -tj)  Q 


h  c 


(40) 


which  is  seen  to  correspond  exactly  to  the  analytical  result  of  equation  (5). 

Thus  it  has  been  shown  that  the  difference  between  the  temperature 
distributions  calculated  by  the  explicit  difference  technique  at  any  two 
time  steps  agrees  with  the  analytical  result  after  a  certain  time  (say  corres¬ 
ponding  to  100  steps)  has  elapsed.  This  shows  that  there  is  no  accumulation 
of  truncation  error  despite  the  coarse  space  mesh  regardless  of  the  number 
of  time  steps  computed. 
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3.  PROBLEM  B 


Coniider  next  the  case  where  no  heat  flows  through  either  wall  but 
where  there  is  a  nonzero  temperature  distribution  at  some  time  to. 


The  argument  used  in  this  case  to  prove  that  truncation  error  does 
not  accumulate  aifter  a  large  number  of  time  steps  is  based  simply  on  energy 
conservation.  First,  it  is  known  that  the  solution  to  the  difference  scheme 
(which  is  simply  system  (6)  with  no  Q  term)  approaches  the  flat  or  constant 
temperature  condition  with  increasing  time  and  that  this  uniform  distribution 
is  the  correct  one  for  large  time.  The  only  question  remaining  is  therefore 
that  of  a  possible  constant  vertical  displacement  between  the  exact  and  the 
difference  solutions.  Once  again  round-off  error  is  not  considered. 
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It  ia  this  shift  which  can  be  shown  not  to  accumulate  by  virtue  of 
accumulated  truncation  error. 

Consider  the  energy  content  or  simply  the  integrated  temperature  for 

the  initial  distribution, 
h 

T  (x,  to)  dx  (41) 

This  is  approximated  by  the  trapezoidal-rule  quadrature 

n  “o 

w  T  °  Tw  /»o  no 

T(x,  to)  dx  *  — 2—  + -  +  IT,  +  Tj  + 

2  2 

where  . 

Define  the  scalar  quantity  I*'  by 

l"  *  *  TM,-  +  +  T?  ♦ - +  (45) 
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n+l 

Now  consider  the  expression  for  I  utilizing  the  homogeneous  form 


of  system  (6). 

n+l  n+l  j 

.2^ .  .  (Tr*  ^  Tr ‘  ^ . 

*  T  (to"  +  2  P  (T?  -  To"))  +  J  (tJ,  +  2  P  (t"  (44) 

Jill/  N 

+  >  (Tq.  +  p  (T„.,  +  -  2  T„  )^ 

m=  I  . 


Expanding  and  collecting,  this  becomes 

n+l  I  ,  n 


n+l  I  ,  n  ,_n  _n  ^  in  ,  n  n 

I  =  j  (Tot2?CI^  -To)  +J  (  T^  +  2P  (T,^.,  -  T^) 

+  t"  +  p  +  Tj"  -  2  T?) 


n 

+  Tj  + 


P  (t? 


n 

+  T, 


2  T 


?) 


n  /  n  n 

+  ,  +  pfT^,  ,  +  ,  -  2  T 

M-2  M-3  M-l 


M-2  I 


M-l 

+  ^ 

^M-2 

+  T 

M 

r*.n 

-  2  T 

M 

1  ™,n 

1 

rr.  n 

/'n 

n 

-  T 

2  o 

"I 

T 

M 

+  1  Tj  +  T,  + 

"■"m-J) 


=  r 


(45) 
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Thu*,  neglecting[  round-off,  the  trapezoidal-rule  quadrature  i* 

XI 

"conserved".  That  is,  the  value  of  I  is  the  same  as  the  value  of  I 
for  all  n  greater  than  hq.  Any  error  (aside  from  round-off)  in  the  trape- 
zoidal-rule  approximation  to  I  Tdx  appears  directly  as  a  quadrature 
error  at  initial  time  tg,  and  never  changes.  This  rules  out  the  possibility 
of  an  increasing  displacement  between  true  and  difference-calculated 
temperatures  due  to  accumulated  truncation  errors  for  large  times. 


Stated  differently:  the  analytical  or  exact  solution  conserves  the 


integral 


,  0)  dx  for  all  time  after  Q  vanishes.  The  ditference 


solution  conserves  i  which  is  the  trapezoidal  approximation  to  thir  integral. 
Neglecting  round-off,  if  T(x,  t)  is  the  exact  solution  with  initial  values 


T  (x,  to)  and  if 


["o= r 

Jo 


T  (x,  to)  dx  + 


(46a) 


Then 


where  Co  is  quadrature  error  at  time  to 


n 

I  =1  T(x,  t)  dx  + 


for  all  n  >nr 


(46b) 


since  the  two  integrals  are  equal  and  =  I  ^  . 

Thus  the  accumulated  truncation  error  cannot  cause  the  increasing  divergence 
of  the  exact  and  difference  solutions;  i.  e. ,  differs  from  the  integral 


I  T  dx  at  time  t  by  exactly 

^h 

I  T  dx  at  time  to* 


the  same  amount  that  I  differs  from 
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APPENDIX 


EXACT  SOLUTION  OF  PROBLEM  (A) 

Now  it  can  be  shown  that  the  general  n^^  eigenrow^orresponding  to 
*  !•  +  2  given  by 


[- 


1.  ..n 


*(n)  *  I  i  •  ®n*  2®n»  cos  30^,  cos  40n,  . . .  ,cos(M-I)0n,-|  cos>^ 

(47) 


where  it  should  be  recalled  that  0^  *  nir/M. 

The  expression  (29)  for  p°  is  obviously  a  special  case  of  this  for 

(n  *0,  0n  *0,  >1).  With  some  trigonometric  effort  it  can  be  shown  that 

,th 


the  n  norm  ^(xi)  given  by 

”(n)  •  (‘’(n)  •  ^n))  *  “ 


for  n  >0,  M 
*  M/2  forO<n<M 


(48) 


Having  this  data,  it  is  now  possible  to  write  down  a  complete  solution 
to  problem  (a). 

From  expressions  (12)  and  (27) 

t“  .  (a"'‘  +  a”'^  + . +  A*  +  A  +1  )b 


(a“-‘  +a“-^  + 


.A*  .A  .1)^  b, 


(49) 


M 


'*(m)  +....+  1)  q 


(m) 


m 
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And,  we  have  finally  for  T“  the  result: 


h  c 


I  ,  ,  7  I 

I  n-l  n-2  n-3  .  , 

n  [  -2(X.,  +\.,  +...+  1)  cos  «/M 

}  cos  Ztt/M 

t  cos  3ir/M 


cos  M  ir/M 


+  2^“'^  +  +. . .  +  ^Fcos  2ir/M  1  .  +^r cos  3tr/M 


cos  4ff/M 


cos  2Mit/m1 


cos  6v/m 


cos  3ir/M 


.  ,  ,,M  /n-l  n-2 


X  »  I  +  2  0  (cos  m  -  li 

“KT 


from  (26) 


^o  *  ^ 


n-l  -n-2 
I  +  X  + 
m  m 


I  -  X, 


)*Ti 
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Thus  (54)  may  be  rewritten  in  the  form 


Til  r  1 

2(1 -K,)  cosir /M  cos  2ir/M 

-7; — cos2ir/M  "T”“7 —  cos4ir/M 
(I  -  ^  "^2 


n.  cos  3it/M 
cos  6ir/M 

1-X. 


Icos  Mir/M 


Icos  2M1r/^f 


cos  3Mir 


+(-i) 


M  (1-X„) 


Recalling  that,  for  p  <  ^  »  j  rnj  ^  ^  ni  >  0,  expression  (56) 

may  be  simplified  for  a  time  sufficiently  large  that  and  hence  , 


m  >1,  are  negligible. 


In  large  enough  so 
that  x"  «1 


cosir/M  +  _ L 


2ir/M 


I  I  p(l-cos^)  cos2ir/M  p(l-cos^)  cos4ir/M 


cos  ir 


cos  2  ir 


cos  3ir/M  + . + 


^(1-cos— )  cos  6ir/M 


cos  3  ir 


2P(l-cos^) 
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In  the  above  expression,  only  the  first  term  is  time -dependent. 

This  corresponds  to  (and  agrees  with)  the  first  term  in  (1).  The  remaining 
terms  must  correspond  to  the  second  and  third  terms  in  (1).  A  numerical 
check  indicates  that  this  is  very  closely  satisfied  for  the  particular  case 
where  M  *  6.  In  fact,  the  error  can  be  shown  to  be  given  by  the  difference 
between  exact  integrals  or  weighted  integrals  of  these  terms  and  seven- 
point  trapezoidal-rule  quadrature  approximations  thereto. 
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C.  THERMAL  STRESSES  IN  SPHERICAL  SHgLLS 
OF  ARBITRARY  THICKNESS  UNDER 
ARBITRARY  AXISYMMETRIC  TEMPERATURE  DISTRIBUTIONS 

1.  INTRODUCTION 

The  present  report  describes  an  analytical  method  for  determining 
the  elastic  stresses  developed  in  a  complete  spherical  shell  of  arbitrary 
but  uniform  thickness  under  axisymmetrical  thermal  loading.  A  related 
problem  has  been  solved  by  McDowell  and  Sternberg  (Ref.  1)  with  the 
restriction  that  the  thermal  loading  be  entirely  steady  (i.e. ,  the  tempera¬ 
ture  satisfies  Laplace's  equation).  It  is  the  purpose  of  the  present  report 
to  extend  their  analysis  to  the  case  where  the  impressed  axisymmetric 
thermal  loading  is  an  arbitrary  one  and  may  have  come  from  transient, 
unsteady  conditions.  In  other  words,  the  temperature  distribution  is  not 
restricted  to  be  harmonic  in  the  present  analysis. 

The  results  should  be  useful  in  two  respects.  First,  they  provide 
estimates  of  thermal  stresses  in  the  nose  region  of  blunt-nosed  radomes 
under  transient  thermal  loading.  Second,  they  should  serve  as  numerical 
checks  on  the  solutions  to  radome  thermal  stress  problems  obtained  by 
thin  shell  theory  in  the  nose  region. 
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2.  ANALYSIS 

In  the  absence  of  body  £orcee»  the  general  equations  of  thermoelasticity 
are  given  (Refs*  1*  3  and  4)  as 

Equilibrium:  ^  0  (0 

Stress-Strain:  T^js2|fc  [*‘J  tl-2v 

8Ui 

Strain-Displacement:  Ze^i  >  _ I 

^  8xj 

where  p  is  the  first  Lame^  constant,  v  is  Pois8on*s  ratio,  a  the  coefficient 
of  linear  expansion,  and  xj  a  rectangular  cartesian  coordinate  system. 

In  the  absence  of  surface  tractions,  the  boundary  conditions  are  given  as 

Tjj  1  j  *  0  on  B  ,  (4) 

where  are  the  scalar  components  of  the  outer  unit  normal  to  the 
boundary  surfaces  B. 

The  solution  to  (1),  (2),  (3)  will  be  obtained  in  the  following  manner: 

The  complete  solution  [s]  will  be  subdivided  into  two  solutions, [s^  and[s*] 
such  that 

[s]  •  [s<5  +  [s<3  .  (5) 

where  js°]  represents  any  particular  solution  to  the  general  equations  (1), 
(2) ,  (3) ,  and  [S*]  is  the  solution  to  the  so-called"residual  problem"  which 
satisfies  (1),  (2) ,  (3)  with  T  *0  and  annuls  the  surface  tractions  which 
(S^  gives  rise  to  on  the  boundary  B. 
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It  may  be  noted  that  a  particular  solution  to  (1),  (2) ,  (3)  may  be 
found  immediately  by  setting 

'^i  *  ^  '  »  (6) 

‘"‘i 

with  =  -i-liiSL  T  •  ko  T.  (7) 

I  -  V 

The  temperature  field  T  in  this  problem  need  not  be  steady-state 
(i.e.  T  0  )  and  therefore  4  need  not  be  biharmonic. 

This  known  axisymmetric  temperature  distribution  T  (r,  cos  0) 
may  be  put  in  the  form 

CD 

T  (r.  cos  e)  =  ^  F„(r)  Pn(cose)  .  (8) 

"n*  0 


where  r,  0  are  the  radial  distance  from  the  origin  and  latitude  measured 
from  the  Xj  axis,  respectively.  (The  cyclic  longitude  coordinate  may  be 
denoted  by  y.)  The  functions  Fn(r)  are  given  by 


F„(r)  = 


Zn-t-l 


T(r,  cos  0)  P„(co8  0)  d(cos  0), 


(9) 


The  particular  solution  to  the  Poisson  equation  (7)  may  be  obtained  by 
first  substituting  functions  of  the  type 

(r,  cos  0)  *  fn(r)  •  Pn(cos  0)  (10) 

into  (6)  where  in  spherical  coordinates  becomes,  for  axisymmetry, 

._£L  .-L J.  .  _L  JL  (u) 

dr*  r  0r  r*  00  r*  80*  * 
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There  rerults  a  differential  equation  on 

£„"(r)  +  A  £„'  (r)  .  yr)  -  k„  F„<r)  . 


for  which  a  particular  solution  may  be  given  in  terms  of  a  Green's 


function  as 


yr)  = 

2n+l 


jn-l 


n+2  n  ,-r.+  l 
&  -  r  ; 


where  r|  is  the  radius  of  the  outer  shell.  (The  radius  of  the  inner  shell 


will  subsequently  be  denoted  by  r^.  .)  And  hence. 


4^(r,  coe0)  = 


-n-l  , n+2 
r  4  -  r 


I 

n^  n+l  ^  P^(cos0) 


is  the  particular  solution  to  (6)  that  is  sought. 
For  notational  purposes  define 


r"  r„(t)  d;  . 


‘  r  ’  • 


A  r  o’’ 

and  denote  cos  0  by  p  and  sin0  by  p.  Then  may  be  written  as 

('^n  -  ^nir))  P^ip)  ,  (17) 

radial  deformation;  u°  = =  -r  j^(n+l)  Gn(r;  +  nHj^(r^  P^tp)  ,  (18) 

o 

meridional  deformation:  w®  =  =  -  P  r  j^G„(r)  -  H^(r)J  P;,(p),  (19) 
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'rr 


YY 


•  jfn+l)  (n+2)  a„(r)  -  n(n-l)  H„(r)  +  F„(r)1  P„(p).  (20) 


^  i  •  -  j(n+l)  G„(r)  +nH„(r^  P„(p)  -p[o„(r)  .H„(r^P;(p)  (2« 

^  =  -j^(n+4Ci,(r)  +nH„(rijp„(p)  -p^(r)  -  H„(r)jp^(p) 

+  d-p")  jG„(r)  -  H„(r)1p;;(p),  (22) 

I  /l  8u°  w°  9w_\  .  r  ”1 

ll;;  T^"  ■  ^^-877' P  ^  H„(r)Jp„'(p).  (23) 


rO 


Tr^  -  2P 


(n+l)  (n+2)  Gn(r)  -  n(n- 


T  *  -  2  p 
YYn 


-1)  H„(r)J  Pn(p)  . 

|(n+l)  Gn(r)  +  nHi(r)  +koFn(r^  P„(p)  +  pjq,(r)  -  H„(r)Jp^(p) 

p>  ^ 

Ijn+l)  Gn(r)  +  nf^(r)  +  koFn(r)Jpn(p)  +  pjciJir)  -  I^(r)j  P,i(p) 

-  (l-p')[G„(r)  -  H„(r)j  p;;(p)| 


(24) 

K  (25) 


r0, 


*2  p  ^  n 


(n+2)  Gjj(r)  +  ( 


L 


n-l)  H„(r)J  Pn' 


(p)  . 


.  (26) 


(27) 


where  the  primes  denote  differentiation  with  x'espect  to  p. 


The  next  step  is  to  obtain  the  solution  to  the  "residual  problem" 
(where  T  *  0)  and  whose  boundary  conditions  are 

T*j,^(rk,  p)  s  -  2p  J^(n+l)  (n+2)  G„(rk)  -  n(n-l)  Hn(rk)J  Pn(p)  , 

'^e^(rk,  p)  »  -  2p  ?  J(n+2)  G^(rk)  +  (n-l)  Hj^(  rk)J  Pn(p)  , 


(28) 

(29) 
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where  k«  0,  1.  In  view  of  (15)  and  (16),  on  the  outer  surface 


(r|.p)  *  0  , 


while  on  the  inner  surface  . 


Tjr^  (ro.  P)  «  -  2  ,1  |(n-rt)  (n+2)  G„(ro)  -  n(n-l)  Hn(xb)J  Pn(p) .  (32) 

^r0n  $|(n+2)  G„  (rc)  +  (n-l)  Hn(ro)J  Pn(p)  •  (33) 

The  solution  of  the  residual  problem  of  Elasticity  Theory  can  be  traced 
through  Refs.  6,  5,  7,  2  ,  and  finally  I  ,  where  McDowell  and  Sternberg 
have  presented  an  explicit  solution  in  the  form  of  a  linear  combination  of 
four  separate  solution  fields.  That  is. 

[®nj  *  •n[^i:^  ^n  nj^-n-lj  **0  |®-n-lj  » 

where  for  : 

2 

nr.  -  n 

2jjiu*— f—  (rj)  P„(p)  ,  (35) 

2  p  w  X  -  P  (f^)"  P^(p)  ,  (36) 

T„.  n(n.l)  (^)*(i)"  P„(P)  .  (37) 


- 


Pn'Pi  -P<f,)"  Pn’trt 
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and  for 


["'  ‘ri"  Pn(l»  -  P  (f.>"  P,;(rtl. 

(39) 

Trt  •  -  (-«  (^'  $  (-^t”  P„’  (P)  . 

(40) 

[Bn]  : 

2  p  u  *  -  (n+l)  (n-2  +  4v)  r  (■^’^  P^(p)  , 

(41) 

2  p  w  =  (n+5  -  4v)  r  p  (^“  P^(p)  , 

(42) 

Trr  =  -  (n+l)  [(n+l)  (n.2)  -  2  v J  (•^)"  P„(p)  , 

(43) 

*  -  (n+l)  [n-2  -  2  (2n+l)  vj  (-f^"  P„(p)  +  (n+5-4  v)  p(^"^’(p)  ,  (44) 

Tgg  *  (n+l)  (n*  +  4n  +  2  +  2v)  Pn(p)  -  (n+5  -4v)  p  {^“P^(p)  (45) 


Tre  *  (n*  +  2n  -  I  +  2v)  p  (-fj)”  Pn(p)  • 


The  eolutione  are  obtained  by  replacing  n  by  -n-l 

wherever  it  appear*  in  equation*  (35)  through  (46).  Hence  for 


n-ll  noting  that  P_n-l(p)  *fVx(p)  ): 

-*  y  2  .1 

2  pu  «  -  (n+l)  -p  (t*)  ^  (p)  . 

(47) 

■»  A  /£l%n+l  Tji/  » 

2pw  s  -  p  ("T)  Pn(p)  , 

(48) 

Tyy  .  (n+l)(n+2)  (.^)*  (^)"‘^^  P„(p)  , 

(49) 
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.  -  (^)*  j^(n+«  P^(p)  +  p  P^(p)J  ,  (50) 

Tge  =  -  +  2n+l)  P„(p)  -  p  P^(p)J  ,  (51) 


T"  n  * 
rO 

(n+2)  (^"^  p 

PntP)  • 

(52) 

for 

-n-l 

■ 

2  pu  = 

-  n(n+3  -4v)  r 

Pn<Pl  . 

(53) 

2  pw  s 

-  (n- 

■4+  4v)  rp(J^)*"‘^^P^(p)  . 

(54) 

Trt' 

n  ^n(n+3)  -  2sJ  (■^)"'*^* 

Pi.(P>  . 

(55) 

* 

-  n 

l^n+3  -  2  (2n+l)  vj 

(-^)  Pn^P^ 

-  (n-4  +  4v)  p  (■^)”^^P^(p),  (56) 

■"! 

n*  -2n-  l  +  2vj  (Y^”^^Pn(p)  +  (n 

-4  +  4v)  pf^"'^^P„'(p)  ,  (57) 

<n' 

-  2  +  2.)  5 

Pn'{p)  . 

(58) 

The  coefficients  a^,  b^,  Cjj,  are  determined  from  the  boundary 
conditions  on  T?r  and  .  Combining  (37),  (43),  (49),  (55)  and  equating 
to  (30)  and  then  (32);  and  then  combining  (40),  (46),  (52),  (58)  and  similar ly 
equating  to  (31)  and  then  (33)  will  yield  n  sets  of  four  simultaneous  equations 
for  the  determination  of  a„,  bj^,  c^^,  dj^. 
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These  equations  are 

n(n-l)  an  -  (n+l)  ^(n+l)  (n-2)  -  bn+  (n+l)(n+2)  €„  + n|n* +3n -2vJ  =  0  , 

(59) 

-(n*  1)  an  +  (n*  +  2n  -I  +  2v)  bn  +  (n+2)  Cn  +  (n*  -  2  +  2v)  dn  *  0  ,  (60) 


n(n-l)  ^an“  (j^+O  j^(>^+l)  (n-2)  -2v  p"bn  +  (n+l)  (n+2)  p*^"^Cn 


(61) 


+  n  n*+ 3n-2vj  p“"^dn*- ZpiRn+l)  (n+2)C:h(ro)  -  n(n-l)  Hn(rJ]  , 


-(n-l)  p"  ^an  +  (n*+2n-l +2v)  p”bn  +  (n+2)  p*“"^  Cn  +  (n*  -  2  +  2v)  p”"^dn 

r  -j  (62) 

*  -  2pL  j^(n+2)  G^(ro)  +  (n-l)  Hn  (ro)J  , 

where  p  *  <  I  .  Note  that  the  terms  in  |snj  containing 

diQ,  a| ,  dg  will  not  contribute  to  the  stresses  and  therefore  represent 

rigid  body  motions  and  will  be  indeterminate. 
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APPENDIX  I 


SIMPLIFICATION  OF  AND 


The  dependence  of  e^.  (£q«  22)  and  (Eq.  26)  upon 

eOn  w^n 

Pq(p)  may  be  removed  in  the  following  manner: 

The  differential  equation  satiefied  by  Pn(p)  i*  well  known  to  be: 


(I  -  pM  p^'*- (n*  +n)  P„  +2pp;i 


(A-1) 


And  therefore  the  original  expressions  can  just  as  well  be  written  as 


o 

*een‘ 


r<n+l)‘  ■  rf  P„  (p)  +  p  [g„  -  H„j  P^(p)  . 


(A -2) 


T”  « -  2  p 

00n 


1 

(n+l)*G„.n"Hn  +  koF„Jp„(p)  -pjG„-H„jp^{p)l  (A-^) 
.  J 
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APPENDIX  n  -  TEST  CASES 


C»»e  1  ;  T  (r,  p)  *  Tq  (constant  throughout  shell) 

This  case  should  give  a  trivial  stress  distribution. 
By  inspection  FQ(r)  *To  ;  Fn(r)  *0  ,  n^O. 


Oo('l  •  r”  To  dt  .  I  ko  To  -  ij 

Ho(r)  •  koj  r-'  ;  To  dl  .  j  k„  T„  -  ij 

The  particular  solution  gives  >  2p  |2  CV>(>^)J  *  “f"  J^oTojt^)  -  ij 

and  ^80  '  "  J'o^ot')]  *  --X  *^YY 


and  T 


'^ra'  0 


For  the  residual  problem  the  simultaneous  equations  reduce  to 


2(l+v)  bo  +  2co  =  0 


2(1+ v)  bg  +  2  Co  p"*  =  -  2p  1^2  Go(ro)] 


which  yield 


Cq  =  "  j  p  J^o^o 


^o  =  -f" 


2  p  hoT  < 


(l+v) 


T*.  =  2(l+v)  b„  t  2  {^’  '0  =  4^  koTo  fl  -  l4)’j 
=  2(l+v)  lb  -  (4)’  Co  •  koTo  [2  +  J 


^re=  » 
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and  hence  ^Vr  *  *^rr  '^fr 

Tbq  =  '’bo  '’^00 

'^yy  ‘  '^99  ~  ° 

=  0 


which  is  the  correct  resuLt. 


»  0 
=  0 


Case  2  .  T(r,  p)  =  Tq  ("^  cos  0  =  Tq  (-^  Px(p) 

This  (linear  temperature  distribution)  case  should  also  give 
trivial  stresses 

By  inspection  Fi(r)  =  To(-^)  ;  0  (n^  1) 

r, 

ko  r  r*4  .4 


G,(r)  = 


Hi(r)  = 


/  4 1 


The  particular  solutions  are 

T°,  =2|i‘  6  Gx(r)  cos  0  * koTj, 


cos  8 


T®q  =  6p  Gi(r)  sin  0  =-f  .ii>  0 

T00  =  -  2jij3Gi<r)  -  ko  F,(r)j  cos  0  =-^  koToji^^+4t^j  cos  0 


T  -  T 

'yy  -  '00 


ilkDD  TR  59-22 


764 


The  simultaneous  equations  become 

2(2+2v)  b|  +  6  Cj  +  (4-2y)  dj  *  0 
(2+2v)  bj  +  3  Cl  +  (-l-t-Zv)  di  *  0 
2(2+2v)  p  bi  +  6p-^  Cl  +  (4-2v)  p-^  di  =  -2  |i-  6Gi{r^ 

(2+2v)  p  bi  +  3p~^  Cl  +  (-I'fZv)  p"^  di  *  -  2  ji*  3  Gi(ro) 

from  which  di  *0;  ci  *  -  JjL  k^To  ;  Gi  «  **^0^0 

^  5(l+v) 

and  T*r  =  [4(1+ v)  bi  +  6(^'*  cj  cos  9  *  ^  »^oTo[<|j)  -  (^M^cos  0 
T*e  =  [2(l+v)  (^)  bi  +  3  Cl  J  sin  9  =  -|-  ^koT^  |i|j)  -  (|lj^sin0 

^60  “  ‘j®^^'*''')^W-3^^)^Cijco80  *-|-pkoToj4^+(^^os0 

Hence  T„  x  x  0 

^r0  *  •^'^9=0 

'^99  =  =  '^Qe  *  '^9*9  =  0 

which  is  the  correct  result. 


/AID  TR  59-2? 


765 


Case  3  .  T(r»p)  *  ,0-)  (Surfaces  at  different  uniform  temp.) 

This  is  a  purely  radial  temperature  distribution  for  which  the 
stress  distributions  are  easily  calculable  (References  3  and  4  ). 

By  inspection  Fo(r)  «  T9T**,  (  £i-  -  l  )  ;  F„(r)  *  0  n^O 


Go<') 


HJr)  =  ko 


4 


(f  ■  y  ^  4] 


The  particular  solution  gives  *2p  j2Go(r)J> 


ZpkoTorori  ^  +  2J] 

rj-r©  [3 


*  0 


For  the  residual  problem  the  simultaneous  equations  become 

2(l+v)  bo+  2  Cq  *  0  (1) 


Oq  +  (-l+2v)  bQ  +  2  c©  +  2(-l+v)  d©  »  0 
2(l+v)b©+  2p'’  Co  s  -2p|2Go(ib)J 


(2) 


(3) 


a^p”*  +  (-l+2v)  b©  +  2  p"’  Cq  +  2(-l+v)  p  'd©* 
Only  (1)  and  (3)  need  be  considered;  they  yield 


2p  r2G©(ti>).H,^J 

(4) 


Cq* 


p  k©T©  r©  £  _  j  _i  2 


-cq 

T+T 
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Hence 


T^.O 


V#  •  T^=«‘+->  'h-^'  Co=^!i^sZ£ 


And  adding  and  simplifying: 
T  /»  =  0  and, 


(l-p-^(r| 


l-v  rj-r^  r* 


^»9  * 


aET^  r^ri 


0*^1  I  ,  r.*  r  * 

(ri+  r,ro+  r*)  ^ 

I  -  r©  ir  2  r* 

^  wa 


This  solution  checks  with  the  stress  distribution  calculated  using 
the  simple  schemes  available  (Ref.  3  and4)for  purely  radial  temperature 
distributions. 
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CHAPTER  IV  -  PART  D 

DIGITAL  PROGRAMS  FOR  THE  STRUCTURAL  ANALYSES 
OF  HOMOGENEOUS  CONICAL  AND  OGIVAL  SHELLS 
UNDER  MECHANICAL  AND  THERMAL  LOADING 


By  Daniel  E.  Magnus 
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SYMBOLS 


a 


dm 

h 

i 

m 

n 

P 

W,  V,  u 

y.  z 
A 

Cm 

F 

J 

Ki 

Kz 

K4 

K5 

K6 

K? 


fraction  of  meridian  at  which  truncation  occurs 
metric  tensor 

coefficient  =  0  for  m  =  0  ;  =  L  for  m  >  0 

shell  thickness 

spring  constants  i  =  1,  7 

meridian  length,  extended  to  apex  (L) 

harmonic  number 

number  of  spaces  along  generator 

pressure  (Ib./L*  ) 

physical  deformation  components  (See  Figure  I  ) 
coordinates  (See  Figure  1) 
equivalent  to  /  and  y^,  respectively 
matrix  associated  with  quadratic  form 
coefficient  =  I  for  m  =0,  =  l/2  for  m  >  0 
free  energy 

Jacobian  or  volume  element 

edge  moment  spring  constant  at  base;  k^/c 

edge  shear  spring  constant  at  base;  kj  /c 

transverse  displacement  spring  constant  at  truncation;  k,  /c  i 

bending  spring  constant  at  truncation;  k^/el* 

edge  moment  spring  constant  at  truncation;  ks/cl^ 

edge  shear  spring  constant  at  truncation;  k^/t 
axial  stretching  spring  constant  at  truncation;  k^/tf 
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Ri ,  Rj  principal  radii  of  curvature;  also  R  s  Rj 

T  temperature;  reference  temperature 

W  virtual  work 

a  coefficient  of  linear  expansion 

P  Local  meridian  -  curve  slope 

Po  P  at  apex 

c  Young's  modulus;  =  reference  modulus 

V  Poisson's  ratio 

0  Longitude  angle  = 

X  cone  semi -angle 

stress  tensor 

T]..  strain  tensor 

'ij 

physical  stress  in  0  direction;  k  =1,2,3  for  z  =  -h/2,  0,  h/2 
Cj?  physical  stress  in  y*  direction;  k  =  I,  2,  3  for  z  =  -  h/2,  0,  h/2 

physical  shear  stress;  k  =  I,  2,  3  for  z  =  -  h/2,  0,  h/2 
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D.  DIGITAL  PROGRAMS  FOR  THE  STRUCTURAL  ANALYSES  OF 
HOMOGENEOUS  CONICAL  AND  OGIVAL  SHELLS  UNDER  MECHANICAL 

AND  THERMAL  LOADING 

1.  INTRODUCTION 

In  this  part,  the  digital  programs  for  the  structural  analysis  of  the 
homogeneous  conical  and  ogival  shell  under  mechanical  and  thermal  loading  are 
described.  Both  programs  are  based  upon  the  same  analytical  approach,  but 
as  will  be  seen,the  two  programs  use  different  numerical  procedures  in 
formulating  the  energy  integrals.  In  general,  the  numerical  techniques  used 
in  the  cone  program  will  be  explained  in  detail  and  then  supplemented  by  any 
variations  used  in  the  ogive  program. 

The  cone  program  is  designed  for  the  structural  analysis  of  conic  shells 
with  closed  or  trimcated  nose  conditions.  By  merely  changing  the  problem 
inputs  and  machine  setting,  the  analysis  can  be  changed  from  one  condition  to 
another.  In  this  regard,  it  should  be  emphasized  that  the  closed  cone  is  not 
analyzed  as  a  specialization  of  the  boundary  condition  used  for  the  truncated 
cone.  Rather,  independent  requirements  are  imposed  for  the  closed  cone; 
a  factor  which  adds  considerably  to  the  internal  complexity  of  the  program. 

The  ogive  program  is  for  the  closed  nose  shell  only. 

Both  programs  are  designed  for  the  analysis  of  any  general  loading  where 
temperature  and/or  pressure  are  involved.  In  order  to  analyze  a  shell  with  a 
generic  loading,  the  load  must  first  be  expressed  as  a  Fourier  series  in  the 
0  direction.  (See  Figure  L  for  nomenclature. )  Each  of  the  Fourier 
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harmonic*  must  be  separately  analysed  by  the  digital  program,  and  since  the 
analysis  is  based  upon  linear  theory,  the  final  deflections  and  stresses  can  be 
found  by  superimposing  the  resiilts  from  each  harmonic.  Taking  accoimt  of  the 
lack  of  precision  of  most  aerodynamic  and  thermal  input  data  derived  from 
unsymmetric  loading  conditions,  it  is  probable  that  no  more  than  two 
harmonics  will  be  required  to  adequately  describe  a  realistic  loading. 

The  analysis  of  the  shell  geometries  is  based  upon  a  variational  principle 
involving  the  "free  energy",  F,  of  the  system  and  the  external  work,  W,  of  the 
mechanical  loading.  The  concepts  and  analytical  details  of  the  analysis 
appear  in  Part  A.  It  is  found  that  the  nth  harmonic  component  of  the  free 
energy  can  be  written  in  the  following  form,  v^ich  is  suitable  for  both  the  cone 


and  ogive  geometries: 


2irC„^l*- jj  f  dy»dy|(4_) 

o  i  ^  J 


/«  \  2o(T-Tq)  ^ 


where 


As  will  be  shown  in  Section  Z,  the  numerical  method  of  evaluating  the 
energy  of  the  cone  is  not  the  same  as  the  method  for  the  ogive.  For  the  cone, 
the  integration  through  the  thickness  (y  )  is  carried  out  analytically, and  the 
integration  along  the  generator  is  done  numerically.  For  the  ogive  both  the 
integration  through  the  thickness  and  along  the  generator  are  done  by 
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numerical  methods;  the  integration  through  the  thickness  being  carried  out 
using  a  Gauss  quadrature.  In  both  programs  the  integral  along  the  generator 
is  expressed  in  finite  difference  form  as  a  fimction  of  the  w,  v,  and  u  dis¬ 
placements  of  the  mesh  points.  Then,  by  taking  the  variation  of  the  energy 
expression,  a  set  of  linearly  simultaneous  equations  for  w,  v,  and  u  are 
obtained  for  the  n  stations  along  the  generator.  The  system  of  equations  have 
the  very  important  property  of  being  symmetric  positive  definite,  which  is  a 
direct  consequence  of  the  energy  formulation.  As  will  be  shown  in  Section  4, 
a  positive  definite  matrix  is  a  sufficient  condition  for  the  block  tri-diagonal 
method  of  solution  of  the  equations. 

By  way  of  contrast,  an  alternate  method  of  obtaining  the  simultaneous 
equations  will  be  described,  and  the  inherent  disadvantage  noted.  The 
variation  of  the  energy  expression,  equation  (1),  could  be  taken  in  analytical 
form  thus  resulting  in  the  Euler  eqtiations  and  the  associated  boundary 
conditions.  These  differential  equations  would  then  have  to  be  expressed  in 
finite  difference  form.  Unfortunately,  the  final  simultaneous  algebraic 
equations  would  not  necessarily  be  symmetric, and  the  coefficient  matrix  could 
appear  as  either  positive  definite  or  indefinite  depending  on  the  method  of 
differencing.  Furthermore,  the  solution  of  the  equations  would  be  considerably 
more  involved.  This  disadvantage  is  particularly  significant  since  the  matrix 
under  consideration  can  be  of  the  order  450. 

As  yet,  no  mention  has  been  made  of  the  boxmdary  conditions  to  be 
imposed.  In  general,  the  boundary  conditions  are  expressed  in  terms  of 
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springs,  ^ichovary  in  magnitude  from  zero  to  infinity.  A  zero  spring  implies 
the  shell  is  free  of  a  particular  constraint.  For  finite  values  of  the  springs, 
the  energy  contribution  is  computed  and  added  directly  to  the  energy  matrix. 
An  infinite  spring  represents  a  rigid  constraint,  and  in  order  to  avoid 
numerical  difficulties  in  the  limit,  special  procedures  are  employed. 

In  Figure  2  ,  the  various  types  of  spring  configurations  are  shown. 

The  base  conditions  for  the  cone  or  ogive  are  described  by  springs  and  K2. 
The  other  5  springs  shown  in  Figure  2  are  applicable  only  to  the  truncated 
cone.  With  these  springs,  it  is  possible  to  impose  any  boundary  condition  at 
the  truncation  including  the  conditions  of  the  unified  radome  with  a  nose  boom. 
Regarding  the  closed-nose  cone  or  ogive,  independent  closure  conditions  are 
imposed  depending  on  the  harmonic  n\imber.  The  methods  of  introducing  the 
spring  and  closure  conditions  are  described  in  Section  3. 

Upon  solving  the  system  of  simultaneous  equations,  the  displacements 
w,  v,  and  u  are  found  for  n  points  along  the  generator  of  the  shell.  These 
displacements  and  the  first  and  second  derivatives  with  respect  to  y^  are  the 
necessary  variables  for  determining  the  state  of  stress,  which  is  discussed 
in  Section  5  below. 

The  prograims  have  been  written  for  the  IBM  704  digital  Computer  with 
a  core  capacity  of  8K  and  at  least  2  magnetic  tapes  for  the  cone  program  and 
3  magnetic  tapes  for  the  ogive  program.  In  addition,  the  ogive  program 
reqmres  2  physical  drums.  The  programs  can  be  read  into  the  computer 

I 

using  either  binary  cards  or  binary  magnetic  tape.  The  latter  method  of 
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operation  is,  of  course,  desirable  for  economy  reasons. 

When  using  the  maximum  number  of  spaces  (n  =  72)  along  the  generator 
of  the  cone,  the  total  running  time  for  the  zeroth  harmonic  is  less  than  3 
minutes;  for  the  first  harmonic  with  216  simultaneous  equations,  the  program 
reqmres  approximately  5  minutes.  Approximately  75  seconds  of  this  time 
are  required  to  solve  the  216  simultaneous  equations. 

In  the  ogive  program,  the  maximum  number  of  points  along  the  generator 
was  increased  to  151,  thereby  requiring  the  solution  of  453  simultaneous 
equations  for  the  first  harmonic.  As  yet.  a  complete  problem  has  not  been 
run  using  the  ogive  program,  because  part  of  the  stress  calculation  is  still 
being  modified.  Hence,  an  accurate  timing  of  the  program  is  not  available, 
but  it  is  estimated  that  the  total  runxiing  time  will  be  10  minutes. 

Numerous  test  problems  have  been  analyzed  using  the  cone  program,  and 
5  of  these  problems  are  discussed  in  Section  6  of  this  report.  Most  of  the 
problems  are  for  the  zeroth  harmonic  since  these  problems  are  more 
amenable  to  theoretical  analysis,  but  some  discussion  is  included  for  a  cone 
with  a  first  harmonic  thermal  loading.  The  ogive  program  has  been  used  to 
compute  the  displacements  w  and  v  for  the  zeroth  harmonic  with  constant 
pressure  loading.  These  results  are  presented  in  graj^ical  form  and  are 
compared  with  theoretical  results  based  upon  membrane  theory. 
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FORMULATION  OF  THE  ENERGY  MATRIX  AND  LOADING  VECTOR 


For  the  cone  and  ogive,  e^ich  are  axi symmetric  shells  with  a^^  °  °  ^ 

J 

for  i  ^  J.  the  free  energy  part  of  equation  (1)  can  be  written  in  the  follewixig 
form: 


-  2a(l4-v)  (T-To) 


a^^0,a22^ 


'll 


'12 


L^22j  J 
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aiyi  0  valla22 
0  2a^a22  o 
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^12 

^122 


(2) 


Equation  (2)  is  in  dimensionless  form  with  all  physical  dimensions  scaled  by 
the  length  of  the  generator  from  the  apex  of  the  shell  to  the  base  mounting. 

The  above  energy  integral  can  be  numerically  evaluated  in  finite 
difference  form  using  two  different  methods.  One  approach  is  to  e3q»ajui 
analytically  the  in  terms  of  the  displacements  w.  v  and  u  and  perform 
the  indicated  matrix  multiplications  with  the  metrics.  The  integrand  of  the 
energy  becomes  a  complicated  expression  involving  various  powers  of  the 
displacements  and  their  derivatives.  These  must  be  integrated  in  the  y  and 

3 

y  directions.  The  integration  through  the  thickness  can  be  carried  out  by 

2 

analytic  procedures  leaving  only  a  single  integral  in  die  y  direction.  The 
integrand  of  the  final  expression  is  expressed  in  finite  difference  form  and  the 
integration  replaced  by  a  summation.  This  approach  for  the  evaluation  of  ^e 
free  energy  expression  is  used  in  the  cone  program  vdiere  the  final  analytical 
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expressions  before  introducing  the  rumination  procedure  are  not  too  com¬ 
plicated. 

2 

For  the  ogive,  the  expression  to  be  integrated  in  the  y  direction  was 
found  to  be  complicated  and  would  require  extensive  programming  effort. 

By  using  a  slightly  different  approach,  a  simpler  programming  procedure  was 
found.  This  method  consists  of  expressing  the  in  finite  difference  form 
and  introducing  them  directly  into  equation  (2)  thereby  eliminating  the  analytic 
steps  of  the  procedure  used  on  the  cone.  The  matrix  multiplication  is  done 
numerically  and  using  a  Gauss  quadrature,  the  matrix  is  numerically 

2 

integrated  through  the  thickness.  Finally,  the  matrix  is  integrated  in  the  y 
direction.  This  procedure  is  specifically  applied  to  the  quadratic  form  in 
equation  (2);  the  linear  form  involving  the  temperature  terms  is  not  complicated 
and  can  be  integrated  analytically  through  the  thickness  as  was  done  with  the 
cone. 

After  the  numerical  integration  by  either  of  the  above  methods,  the 
free  energy  and  external  work  term  can  be  written  as: 

1/2  X' AX  =  (W-T)X  (3) 

where  X  is  the  displacement  vector  with  components  w.  v.  and  u  at  the  various 
mesh  points  along  the  generator.  XAX  is  the  quadratic  form  of  the  energy 
expression  and  (W-T)X  is  the  linear  form  of  the  work  term  W  combined  with 
the  linear  form  of  the  temperature  terms  T  appearing  in  the  free  energy 
I  expression.  By  taking  the  variation  of  eqviation  (3).  the  following  equation  is 
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obtained  for  the  dieplacementa: 

AX  =  W-T  (4) 

'Where  the  matrix  A  is  of  S'ymmetric  positive  definite  form.  For  simplicity 
in  sol'ving  the  simultaneous  equations,  the  rows  and  columns  of  matrix  A  are 
arranged  so  that  all  non- zero  elements  are  contained  in  15  diagonal  rows. 
With  this  arrangement  of  the  matrix,  the  order  of  the  components  of  the  dis¬ 
placement  v  ector  is: 


^n+1 

’^n+l 

'‘n+l 


\diere  n+1  is  the  total  nxunber  of  points  along  the  generator.  The  arrangement 
of  the  X  vector  is  the  same  in  the  cone  and  ogive  program  except  for  fictitious 
elements  which  are  placed  before  w^  and  after  \^^in  the  ogive  program.  These 
elements  will  be  discussed  in  the  latter  part  of  this  section, 
a.  Energy  Expression  and  Grid  Spacing  for  the  Cone 

Omitting  the  rather  lengthy  algebraic  details,  the  free  energy  for 


the  cone  can  be  written: 
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Ulcere  A  =  y  lin  X  -  zco*  X  *  /"in” 


In  many  shell  investigations,  various  approximations  are  introduced 
into  equation  (6)  to  reduce  the  equation  to  a  simpler  form.  For  example. 
Love's  first  approximation  (see  Reference  1)  drops  the  a  in  the  Aterms. 


WADD  ra  59-?2 


782 


Still  other  methods  expand  the  A~^  terms  as  an  infinite  power  series  and  from 
order  of  magnitude  studies,  only  certain  terms  in  the  series  are  retained. 

For  the  digital  program,  however,  there  are  no  reasons  for  introducing 
approximations  into  equation  (6).  As  will  be  seen,  the  terms  at  worst 
introduce  logarithmic  terms  upon  integration  through  the  thickness,  and  these 
terms  are  certainly  not  difficult  to  evaluate  in  the  digital  program.  Hence, 
equation  (6)  is  integrated  assuming  a  uniform  wall  thickness  and  upon 
rearrangement  becomes: 
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where  the  logarithmic  terma  are 
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e  =  — = -  in  ■  .y. 
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and  the  temperature  terms  are  in  the  integral  form 


The  quantities  T^^^and  T^jwhich  can  be  the  results  of  a  heat  condition 

analysis,  are  the  necessary  thermal  inputs  to  the  digital  program  and  must 
be  specified  at  the  mesh  points  along  the  generator.  Since  the  coefficient  of 
linear  expansion  a  appears  imder  the  integral  sign  in  equation  (9),  the  cone 
program  can  be  used  for  constant  or  temperature  dependent  a  •  Regarding 
the  elastic  material  properties.  Young's  modulus  c  has  been  assumed 
constant  and  factored  out  of  the  integral  in  equation  (7);  it  is  not  a  required 
input  to  the  program  but  is  used  to  scale  the  pressure  inputs  as  shown  in 
equation  (1).  Poisson's  ratio  v,  which  is  assvimed  constant,  appears  in 
equation  (7)  and  is  a  required  input  to  the  program.  The  other  program  inputs 
are:  a,  the  distance  from  the  apex  to  the  location  at  the  truncation: 

X  ,  the  semi-angle;  h,  the  thickness;  and  n,  the  number  of  spaces  along  the 
generator  between  a  and  the  base  mounting.  The  number  of  mesh  spaces  n  can 
be  varied  from  12  to  72  in  increments  of  12. 

As  shown  in  Figure  3  ,  the  program  does  not  assign  a  uniform  mesh 
to  the  cone;  rather,  a  uniformly  fine  mesh  is  used  near  the  edges,  and  a 
coarse  mesh  is  used  in  the  middle.  The  two  edge  areas  with  closely  spaced 
points  will  be  called  the  boxmdary  region  awhile  the  area  between  the  boundary 


region  it  celled  the  membrane  region.  The  program  atiignt  n/3  spaces  to 
the  membrane  region  and  the  remaining  2n/3  spaces  are  divided  equally 
between  the  two  boundary  regions.  Further,  the  program  assumes  the 
distance,  6,  between  stations  in  the  membrane  region  is  4  timet  at  large  as 
the  spacing  in  the  boundary  region.  This  is  sufficient  information  to  com¬ 
pletely  establish  the  mesh  points  along  the  generator. 

To  express  equation  (7)  in  finite  difference  form,  the  formulas  shown 
in  Appendix  A  were  used.  Note  that  special  formulas  were  required  at  the 
junction  points  between  the  membrane  and  boundary  regions.  For  the  matrix 
of  the  quadratic  form,  the  difference  formulas  are  applied  to  the  terms 
multiplied  by  h^/l2  in  equation  (7)  and  integrated  over  the  mesh.  This 
results  in  an  n  x  n  matrix,  which  is  added  to  a  null  matrix  Aj^.  Then  the 
terms  multiplied  by  f,  e,  and  h  are  evaluated  in  the  same  manner  and  are 
added  directly  to  the  matrix  A^  as  the  computation  proceeds. 

b.  Energy  Expression  and  Grid  Spacing  for  the  Ogive. 

For  the  ogive,  the  quadratic  portion  of  the  free  energy  expression 
(equation  2)  is  numerically  integrated  in  both  the  y^  and  y^  directions.  To 
avoid  some  of  the  differencing  problems  which  were  manifest  in  the  cone,  the 
terms  in  the  positive  definite  symmetric  matrix  differenced 

by  first  using  a  two-point  forward  first  derivative  and  a  three-point  central 
second  derivative.  The  expression  is  then  reformulated  by  replacing 
all  first  derivatives  by  a  two  -  point  backward 
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first  derivative.  For  the  uniform  mesh  the  resiilting  quadratic  forms  are  then 


combined  in  the  following  manner: 


’>'Tn|  =  i[v 


backward 


(b-l) 


where  T  ■  T 


afiy6 


The  energy  term  is  then  numerically  integrated  in  the  y  direction 

2 

by  uiing  a  five-point  Gauss  qiiadrature}  that  is,  at  any  given  station  y^  , 


where  the  are  h/2  times  the  normalised  five-point  Gauss  quadrature 
coordinates,  and  the  a^  are  h/2  times  the  normalised  five-point  Gauss 
qiiadrature  weighting  factors.  The  grid  spacing  throu|^  the  thickness  of  the 
shell  consists  of  the  five  mesh  points  whose  positions  are  dictated  by  the  Gauss 
quadrature  (see  Appendix  B).  Note  that  since  the  normalised  Young's  modulus 
and  Poisson's  ratio  appear  under  the  integral  sign  of  equation  (b-2),  the 
program  will  permit  the  analysis  of  non-homogeneous  but  isotropic  materials. 
Specifically,  the  program  has  been  written  for  v  and  f/c^  varying  thro\igh  the 
thickness  (y^)  but  uniform  in  the  y  direction. 

Integration  in  the  y^  direction  is  then  achieved  by  svunming  over  all 
stations  along  the  generator.  The  mesh  spacing  along  the  y^  can  be  chosen  by 
the  operator  to  be  either  a  uniform  mesh  of  n  -f  1  points  along  the  entire  length 
of  the  shell  or  a  non-uniform  mesh  representing  a  membrane  region  for  the 
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first  two-thirds  of  the  shellt  and  a  boundary  layer  region  for  the  remaining 
one -third  of  the  shell  near  the  base.  For  the  latter  case  the  spacing  of  the 
boundary  layer  region  is  taken  to  be  one-fourth  that  of  the  membrane  region 
and  the  total  nxunber  of  spaces  n  must  be  divisible  by  three.  Special 
differencing  formula  must  be  used  when  at  the  jxmcture  point  and  note,  in 
particular,  that  equation  (b-1)  will  be  replaced  by 


=  a 


J  juncture 


,  Ti'Tn , 

•  • 


+  b 


forward 


•  • 

<  n'Tti. 

•  .. 


backward 
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where  a  and  b  are  appropriate  weighting  factors. 

In  order  that  the  forward  and  backward  difference  schemes  may  be 

applied  to  all  stations  along  the  generator  including  the  first  and  last  point, 

it  has  been  necessary  to  create  two  fictitious  points,  number  zero  and  n  +  2. 

These  points  are  subsequently  eliminated  by  defiiting  their  contributions  as 


Wg  =  2wj 

-  W2  . 

^n+2 
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-  ^n 

Vq  =  2vi 

-  ^2  ’ 
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=  2^n+l 

■  ^n 

Uq  =  2uj 

-'^2  * 

^"0+2 
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when  no  other  boundary  or  closure  conditions  take  precedence. 

As  has  been  already  discussed,  the  program  for  the  radome  ogive  is 
in  some  respects  more  general  than  the  program  for  the  cone  in  that  the  shell 
material  need  not  be  restricted  to  a  homogeneous  one,  and  the  operator  may 
select  either  a  completely  uniform  mesh  or  one  representing  the  membrane  and 
boundary  layer  regions  of  the  shell. 

In  addition  to  these  options,  the  operator  may  also  select,  by  proper 
choice  of  sense  switches, Love's  first  approximation  to  the  ogive  analysis. 
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This  approximation  eliminates  all  a/R^  and  z/R  terms  in  the  formulation  of  the 
problem  (see  Reference  1).  Another  option  available  to  the  operator  permits 
the  solution  of  a  non-truncated  cone  problem. 
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3.  BOUNDARY  CONDITIONS 


The  energy  formulation  of  the  previous  section  results  in  a  set  of 
simultaneous  equations,  and  for  n  -f  1  grid  points  along  the  generator,  there 
are  3(n  +  1)  equations  for  the  w,  v  and  u  displacements.  A  special  case  is 
the  zeroth  harmonic  whe  re  by  symmetry  all  the  u  displacements  are  zero. 
The  number  of  equations  is  then  reduced  to  2(n  +1).  In  the  following 
discussion,  the  general  case  will  always  be  assumed. 

The  matrix  A  of  the  simultaneous  equations  is  singular  and  the  rank 
is  3  n  4-  1.  By  physical  reasoning,  it  can  be  seen  that  such  must  be  the  case 
because  the  cone  is  free  in  space  and  an  arbitrary  rigid-body  motion  can  be 
added  to  the  displacement  vector.  Hence,  the  displacement  vector  is  not 
unique.  This  condition  is  easily  changed  by  specifying  the  base  as  a  plane  of 
reference  for  the  v  and  u  displacements.  Hence,  the  Vq^|  and  rows  and 
columns  of  the  A  matrix  are  eliminated.  The  displacement  vector  of 
equation  (5),  which  is  specifically  written  for  the  cone,  then  becomes; 
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For  both  the  cone  and  ogive,  the  boundary  constraints  are  imposed  using 
an  energy  technique:  that  is,  the  various  boundary  conditions  are  represented  by 
springs  loading  the  edges  and  the  strain  energy  of  the  springs  is  added  into  the 
energy  matrix  for  the  shell.  The  various  springs  are  shown  in  Figure  2 
and  the  method  of  applying  these  springs  will  be  demonstrated  for  the  base 
spring  kj.  Specifically*  the  cone  geometry  will  be  considered  in  this 
discussion  unless  indicated  otherwise. 


The  numerical  value  of  each  spring  constant  determines  the  method  of 
modifying  the  A-matrix.  If  the  spring  constant  is  zero,  the  digital  program 
makes  no  modification  to  the  A-matrix.  For  finite  values  of  the  spring,  the 
energy  contribution  is  computed  and  added  to  the  A-matrix.  For  exaunple. 
the  energy  contribution  for  k^  is 


Jv^-i 


cos^  (mO) 


=  2irC, 


Or  in  finite  difference  form  equation  (10)  becomes: 


Fj‘)  -  iiiiiEsslU  (w„.i  -  wj 
6^ 


e2 


^Wn,  Vu 


•“n-^n+1 


10  0-1 

0  0  0  0 

0  0  0  0 

-10  0  1 


"n 

^n+1 
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ix^ere  K|^=  k|^/ci  ^  and  6  is  the  spacing  of  stations  in  the  membrane  region. 
Hence,  the  matrix  in  equation  (11)  is  added  directly  to  the  *  'iWn+l 

elements  of  the  A-matrix.  If  y  is  defined  as  the  4  x  4  matrix  in  equation  (11), 
the  spring  energy  can  be  written  as  a  partitioned  matrix  in  terms  of  the 
displacement  vector  X: 


Of  course,  before  making  an  energy  addition,  any  constants  factored  out  of 

matrix  A  must  likewise  be  factored  from  the  matrix  .  In  the  cone  program, 

vC  c  i  ^ 

the  factored  constant  is  - . 

l-v2 

The  above  procedure  applies  to  the  other  base  spring  constant  and 

all  the  spring  constants  which  are  finite  at  the  location  of  the  cone  truncation. 

In  Table  I  of  Appendix  C,  the  y  matrices  are  given  for  each  spring  constant 
K^.  Also  given  in  this  table  are  the  definitions  of  the  dimensionless  and 
the  assumed  form  of  the  displacement  vector. 

In  theory,  the  procedure  of  adding  the  energies  F^^ )  to  the  A-matrix 
still  applies,  when  the  constraint  becomes  infinitely  rigid;  that  is,  K — ¥  oo. 
However,  numerical  difficulties  occur  when  very  large  values  of  are  used. 
In  order  to  avoid  the  difficulties  special  consideration  is  given  to  the 
infinitely  rigid  constraints  by  using  connection  procedures.  For  example, 
when  Kj  is  to  be  considered  infinitely  rigid  equation  (10)  implies: 

^1=0  .  (13) 
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Or  in  finite  difference  form  equation  (13)  becomes 

^n+1  =  ^n 

To  impose  the  condition  of  equation  (14),  a  connection  matrix  C  which  is  a 
rectangular  matrix  of  size  n+1  by  n  is  defined  as  follows: 

X  =  cx 


where 


CX  = 


I  0 


0  Yc 


(15a) 


1  0  0 
0  1  0 
0  0  1 
1  0  0 


(15b) 


The  coiuiection  matrix  C  is  used  to  eliminate  from  equation  (3)  the  dis¬ 
placement  Wj^^j  in  terms  of  Wjj  in  the  following  manner: 

x'  AX  =  (W-T)X  (3) 

Now  by  equation  (15) 

i  x'  (C*  AC)X  =  (W-T)CX  (16 


and  adter  taking  the  variation  of  equation  (16) 
(C'AC)X  =  WC-TC  . 
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Equation  (L7)  expresses  a  set  of  sunuLtaneous  equations  for  w,  v  and  u  with  the 
required  boundary  constraint  of  equation  (14)  imposed.  The  set  of  simultaneous 
equations  can  be  solved  for  X,  and  the  desired  displacement  vector  X  obtained 
from  equation  (15). 

In  Table  II  of  Appendix  C,  the  submatrices  ^  connection  matrices 
are  given  for  the  various  combinations  of  the  infinite  constraints  K^.  A  s  shown 
in  the  table,  only  certain  spring  constants  are  associated  with  a  specified  value 
of  the  harmonic  number  m.  If  m  =1,  only  and/or  are  pertinent  to  the 
connection  conditions  at  the  trvincation,  and  the  program  ignores  the  springs 
Kj  ;  j  =  5,  7.  For  certain  combinations  of  rigid  constraints,  the  implied  con¬ 
nection  conditions  are  trivial  such  as  U|^  =U2  =0  for  =00  .  Under  such 

circvim stances,  a  connection  matrix  is  not  defined  and  the  conditions  are 
imposed  upon  the  matrix  and  forcing  vector  by  eliminating  the  required  rows 
and  columns. 

As  yet,  no  mention  has  been  made  of  the  closed  cone  boundary  conditions, 
auid  how  they  are  imposed.  The  base  springs  K|^  and  K2are  handled  in  the 
usual  manner.  For  the  nose  of  the  cone,  special  boundary  conditions  were 
derived  in  Part  A.  These  conditions  are  expressed  as  connection  matrices  for 
the  digital  program  and  used  in  the  manner  already  described.  The  connection 
matrices  are  different  for  different  harmonic  numbers  and  are  summarized  in 
Table  III.  As  a  passing  parenthetical  remark,  notice  that  for  m  =0,  the 
connection  matrix  for  the  closed  cone  in  Table  III  is  identical  with  the  con¬ 
nection  matrix  of  Table  II  for  sKg  >  oo.  This  implies  that  the  cone  with  a 
clamped  boundary  and  a  very  small  truncation  behaves  like  the  closed  cone. 
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This  is  a  coincidence,  since  the  method  of  deriving  the  closed  cone  condition 
was  independent  of  the  spring  conditions.  As  shown  in  the  tables  for  m  >  0, 
the  closure  condition  does  not  correspond  to  any  of  the  spring  constraints. 

In  the  cone  program,  all  of  the  springs  shown  in  Figure  2  ,  are 
applicable  to  the  analysis.  However,  the  cone  program  requires  al¬ 
ways  to  be  infinite.  For  the  ogive  program,  only  springs  and  K2  at  the  base 
are  applicable,  since  the  shell  is  assumed  to  be  closed  at  the  nose.  In 
Appendix  C,  Tables  IV  and  V  summarize  the  energy  matrices  Y  and  connection 
matrices  Yc  for  the  spring  constants  K|  and  K^,  and  Table  VI  shows  the 
closure  matrices.  Particular  note  should  be  taken  of  the  closure  matrices 
which  must  not  only  impose  the  closure  conditions  on  the  elements  of  the  first 
few  mesh  points  but  also  remove  the  fictitious  elements  w^,  Vq,  and 
As  mentioned  in  the  previous  section,  the  fictitious  mesh  points  at  the  nose 
and  base  were  used  as  a  convenience  in  differencing  the  energy  expression. 

To  remove  these  points,  the  following  relation  is  used: 

fo  =  2fi  -  f2  (18) 

where  f  can  be  w,  v,  or  u.  Equation  (18)  has  been  included  in  the  closure 
matrices.  Likewise  at  the  base^the  fictitious  point  n-fZ  is  removed  by 


^n+2  ■  ^^n+1  ’^n  * 

Equation  (19)  has  been  included  in  the  constraint  matrices  for  Kj  and  K2. 


(19) 
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4.  SOLUTION  OF  THE  LINEAR  SIMULTANEOUS  EQUATIONS 

For  the  cone,  the  maximum  number  of  eimultaneous  eqiiations  is  216, 
while  for  the  ogive  the  nvunber  of  equations  increases  to  453.  With  systems 
of  this  size,  the  method  of  solution  mast  be  selected  with  care,  if  the  com¬ 
putational  time  and  costa  are  to  be  within  reason.  Standard  methods  of 
elimination  for  the  solution  of  the  simultaneous  equations  require  approxi- 
mately  i  /3  operations  where  i  is  the  number  of  equations  and  an  operation 
is  considered  to  be  a  multiplication  or  division.  Indeed,  the  computational 
time  would  be  prohibitive  despite  the  high  speed  computers  of  today.  For  i=450, 
approximately  5  hours  would  be  required  to  carry  out  just  the  required 
multiplications  and  divisions  on  the  IBM  704.  Obviously,  solutions  must  be 
made  using  mt^thods  which  take  advantage  of  the  convenient  form  of  the 
equations. 

In  Section  2,  it  was  indicated  that  by  properly  arranging  the  rows  and 
columns  of  the  A-matrix  all  non- zero  elements  can  be  located  in  15  diagonal 
rows.  Such  a  matrix  is  in  a  convenient  form  for  solution  by  the  block  tri¬ 
diagonal  method  vhich  permits  skipping  the  operations  on  the  zero  elements. 

13  2 

The  operation  count  for  the  solution  is  then  approximately  ij  where  j 
is  the  size  of  square  submatrices  to  be  defined  subsequently.  By  taking  the 
ratio  of  the  operation  coimt  of  the  block  tri-diagonal  method  and  the 
standard  elimination  method,  we  find: 
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A*  mentioned  earlier,  the  A>matrix  is  symmetric  positive  definite  and  as  will 
be  shown,  tiiis  is  a  sufficient  condition  to  assure  that  is  non- singular.  The 
proof  follows  from  the  fact  that  all  diagonal  minors  of  a  positive  definite  matrix 
are  non- singular, or  if 


.  \ 
%  ' 


«  % 


\  C,-! 


A'  '  ^  0  for  f  a  1,  2,*  •  •  k 


Now,  proceeding  by  induction  assuming  aj(j=l,  •  •  *1-1)  exists,  we  compute  oj 
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Therefore 


I 


and  is  non- singular  for  all  i. 

When  solving  a  large  system  of  equations,  the  growth  of  round-off  error 
and  inherent  errors  is  a  familiar  difficulty  and  has  provided  the  field  on 
numerical  analysis  with  a  bottomless  problem  which  has  absorbed  tremendous 
amounts  of  hxunan  effort.  This  error  difficulty  persists  and  has  not  been 
resolved  by  large  computers;  rather,  it  has  been  reduced  to  a  lower  order  of 
magnitude  by  the  use  of  double  or  triple  precision  operations.  However,  the 

i 

numerical  analyst  must  always  concern  himself  with  the  fact  that  some  matrix 
will  be  encountered  that  will  cause  unpermissible  error  in  the  result.  Being 
aware  of  these  facts,  a  method  of  improving  the  solution  was  devised  for  the 
digital  program.  The  method  consists  of  computing  the  first  solution  and 

i 

using  this  to  compute  a  residual  R^,  which  in  turn  is  used  to  obtain  a 

I 

correction  to  the  original  solution.  The  procedure  is  repeated  vintil  a  specified 

i 

I 

number  of  significant  figures  is  achieved  in  each  component  of  the  i-dimensional 
solution.  Symbolically  the  method  is  as  follows: 

x„  =  A-‘b 

Ro  =  b  -  AXo 

I 

AX^  .  A‘'Ro 
Xl  =  X„+AX„ 
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5. 


EVALUATION  OF  STRESSES 


At  each  mesh  point  along  the  generator,  the  stresses  of  eqiiation  (27) 
are  evaluated  at  the  outside  surface,  the  znid-plane,  and  the  inside  surface 
(z  s  'h/Z,  0,  h/2).  In  order  to  evaltiate  the  stress  at  a  particular  point, 
the  product  of  the  temperature  difference  and  a  must  be  provided  as  an  input. 
An  exception  for  'which  the  program  does  not  require  an'y  temperature  inputs 
is  the  isothermal  case  with  pressure  loading.  The  stresses  are  also  functions 
of  the  displacements,  which  are  obtained  from  the  solution  of  the  simultaneous 
equations,  and  the  derivatives  of  the  displacements  which  are  computed  b'y 
finite  difference  techniques.  It  should  be  recalled  that  the  cone  has  two 
boundar'y  regions  and  membrane  region.  The  computation  of  the  derivatives 
at  junction  points  between  these  regions  requires  special  consideration. 

For  the  first  derivativea  a  standard  3-point  centered  difference  form’xla 
is  used  for  all  internal  points  removed  from  the  boundaries  and  the  junction 
points.  At  the  edges  of  the  cone,  a  3-point  backward  or  forward  differencing 
formula  is  used  to  compute  the  first  derivative.  The  forward  difference 
would  be  used  at  the  left  edge,  and  the  backward  difference  at  the  right  edge 
or  base.  To  compute  the  first  derivatives  at  the  junction  points  of  the  two 
regions,  a  3-point  for'ward  or  backward  difference  formula  is  again  used  and 
the  computation  is  al'ways  based  upon  data  in  the  fine  mesh  region. 

The  second  derivatives  are  computed  in  a  nnanner  similar  to  the  first 
derivatives.  A  3-point  centered  difference  formula  for  the  second  derivative 
is  used  at  all  points  removed  from  the  edges  or  junction  points.  A  4-point 
forward  or  backward  differencing  formula  is  used  at  the  edges  and  junctions 
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with  the  c<»nputation  based  on  data  from  the  fine  mesh  region. 
The  physical  stresses  for  the  ogive  are  given  by: 

CT  ^  .nnD  I  ^  x  >  _ 2  ^ 


*  o  *  < 


cos  mO 

r 

1 

1  - 

Rj  cosp 

“”K‘- |W  - 


+  v  -z  -  a  (l+v)  (T-To) 

R  '  R  '  8y  ay* 


O'  e  cosmOj  v  \  /  ,  nr  z  \  •  ^ 

- - — ^ -  U" — 1  mu  +  Bind  (I - iv-icosd - Vw-Bin<3z  — 

«o  «o  I-  V*  I  R,co8P|j-  r/  ^  R'  \  Rjcos.d/  9y 

4-  ^  M  I  -4)  .  1^1  -  .  (  1  +  V)  (T-T„)l 


g-  €  sin  m9  )  cos^q  sinp 

«o  *o  2(l+v)  |co8*PR/  rJ  cos^P 


+  I  4.  dw 

R]  8y  Rjcosp  8y 


+ - ^ (l  sin^  u  -  m  (l-  v  +  ^  8inft  ^  8^  | 

Ricosp'^  R/  R'  R,co8p  8yj 

(28)'* 

The  finite  difference  formulas  for  the  derivatives  in  the  above  expressiom 


are  the  same  as  for  the  cone  program.  The  stresses  themselves  are  calculatec 
at  the  mid-plane  and  outer  and  inner  surfaces  and  require  values  of  Poisson's 
ratio  and  a  normalized  Young's  modulus  at  these  points. 
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6. 


NUMERICAL  RESULTS 


Numerous  problems  have  been  analyzed  using  the  digital  programs,  and 
for  problems  which  have  a  theoretical  solution  the  analytical  and  numerical 
results  have  been  compared.  Theoretical  results  are  available  for  comparison 
with  several  interesting  cases  for  tlie  zeroth  harmonic  with  pressure  or  thermal 
loading.  Five  of  these  cases  are  discussed  in  this  section.  For  harmonic 
numbers  greater  than  zero,  the  number  of  theoretical  solutions  are  limited 
and  in  this  report  only  one  problem  with  m  ~  1  is  discussed. 

For  the  cone,  the  following  problems  are  discussed: 

Problem 

No.  m  X  a  Loading 

1  0  .  1489  -  Closed  cone  -  constant  temperature  increase 

2  0  .  1489  . 1  Truncated  cone  -  constant  pressure 

3  0  .  1489  . 1  Truncated  cone  -  constant  temperature  increase 

4  0  .4636  .  1  Truncated  cone  -  Linear  temperature  along  ^ 

5  1  .6435  .1  Truncated  cone  -  Linear  temperature  along 

diameter 

In  all  of  the  above  problems,  the  (non-dimensional)  wall  thickness  and  Poisson's 
ratio  are  assumed  constant  and  equal  to  0.  00206  and  0.  25,  respectively. 
Furthermore,  only  the  results  with  72  spaces  along  the  generator  are  con¬ 
sidered. 

Problem  1  is  the  analysis  of  a  closed  cone  uniformly  heated  with  a 
temperature  change  of  lOOO^F.  The  base  is  assumed  to  have  no  edge  moment 
but  edge  shear  prevents  any  radial  growth.  This  is  equivalent  to  =  0  and 

'ik-JD  7B 
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K2  =  00.  Aa  ahown  in  Figure  4 ,  the  deflectiona  w  end  v  are  amoothly  varying 
curvea  over  the  73  meah  pointa.  The  abaciaaa  haa  been  labled  with  two  acalea. 
The  upper  acale  indicatea  the  y  diatance  along  the  generator  while  the  lower 
acale  indicatea  the  meah  point  number.  The  firat  meah  point  ia  located  at 
y  s  .  0069  which  ia  alightly  larger  than  the  diatance  to  the  inaide  apex.  The 
firat  bovindary  region  enda  at  atation  25  and  the  aecond  boundary  region  atarta 
at  atation  49..  The  membrane  region  between  atationa  25  and  49  haa  a  meah 

apacing  of  6  ,  u^iere  6  =  .0276.  Aa  mentioned  in  Section  2,  the  pointa  in  the 

i 

boundary  regiona  have  a  uniform  apacing  of  6/4. 


If  the  baae  of  the  cone  ia  completely  unreatrained,  the  cloaed  cone  with 
uniform  heating  ia  theoretically  in  a  zero-atreaa  atate.  Further,  the  dia- 
placementa  are  found  to  be: 


w  =  -  aT  tan  \  =  -  .  15001024  x  10 

-3 


-2 


(29) 


V  =  -  aT  (l-y)  =  -  10  X  10  (l-y) 

where  a  =  10~^,  the  value  uaed  in  the  digital  program.  The  nmnerical  reaulta 
from  the  program  agree  well  with  the  theoretical  valuea  from  equation  (29), 
for  all  valuea  of  y  away  from  the  baae  region.  The  miudmum  error  waa  foimd 
to  be  3  parta  in  1500  or  0.2%.  In  the  vicinity  of  the  baae  mounting,  the 
theoretical  aolution  ia  not  applicable  becauae  the  baae  ia  aaaumed  to  be  free 
whereaa  in  problem  I  the  baae  haa  aome  conatraint.  The  numerical  aolution 
with  the  damped  oacillationa  in  w  ia  typical  of  the  edge  phenomexton  of  ahella. 

Alao  ahown  in  Figure  4  ia  the  atreaa  along  the  generator.  The 

atreaa  ia  in  dimenaionleaa  form  and  muat  be  multiplied  by  Young'a 
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modulu*  <  to  obtain  the  phyeical  streas.  The  atreea,  like  the  w  displacement, 
has  a  damped  oscillatory  character  in  the  vicinity  of  the  base,  and  away  from 
the  base  the  stress  varies  from  30  to  300  psi  for  c  =30  x  10^.  In  theory, 
the  stress  should  become  zero  away  from  the  base.  In  the  immediate 
proximity  of  the  nose,  the  stress  has  not  been  plotted  because  the  structure 
violates  the  shell  assiunptions  and  behaves  like  a  3 -dimensional  body. 

The  mid-plane  stress  and  the  stresses  off  the  mid-plane  are 
similar  in  nature  to  0^  ;  that  is,  the  stresses  are  large  in  the  vicinity  of 
the  base  and  become  small  in  the  membrane  region.  These  stresses  are 
found  to  be  in  good  agreement  with  theory. 

For  problem  2,  the  cone  has  a  10%  truncation  and  a  uniform  pressure 
loading.  The  base  boundary  conditions  are  =  0  and  oo, while  at  the 
location  of  the  truncation  no  constraint  is  imposed;  that  is,  0. 

The  deflections  w  and  v  for  problem  2  are  plotted  iu  Figure  5  as  the  solid 
curves.  The  circled  points  represent  the  theoretical  solution  to  the  same 
problem  based  upon  the  following  equations  from  membrane  theory: 

p  tan*  X  r  3  2  .  /.i  .  I  .  2  .  1 

w  =  —  +  V  (a  -1)  +  --  +  a  In  yj 

r  ,  -1  (30) 

V  =  -^  H  V  -  -  )  (y*  -  1)  +  a*  In  yJ  . 

For  the  same  problem,  the  mid-plane  stresses  ond  are 

plotted  in  Figure  6,and  the  stresses  (7*^  and  (as  -  h/2  )  are  plotted 

in  Figure  7.  Again  the  circled  points  are  the  results  of  analytical 
calculations  for  selected  points  along  the  generator.  The  agreement  is  seen 
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to  be  good.  For  the  streeeee  0*^  ^3  ^ich  have  not  been  plotted,  the 

numerical  and  analytical  results  were  found  to  be  in  good  agreement. 

In  Appendix  D,  the  complete  set  of  inputs  and  outputs  is  shown  for 
problem  2.  Also  included  in  Appendix  D  is  the  program  outputs  for  problem  3. 
Both  these  problems  have  the  same  cone  geometry  but  different  loads  and 
boundary  conditions. 

In  problem  3,  the  truncated  cone  is  analyzed  for  a  uniform  increase  of 
temperature.  At  the  base  =  0  and  K2  =  oo,  while  at  the  trvmcation  the  edge 
is  ass\uned  to  have  no  radial  growth  or  change  of  slope.  This  corresponds  to 
Kg  s  s  00.  In  Figure  8  .  the  displacements  w  and  v  are  shown, and  the  edge 
phenomenon  is  again  noted.  In  the  membrane  region,  the  theoretical  dis¬ 
placements  are  computed  using  equation  (29)  and  the  error  in  w  found  to  be 
roughly  3  parts  in  1500.  Theoretically,  the  stresses  should  be  zero  in  the 
membrane  region  and,  as  shown  by  the  tabulation  in  Appendix  D,  the  computed 
stresses  are  small,  varying  from  30  to  300  psi  for  c  =  30  x  10^.  The  agree¬ 
ment  is  considered  satisfactory. 

In  problem  4,  a  truncated  cone  is  thermally  loaded  with  a  temperature 
distribution  which  is  linear  along  the  axis  of  the  cone.  If  the  edges  of  the  cone 
are  free  from  constraints,  the  theoretical  displacements  are  found  to  be: 


.  _  m  tan  V  ,,  2  , 

w  =  -  oT  — - —  (I  -  y  ) 

v  =  -  aT  (y  -  l)V2  • 


(31) 
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Furthermore,  the  cone  is  found  to  be  in  a  zero  state  of  stress. 

The  displacements  vr  and  v  from  the  digital  program  are  shown  in 
Figure  9  .  The  circled  points  are  obtained  from  equation  (31),  and  are  seen 
to  agree  well  with  the  numerical  results  from  the  program.  The  stresses  at 
the  mid-plane  and  off  the  mid-plane  were  found  to  be  small  and  in  satisfactory- 
agreement  with  the  theoretical  zero  state  of  stress. 

The  last  example  to  be  considered  for  the  cone  is  problem  S  where  the 
first  harmonic  of  the  loading  is  used.  The  loading  is  a  linear  temperature 
distribution  along  a  diameter  of  the  cone.  The  cone  is  free  at  the  truncation 
and  constrained  from  radial  growth  at  the  base  (K2  =  00).  Because  of  these 
boundary  conditions,  it  is  not  possible  to  obtain  an  amal-ytical  solution  which 
can  be  compared  with  the  w,  v  and  u  displacements  shown  in  Figure  10  . 
However,  the  stresses  in  the  membrane  region  shovild  be  zero  or  at  least 
small.  The  stresses  from  the  digital  program  varied  from  30  to  300  psi 
(  c  =  30  X  10^)  in  the  membrane  region  and  are  considered  to  be  in  satisfactory 
agreement  -with  theory. 

One  of  the  first  problems  analyzed  with  the  ogive  program  was 
that  of  a  uniform  pressure  loading  on  a  shell  with  free  edge  constraints. 

A  uniform  mesh  ha-ving  73  points  was  selected.  The  folio-wing  table  summarizes 
the  inputs  to  the  problem;  note  that  all  lengths  have  been  scaled  by  the  length 
of  the  generator: 
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.25  radian 

»<»)  = 

.  25  for  all  z 

4.0 

«  (a)  /«o  * 

1 . 0  for  all  s 

.00125 

P/«o  = 

.25  X  10"® 

0 

Ki  =  K2  = 

0 

m  =  0 

The  reeiilting  deflections  are  shown  in  Figure  LI.  As  a  comparison  the 
membrane  solution 


=  <P/«o> 


4  h 


cot  ^  ain  |cot^  cac^  in  (cot4/2)J 


+  T  •""  ♦o 


cot  4  csc^^  -f  2  cot  ^1  >•  sin  4  ain  4| 


► 


W  =  V  cot  4  ” 


(pAq)  R* 


I  T  7 


L' 


i  ■ 


sin  4  -  sin  4cg 

sin  4 


♦  =  T 


has  been  calculated  for  several  values  of  4  and  as  seen  in  the  graph,  the 
agreement  is  good.  It  is  to  be  particularly  noted  that  the  maximum  v  and 
minimum  w  do  not  occur  at  the  extremity  of  the  nose,  but  are  located  at  about 
station  number  4.  This  phenomenon  is  predicted  by  theory.  The  ogive 
program  has  not  been  used  to  compute  stresses  since  that  part  of  the  program 
is  undergoing  modification.  Stress  results  will  be  described  in  subsequent 
scientific  reports. 
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7.  CONCLUSIONS 

The  digital  programs  enable  the  design  engineer  to  make  detailed 
investigations  of  the  structural  requirements  for  ogival  and  conical  radomes 
under  thermal  and/or  pressure  loading.  The  programs  can  be  used  to  analyze 
the  radome  for  any  flight  mission,  provided  the  required  thermal  and  pressure 
loading  is  known.  For  a  defined  load,  the  program  provides  the  design 
engineer  with  (1)  the  displacements  w,  v,  and  u  at  the  mid- surface  of  the 
shell,  and  (2)  the  normal  and  shear  stress  at  z  =  -h/2,  0,  and  h/2.  This 
information  is  given  for  n  1  points  along  the  generator. 

Besides  investigating  various  loading  conditions,  the  programs  can  be 
used  to  investigate  different  methods  of  mounting  the  radome.  In  the  cone 
program,  no  radial  growth  at  the  base  is  allowed,  but  any  degree  of  edge 
moment  constraint  can  be  specified  (K^  &  0,  finite  or  oo  and  K2  =  00).  At 
the  triincation  of  the  cone,  any  arbitrary  constraint  can  be  imposed  by 
specifying  the  proper  combination  of  spring  constants.  In  all,  there  are  5 
springs  at  the  truncation)  each  of  which  can  be  specified  as  zero,  finite  or 
infinite.  Hence,  a  conical  radome  with  a  nose  boom  can  be  investigated  for 
various  amoimts  of  flexibility  in  the  boom.  In  addition,  the  cone  program 
can  be  used  for  the  analysis  of  the  closed  nose  structure.  The  ogive 
program  has  been  constructed  for  the  closed  nose  condition  only.  At  the 
base  of  the  ogive,  completely  general  constraints  are  permitted;  that  is, 
both  Kj  and  K2  can  be  zero,  finite  or  infinite. 
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All  of  the  numerical  results  discussed  in  Section  6  are  for  the  spring 
constants  zero  or  infinite.  Problems  with  finite  springs  at  the  edges  are 
omitted  since  suitable  analytical  solutions  are  not  available  for  comparison 
purposes.  However,  the  programs  have  been  used  to  analyze  systems  with 
finite  springs,  and  as  the  value  of  the  spring  constants  was  varied  from  zero 
to  larger  and  larger  values,  the  results  were  observed  to  change  in  a 
consistent  manner  from  the  zero  constraint  case  to  the  infinitely  rigid  case. 

The  numerical  results  indicate  that  the  digital  program  for  the  cone 
accurately  computes  the  displacements  and  stresses  for  pressure  and/or 
thermal  loading.  Any  discrepancies  between  the  numerical  and  theoretical 
results  are  indeed  small  and  well  within  the  requirements  of  the  design 
engineer.  Similar  remarks  are  applicable  to  the  ogive  program  regarding 
the  displacements. 

It  is  concluded  that  the  combination  of  the  above  cone  and  ogive  programs 
make  possible  a  detailed  examination  of  any  conical  or  ogival  radome  for 
all  nose  and  base  boundary  conditions,  material  properties  and  load  dis¬ 
tributions  of  practical  interest. 
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APPENDIX  A 


FINITE  DIFFERENCE  FORMULAS  FOR  THE  CONE  PROGRAM 


The  free  energy  expression  as  shown  by  equation  (7)  involves  zeroth, 
first,  and  second  derivatwes  of  the  displacements  and  products  thereof.  For 
all  of  the  terms  not  involving  a  second  derivative,  the  following  difference 
approximations  are  used: 


r  N 

J  (r')^dy  =^(ri+i-ri)2  .  -j- 

r  .N_  ^ 

I  r's'dy  ^  (n-l-r^iX^i+r^i) 

J  I  i 

Jrs'dy  =  7^(ri+i+ri)(Si^l-Si) 


where  6^  is  the  spacing  of  the  stations  at  the  ith  mesh  point,  r  and  s  are  any 
of  the  displacement  w,  v,  or  u,  and  N  is  the  number  of  spaces  along  the 
generator.  The  above  fornulas  apply  without  any  modification  to  a  uniform  or 
variable  mesh  size. 


For  terms  involving  the  second  derivative  of  a  displacement,  special 
consideration  must  be  given  to  the  junction  points  between  the  membrane  and 
boundary  regions.  Let  the  junction  points  be  called  N1  and  N2  for  the  left  and 
right  boundary  regions,  respectively.  The  finite  difference  forms  are  then 
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written  ai: 


.2  .  <*i+l- 2ei  +  ei.l)*  .  8  I  _  ,  ,2 

q  ■•■  Y  'ni  ®*Ni  ■’■  "ni+i* 


r(ei»)  dy  rj 


^  T  6»  'N2  ^‘nZ-I  "  ^■n2'*‘ ^'nZ+I  ^ 


J  2  6  *  26"  '  Nl+l  Nl  Nl-f'^l-l  Nl  Wl' 


t_ 

26* 


+  -  ^^'^NZ+l  ■  ^'N2  ■*■  'nI-1  ^  ^“NZ-I  “  ®*N2  ■*■  ^“nZ+I  ^ 


where  a  and  b  have  the  following  valuee  for  the  various  regions 


left  boundary  region 
membrane  region 
right  boundary  region 


a  s2  ;  b  -Nl-l 
a  3 Nl+l;  bsN2-l 
a  =N2+lf  b  =  N  . 
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APPENDIX  B 


FIVE-POINT  GAUSS  QUADRATURE  COORDINATES  AND  WEIGHTS 


s 

-.906179846 

=  .236926885 

qz 

= 

-.538469310 

a2  =  .478628670 

<13 

s 

0 

a3  =  .568888889 

= 

.538469310 

a^  =  .478628670 

*15 

= 

.906179846 

ag  =  .236926885 
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APPENDIX  C 

table  I 

Enersy  Contribution  from  the  Boundaries  of  the  Cone 


k2 

K2=7- 


^3 


K4=73 

<1 


^4  16(1-1^ 


ks 

K5=-^ 

cl 


(l-v‘')a  sin  X  j 


cos**X  -sinX  cosX  W|^ 

-  sinX  cosX  sin^  i 

‘ .  .  _  ^  I 

sin^X  sinX  cosx]  fw 

sinXcosX  co»\  V 

UiKiYi 
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VA  I  A  V 


TABLE  n 

Connection  Conditions  for  Variou«  Combliuittlon»  of 
Infinite  Sprins  Constanta  for  the  Cone 
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TABLE  n  (Continued) 
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TABLE  m 


Nose  Closure  Matrices  for  Cone 


Vg  =  0) 


Yc  (nx  .=  1) 


■ 

r  ■' 

r 

r 

0 

1 

“1 

(1-q) 

0 

q  cosx 

0 

COtX 

^2 

tanX 

0 

-  cscX 

1 

0 

i 

V2 

► 

q  cosx 

0 

(1+q  cos^X) 

* 

0 

1 

"2 

1 

0 

0 

• 

0 

1 

0 

s 

s 

0 

0 

I 

. 

V 

""2 

'^2 

U2 

W3 


O 

where  q  =  —————  a.nd  yj  =y  at  station  2. 
4  yj  sin  X 


TABLE  IV 

Energy  Matricea  for  the  Bm«e  Connection  Springs 
Moment  Spring,  modulus 


Radial  Growth  Spring,  modulus  K2 


K, 


0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 


0 

0 

0 

0 

0 

0 


w 

0 

0 

1  1 
~Z^  2R6 

0 

0 

0 

0 

r 

'^n 

0 

0 

0 

0 

0 

0 

0 

0 

0 

^n 

0 

0 

0 

0 

0 

0 

0 

0 

0 

i 

!“n 

1 

26^ 

0 

0 

^  0 
■p- 

0 

1 

2?" 

0 

0 

^n+1 

1 

2l^ 

0 

0 

0  ^ 

0 

1 

2R6 

0 

0 

< 

V  , 

n+1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

"n+1 

0 

0 

0 

1  1 
26*  2R6 

0 

1 

26* 

0 

0 

^n+2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

^n+2 

J 

0 

0 

0 

0 

0 

0 

0 

0 

0 

,  r 

■%+2| 
•-  J 

0 

0 

0 

0 

0 

0 


r 

^n+1 

'"n+l 

^n+2 

r'n+z 

“n+2 


j 


These  are  for  m  2!  1.  When  m  =  0  delete  rows  and  columns  numbered 

3,  6,  9  in  the  first  matrix  and  rows  and  columns  numbered  3  and  6  in  the 
second  matrix. 
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TABUS  V 


Base  Connection  for  the  Radome  Oaive 


Spring  Combination 


(m  s  0) 


(m  >  1) 


K|  Finite. 
K2  Finite 


^n+1 


10  0  0 

0  10  0 

0  0  10 

0  0  0  1 

0  0  0  0 

0  0  0  0 

10  0  2 

0-100 
0  0-10 


“n 

^n+l 


Infinite. 
K2  Finite 


K|  Finite. 
K2  Infinite 


Infinite, 
K2  Infinite 


^n+1 


{'] 


10  0  0 
0  10  0 
0  0  10 
0  0  0  1 
0  0  0  0 
0  0  0  0 
10  0  0 
0-100 
0  0-10 

1  0  0 
0  1  0 
0  0  1 
0  0  0 
0  0  0 
0  0  0 
-10  0 
0-10 
0  0-1 

~  1  0  o' 
0  1  0 
0  0  1 
0  0  0 
0  0  0 
0  0  0 
1  0  0 
0-10 
0  0-1 


“n 

^n+1 
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TABLE  VI 


Nose  Closure  for  the  Radome  Ogive 


0 

•1 

0 

0 

0 

0 

0 

1 

0 
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APPENDIX  D 

DIGITAL  PROORAM  OUTPUTS  POR  PROBLEMS  2  AND  3 

GOSL'  01.24-3000  PRGB  2  ^.'24-50 

1  7 

INPUT 

CLRSS  8  H#T2  HPPMCfUC*  0.  PINGS  0.I480E-00  Rt  1.0000E-01  THICK.*  0.2060E-02  NU«  0.2500E-00 

SPRING  CONS 

0.  0.  0.  0.  0.  0.  0.90E  04 

PRESSURE 

0. 2500E-0S 
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UTPUTS 


MISSILE  AT  LOCflTIOri  0 


TIME#  0 


'CS 

V 

1 

0.1000 

2 

0. 1062 

3 

0.1125 

4 

0.1187 

5 

0. 1250 

6 

0.  1  312 

7 

0. 1 375 

8 

0.1437 

9 

0. 1500 

10 

0. 1562 

1  1 

0.  1625 

1  2 

0.  1687 

1  3 

0.1750 

1  4 

0. 1812 

15 

0. 1875 

16 

0. 1 937 

1  7 

0. 2000 

18 

0.2062 

1  9 

0.2125 

20 

0.2187 

21 

0.2250 

22 

0.2312 

23 

0.2375 

24 

0.2437 

25 

0.2500 

26 

0.2750 

27 

0. 3000 

28 

0. 3250 

29 

0. 3500 

30 

0. 3750 

31 

0.4000 

32 

0.4250 

33 

0. 4500 

34 

0.4750 

35 

0.5000 

36 

0.5250 

37 

0.5500 

38 

0.5750 

39 

0. 6000 

40 

0.6250 

4l 

0. 6500 

42 

0.6750 

43 

0.  700ij 

44 

0,7250 

45 

0. 7500 

46 

0.7750 

47 

0. 8000 

48 

0. 8250 

49 

0 . 85*00 

50 

0.85L--2 

lil 

0.  33746-»48E“05 
0. 341  U'2!6E-05 
0. 3447?228E-05 
0. 3484.3457E-CI5 
0. 352I?319E-0? 
0. 35820049E-05 
0.  360276(i4E-Ci5 
0. 364484 35E-05 
0.  368826 17E-Ci^ 
0.37330427E-D5 
0. 3779204?'E-05 
0.  33267648E-CI5 
0.  387*.7375E-0'5 
0.  31J2l  1  360E-H5 
0. 387V972 IE-05 
0.40312561E-05 
0.4085?’:'75E-rj5 
0.41422040E-CI5 
0.4t89SS76E“05 
0.425?'07l2E-05 
0.43197882E-05 
0.43820152E-05 
0«  44454ii.03E-05 
0.450?'?253E-0’> 
0.4572?'583E-05 
0.43600TI0E-O5 
0«51667250E-Ci5 
0.54l)79C'76E-05 
0.53538727E-05 
0.62'?47294E-05 
0.664056S4E-CI5 
0.70714654E-C5 
0.75:74825E-05 
0. 800867 18E -05 
0.$5150774E-05 
0.  ■■'0467366E-05 
0.T'60?6823E-‘j5 
0.  10185042E-0H 
0.  I079  354  1E-04 
0,  I  r426501E-0«» 
0.  1  2084i;'4rjE-ijH 
0. t2768576E-04 
0.  1  3477723E-0H 
0.  l42l22-:'4E-04 
0. 14072302E-04 
0. 1575774 3E-04 
0.  1ii5'*8605E-0** 
0. t7405T14E-M4 
0 .  18176  35';'E -  Cl** 
0.  t8**9l4l<‘'E-:** 


V 

0.20692?91E-04 
0.2071E:776E-0*» 
0.20739736E-04 
0. 20756079E-04 
0.20768t09E-04 
0.207760S3E-04 
0.207802t5E-04 
0.207806858-04 
0.20777655E-04 
0. 20771 266E-04 
0. 20761 642E-04 
0.20748S9tE-04 
0.207331 14E-04 
0.2CI7IH395E-04 
0. 2069281 6E-04 
0. 20668449E-04 
0.20641 353E-04 
0.2061 1590E-04 
0.20579214E-04 
0. 205*4  4  27 IE -04 
a.20506Sl0E-04 
0. 20466870E-04 
0. 204244ft4E-04 
0, 20379660E-04 
u. 20332367E-04 
0.201  21  S56E-0‘» 
0.  1987561 IE-04 
0. 19594456E-04 
0. 19279 34  It -04 
0. 18931 028E-04 
0. 185501 31E-04 
0. 181 37150E-04 
0.  1  76'.  2503E-04 
0. 17216536F-04 
O.  16.m:i9545E-04 
0.  1617178;E-04 
0.  1  560  34'*- 3E -04 
0. 15004776E— 04 
0.  t4  375S8  3E-C'4 
0.  1  371<;.92‘'E-C'4 
0. 1 302804  IE -04 
0.  I  2309331E-04 
o.  1  15».'089‘^E-04 
0.  1 07828  35E”0‘* 
0. 99752 186E -05 
0.91  3.r  1  226E  -'.•’I' 
0.827161 12E-U5 
0.73757683E-05 
0.6**508583E-05 
0.  -r  ^l  4  *'599t-05 
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OUTPUTS  t'MSSILE  PT  LOCPTION  0.  TIME*  0. 


PCS 

V 

III 

y 

5t 

0.8625 

0.  1S710462E-':i4 

0. 59771 924E -05 

52 

0,8687 

0.  1892133SE-*0^ 

0.57375689E-05 

53 

0. 8750 

0.  r9154988E-04 

0.54961058E-05 

54 

0. 8812 

0. 19 380871 E-Cw 

0.52528153E-05 

55 

0.8875 

0. t?60?109E-04 

0.50077005E-05 

56 

0. 8937 

0.  19S35470E-04 

0.47607522E-05 

57 

0 . 9000 

0.2006  U16E-U-4 

0.451 19458E-05 

5S 

0. 9062 

0. 2028*4466E-04 

0.42612404E-05 

59 

0.*9125 

0. 2050421 OE-04 

u.4n08‘^!^6tE-n5 

60 

0.^187 

0. 2072 3260E-'J4 

0.  375  39466.E-ri5 

61 

0. 9250 

0.  209501 82E-04 

0. 34973493E-05 

62 

0.9312 

0.21 202806E-04 

0. 32389658E-05 

63 

0. 9375 

0.215101 12E-04 

0. 297921 83E-05 

64 

0,  9*4  37 

0.  21 9091 n7E-04 

0.2718373i>E-05 

65 

0. 9500 

0.22*4  7 1057E-04 

0.2459042tE“05 

66 

0.9562 

0.23070183E-04 

0.22009341E-05 

67 

0.9625 

0.  23729^  ?0E-04 

0.  1 9451 722E-05 

68 

0. 9687 

0.  24 1 45927E-04 

0. 1 6904507E-05 

69 

0. 9750 

0.23S14484E-04 

0.  l43l462$E-05 

70 

0.9312 

i:i.2rr'6i226E-o*4 

0. 1 1 56381 6E -05 

71 

0. 9875 

0.  17657584E-04 

0, 84490 962E-06 

72 

0. 9937 

0.  10212593E-i:i4 

0.46903087E-C*b 

73 

t . 0000 

0. 

0. 
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OUTPUTS 


MISSILE  AT  LOCATION  0, 


TIMES  0 


POS  V 


1 

0. 1000 

2 

0.1062 

3 

0. 1125 

4 

0.1187 

5 

0.1250 

6 

0.1312 

7 

0.1375 

8 

0.1437 

9 

0. 1500 

10 

0.1562 

1  1 

0. 1625 

12 

0,1687 

13 

0. 1750 

14 

0.1812 

15 

0.1875 

16 

0. 1937 

17 

0.2000 

18 

0.2062 

19 

0.2125 

20 

0.2187 

21 

0.2250 

22 

0.2312 

23 

0.2375 

24 

0.2437 

25 

0.2500 

26 

0.2750 

27 

0.3000 

28 

0.3250 

29 

0.3500 

30 

0.3750 

31 

0,4000 

32 

0.4250 

33 

0.4500 

34 

0.4750 

35 

0.5000 

36 

0.5250 

37 

0.5500 

38 

0.5750 

39 

0.6000 

40 

0.6250 

41 

0.6500 

42 

0.6750 

43 

0.7000 

44 

0. 7250 

45 

0. 7500 

46 

0. 7750 

47 

0.8000 

48 

0.8250 

49 

0.8500 

50 

0. 8562 

S1GMA-It,t 

-0. 16791 730E-04 
-0. 18084S01E>04 
-0. 1920021 3E-04 
-0. 2033969 lE-04 
-0.21487043E-04 
-0.22633e00E-04 
-0.23779085E-04 
-0.24923291E-04 
-0.26066673E-04 
-0.27209377E-04 
-0.28351505E-04 
-0. 294931 39E-04 
-0.3063434tE-04 
-0.31 7751 97E-04 
-0.32915725E-04 
-0.34055968E-04 
-0.35I95957E-04 
-0.  36335652E-0-. 
-0. 374751 45E-04 
-0.386I49I9E-04 
-0. 3975621 4E-04 
-0.40e99701E-04 
-0. 42039051 E-04 
-0.43148762E-04 
-0.44214238E-04 
-0.48865071E-04 
-0.53417975E-04 
-0.57970725E-04 
-0.625232 1 7E-04 
-0.67075353E-04 
-0.71627212E-CI4 
-0.761 78863E-04 
-0.807303S3E-04 
-0.85281717E-04 
-0.898329796-04 
-0.943841 53E-04 
-0.98935264E-04 
-0. 1034863IE-03 
-0. 10803733E-03 
-0. 1 1258830E-03 
-0. 1 1713924E-03 
-0. 12I69015E-0J 
-0. l2e24IC4E-03 
-0. I3079I92E-03 
-0.I3534278E-CI3 
-0. I 3989350E-03 
-0. I4444332E-03 
-0. 14899776E-03 
-0. 1536281  IE-03 
-0, 1 54733 19E-03 


SlGHA-11.2 

-0.180300076-04 
-0.19340359E-04 
-0.20463409E-04 
-0.215931 19E-04 
-0.227312656-04 
-0.23870723E-04 
-0. 2500991 8E-04 
-0.26148878E-04 
-0.27287702E-04 
-0.28426433E-04 
-0.29565082E-04 
-0. 30703660E-04 
-0.3I842171E-04 
-0.32930643E-04 
-0.341 19057E-04 
-0.35257421E-04 
-0. 36395729E-04 
-0.37533955E-04 
-0.386722216-04 
-0. 398110056-04 


-0.43224227E-04 
-0.44327519E-04 
-0.45385437E-O4 
-0.50052566E-04 
-0. 546061226-04 
-0.591579796-04 
-0.63709S9SE-04 
-0. 68261 730E-04 
-0. 7281 3489E-04 
-0.773651916-04 
-0. 819168456-04 
-0.864684616-04 
-0.910200446-04 
-0.955715916-04 
-0. 100I23I 16-03 
-0. 104674626-03 
-0.1 092261 tE-03 
-0. 1  1 3777586-03 
-0.  1  IS32903E-03 
-0. I  228804 7E -03 
-0. I2743191E-03 
-0. 1 31 98332E-03 
-0. I 3653474E-03 
-0.1^1 08601 E-0  3 
-0. I 4563684E-03 
-0. I  50194866-05 
-0. I  548' 973E-03 
-0. I 5593071E-03 


SIGMA- 11.3 

-0. 19455606E-04 
-0.20771372E-04 
-0.218897O2E-04 
-0.22999214E-04 
-0.241 19069E-04 
-0.25243206E-04 
-0.26369166E-04 
-0.27496471 E-04 
-0.28624961 E-04 
-0.29754475E-04 
-0. 30S34S73E-04 
-0.32016026E-04 
-0. 37117531E-04 
-0.34280220E-04 
-0. 35413093E-04 
-0.36546403E-04 
-0.37680075E-04 
-0. 38814075E-04 


-0.43355S64E-04 
-0.44479796E-04 
-0.45574704E-04 
-0.466231 36E-04 
-0.51300437E-04 
-0.55849425E-04 
-0.60396013E-04 
-0.64943627E-04 
-0.69491953E-04 
-0. 7404081  •8E-04 
-0.785901 ieE-04 
-0.83139764E-04 
-0.876896036-04 
-0.92230856E-04 
-0.96790206E-04 
-0.101340726-03 
-0. 105801386-03 
-0. 1 104*2156-03 
-0. 1 14903036-03 
-0. 1 19543086-03 
-0.124005036-03 
-0. 128646146-03 
-0. I  331 07306-03 
-o. 1 37748526-03 
-0. 1422006*6-03 
-0. 146850836-03 
-0.15141 1826-03 
-0. 15603062E-03 
-0. 1571*7376-03 


-0.3994S536E-04 
-0.410830126-04 
-0.409510766-04  -0.422204026-04 

-0. 420915736-04 


SiaMA'22. 1 

0.1931 19076-06 
-0.825717206-06 
-0.178363066-05 
-0.277914996-05 
-0.373741256-05 
-0.465606096-05 
-0.554193976-05 
-0.640018856-05 
-0.723423466-05 
-0.804701916-05 
-0.884101296-05 
-0.961328806-05 
-0. 103803296-04 
-0.111296336-04 
-0. 118665336-04 
-0.125926676-04 
-0. 133089396-04 
-0.140161166-04 
-0. 147147716-04 
-0.154056326-04 
-0. 16091 3346-04 
-0.167737836-04 
-0.17*755536-04 
-0. 181627126-04 
-0. 188410486-04 
-0.213190546-04 
-0.238621566-04 
-0.263772296-04 
-0.288583236-04 
-0.313118696-04 
-0.337430706-04 
-0.361558826-04 
-0.385533886-04 
-0.409380196-04 
-0.4331 17036-04 
-0.456760106-04 
-0.480322356-04 
-0.503314256-04 
-0.527244636-04 
-0.550620876-04 
-0.573949206-04 
-0.597235026-04 
-0.620*32916-04 
-0.643o9675E-04 
-0 . 666  830266-04 
-0.690036296-04 
-0.713143306-04 
-0. 736093956-04 
-0. 759715626-04 
-0.  765243626-04 


SIGMA-22.2 

0.411934406-08 
-0. 109550226-05 
*0.213169096-05 
-0,312845396-05 
-0.408239686-05 
-0,499926906-05 
-0.588438886-05 
-0.674194676-05 
-0.757542396-05 
-0.838770416-05 
-0.91  SI  26546-05 
-0.995816456-05 
-0. 107201866-04 
-0.114689016-04 
-0. 122055966-04 
-0. 129314566-04 
-0, 136475156-04 
-0. 143546416-04 
-0. 150535916-04 
-0. 157450866-04 
-0. 164298636-04 
-0. 171088716-04 
-0.177829186-04 
-0. 184513146-04 
-0.191 106066-04 
-0.216482656-04 
-0.241995496-04 
-0,267148846-04 
-0.291960146-04 
-0.316496346-04 
-0.340809356-04 
-0.364938626-04 
-0.388914956-04 
-0.412762616-04 
-0,436500876-04 
-0.460145416-04 
-0.483709166-04 
-0.507202606-04 
-0.530634516-04 
-0.554012316-04 
-0,577342206-04 
-0.600629606-04 
-0.623679056-04 
-0.647094*56-04 
-0.670279296-04 
-0.693436506-04 
-0.716562966-04 
-0.739637736-04 
-0.762991 726-04 
-0.768759566-04 


SIGMA-22.3 

-0.23171 1056-06 

-0.140907636-05 

-0.252047606-05 

-0.351592446-05 

-0.446327666-05 

-0.537636736-05 

-0.625889156-05 

-0.71 1420676-05 

-0.794567016-05 

-0.875613606-05 

-0.9548071 16-05 

-0.103235026-04 

-0. 110842006-04 

-0, 1 18317016-04 

-0.125672866-04 

-0.132921276-04 

-0.140072346-04 

-0.147136206-04 

-0.154122226-04 

-0.161037466-04 

-0,167870316-04 

-0.174570646-04 

-0.181078806-04 

-0. 18757D24E-04 

-0. 193967896-104 

-0.219925706-04 

-0.245507326-04 

-0.270652346-04 

-0.295454686-04 

-0.319983616-04 

-0.344290636-04 

-0.368414926-04 

-0.392387086-04 

-0.416231246-04 

-0.439966586-04 

-0.463608676-04 

-0,487170356-04 

-0.510662086-04 

-0.534092556-04 

-0.557469186-04 

-0.58079S13E-0* 

-0.60*03*756-0* 

-0. 6273336 3E-C* 

-0.6505*3586-04 

-0.673732876-0* 

-0.6968S95IE-0* 

-0.720033776-04 

-0. /4322614E-04 

-0.766316026-04 

-0.772323346-04 


MOD  TR  59.22 


•a 


OUTPUTS 


MISSILE  fiT  LOCflTIOM  0, 


TIMES  0. 


POS 

V 

51 

0.8625 

52 

0.8887 

53 

0.8750 

54 

0.8812 

55 

0.8875 

56 

0.8937 

57 

0.9000 

58 

0.9082 

59 

0.9125 

60 

0.9187 

61 

0. 9250 

62 

0.9312 

63 

0.9375 

64 

0.9437 

65 

0.9500 

66 

0.9562 

67 

0.9625 

68 

0.9687 

69 

0.9750 

70 

0.9812 

71 

0.9875 

72 

0.9937 

73 

1.0000 

SIGMn-11.1 

-0. 155841 S5E-03 
-0. 15696640E-03 
-0. 158I0367E-03 
-0. I5924823E-03 
-0. I6039350E-03 
-0. t6153012E-03 
-0. I6284439E-03 
-0. 1S371873E-03 
-0. 16473688E-03 
-0. 16569748E-03' 
-0. 1S683840E-03 
-0. I6787I 84E-03 
-0. 189021 33E-03 
-0. 17103694E-03 
-0. 17414726E-03 
-0. 17868247E-03 
-0. 18449451E-03 
-0. 19032622E-03 
-0. 19297563E-03 
-0. 13650S36E-03 
-0. I8209376E-03 
-0. 10954704E-03 
-0. 3883I540E-04 


SIGMfl-l 1 .2 

-0.1570451tE-03 
-0.1581894e.E-03 
-0. 15930 305E -03 
-0. lfc044t58E-03 
-0. 16I57833E-03 
-0. 162704 17E-03 
-0. 16380789E-03 
-0. 16487895E-03 
-0. 16591480E-03 
-0. 16893483E-03 
-0. I6800074E-03 
-0. 189239048-03 
-0. 170S5266E-03 
-0. 17309892E-03 
-0. I7818I88E-03 
-0. I8005618E-03 
-0. I  840401  IE-03 
-0. t8632895E-03 
-0. t8350527E-03 
-0. I 7039808E-03 
-0. 1 4089864E-03 
-0.90595940E-04 
-0.21999229E-04 


SIGMA -It. 3 

-0. 15828757E-0.3 
-0. 159391 39E-03 
-0.160521 16E-03 
-0. 18165352E-03 
-0. 16278t63E-03 
-0.t63S9656E-03 
-0. 16493960E-03 
-0. 166057248-03 
-0.1671 I066E-03 
-0. 16St3999E-03 
-0. 16933076E-03 
-0. 17082372E-03 
-0. 172701 45E-03 
-0. I75I7442E-03 
-0. I7823413E-03 
-0. 18144771E-03 
-0. 18360373E-03 
-0. 18234968E-03 
-0. t7405236E-03 
-0. 1543C561E-03 
-0. I 1971200E-03 
-0. 71651922E-04 
-0.51689174E-05 


SIGM«-22. 1 

-0.77076984E-04 
-0. 77855 396E-04 
-0.78247329E-04 
-0.78843734E-04 
-0.79480151E-04 
-0.80080215E-04 
-0.80693390E-04 
-0.81235647E-04 
-0.81786008E-04 
-0.321 15758E-04 
-O.S2180347E-04 
-0.8I928806E-04 
-0.81445088E-04 
-0.81 1 31286E-04 
-0.81832263E-04 
-0.85243089E-04 
-0.93366335E-04 
-0. 10339836E-03 
-0. 13I09684E-03 
-0. 15816281E-03 
-0. 1 7834257E-03 
-0. 16795564E-03 
-0. 15532618E-03 


SIGMR'22.2 

-0.77452551E-04 
-0. 7802928SE-04 
-0.78806053E-04 
-0.79t32823E-04 
-0. 79759832E-04 
-0. 803384 19E-04 
-0.80912909E-04 
-0. 81 488509E-04 
-0.82062202E-04 
-0. 82632850E-04 
-0. 83198852E-04 
-0.83760244E-04 
-0. 84320541 E-04 
-0.8488S320E-04 
-0. 85480629E-04 
-0.88123643E-04 
-0.86848702E-04 
-0.87878991 E-04 
-0.88601844E-04 
-0.89523772E-04 
-0.9021 2221 E-04 
-0.9024t753E-04 
-0.87998916E-04 


SIGMfl-22.3 

-0.77832387E-04 
'0.7C407896E-04 
-0.78969458E-04 
-0.79521510E-04 
'0.30063729E-04 
-0.80597205E-04 
-0.81 131 979E-04 
-0.S1695893E-04 
-0.32342S79E-04 
-0.83t53994E-04 
-0.84221 373E -04 
-0.35598272E-04 
-0.87200390E-04 
-0.8S649755E-04 
-0.89083402E-04 
-0.37003670E-04 
-0.80335571E-04 
-0.669641 t3E-04 
-0.481  1 1406E-04 
-0. 20888883E-04 
-0.20849964E-05 
-0. 12529630E-04 
-0.20667869E-04 


MIS  W  59-22 


•27 


^OAiL'  HROB  3 


e/iffc/iy 


INPUT 

CLAbS  to  U 

6PHINO  CONS 


n*«KMONH.= 


mNu-  0«l-*aJt-00 


■»-  i.UOOUL-ui 


TA,.  O.SUt  0<* 

TAU  Uf  TAU  t»  TAU  ^ 

t.UOOUt-Ui  To. 


U.^Uc  Oa 


O.bsSst-OS 


IH10i;.=  O.AOuwt-u*:  niu=  O.baoUc-OU 


O.yot  ON 


KAD  n  59-22 


•28 


TkMH  l-UK  COCATlUN* 


l«U00OE-02 

l*OUOU£-w2 

1.UU0UE-U2 

UUU0UE-U2 

I•UU0UE*U2 

l.U00oe.02 

l«ugooE-o2 

UOOOOE-02 

uuuoue-02 

l.UUOUE-02 

UgU0OE-U2 

l«g00uE-o2 

1.OU0UE-U2 

l«goouE*g2 

i«uuogt-u2 

l.UOOOE-02 

UUOOUE-U2 

l«uOOJE>U2 

UUOOOE-02 

UOOOgE-U2 

1«U00U£-U2 

l*uoouE-a2 

I*uu0uc-u2 

l«U00oE-u2 


i»uggoE-*02 

l»0UOgE»O2 

l*wugwk-02 

l.OUOuE-02 

l*U00gE>02 

l,000uE-02 

l*0OUUE-O2 

l«0000E-02 

l.OOOOE-02 

l.CiU00l:>02 

l*gU0UE-02 

i«guUUE-U2 

l,OO0Uh-U2 

l,0UU0E-02 

l.UOUOE-02 

l«U000t-U2 

l*UU00t-02 

1,OOOUE-U2 

U0U0Uh-U2 

l«UU0uE-U2 

l.UuOUE-02 

l.UOOOE-02 

l,UOUUh-02 

i«go0uk-02 

i«00U0E-02 


i*00uUE*02 

l«OUo0t-02 

l«gguo&-02 

l«UuJuc-02 

l«0UgUt-02 

I»00o0t-02 

1*00UUE-U2 

l«0000E-02 

1«OOOOE->02 

l«00uoE-O2 

UUOuOE-02 

i«Q0o0E-02 

l«UUuOh-02 

l.OUGUc-02 

l«OUOUt-02 

l«UUU0t-U2 

X«UUUOk-02 

l«0UU0t-02 

I.OUUOC-U2 

l*OOUOt-02 

l,O0uUt-02 

l«OOOCb-02 

l«00u0t:>02 

l«OOUUt-02 

l.UOUOc-02 


l«OOOUt-02 

1«OOOUE<-U2 

l»OUUu&-o/ 

l«uuuuc~g2 

l*UU00t-M2 

l«gcuut-02 
1 •uoout^u^ 
l.OOOOt-02 
UOOOOt-02 
l«0000t-o2 
I«ooout-u2 
i  *0000^02 
i.oguoE^u2 
l.UUUOc-02 
I*0o00b-o2 
l«U00Ut-u2 
1.0000c -02 
I«0u00t-u2 
l«OOOOIr-u2 
1 .0000^-02 
1.0000t-o2 
1.0u00t-o2 
l.OOOOt-02 
i«uo00c-u2 


l.OOOOE-02 
l.OOOOE-02 
l*t0UUk-02 
l.OOOOE-02 
i*uOOUc-U2 
i«UU0OE-02 
l.UOuOt-02 
l«0UUOt-02 
UOOOOt-02 
l.UOOOt-02 
1.0000b -02 
l.OOUOb-02 
1.O000E-O2 
l.OOOOt-02 
l.OOOOc-02 
UOOOOt-02 
1.0000t-o2 

l.OOOOt-02 

l.OOOUc-02 

l.oOOOt-02 

l.o000fc-(«2 

.l.UOOOt-02 

l.OOOOt-02 

I.0000t-02 


l.OUOOE-02 
l.OOOOb-02 
1*U000L-U2 
l.OoUOc-02 
l.UOOOE-02 
i*OU00b-0^ 
l*OUO0t:-02 
l.UUOOE-02 
I«0000b-U£ 
l.UuJUt-02 
l.OUOOc-02 
l.OUOOt-02 
i*ui)Ooe-02 
l.OOOOc-02 
I .0UU0L-U2 
1 •UOUOt-02 
i • JOuoc-0^ 
I .OOUDb-O^ 
i.O000t-U2 
1  •OU00c**u2 
l.OOOOt-02 
1 • J000t“02 
1 .00001-02 
1  .0000^02 


1.0000t*02 

l*0000c*02 

lt0u00t-’02 

1.00O0t*U2 

i.OOOOt-02 

I.OOOOE-02 

l.OOOOc-02 

1.0000&-02 

l.OOOOc-02 

l.OUOOt-02 

l.OOOUt-O^ 

l.UOuOc-02 

i.ouuoc-u^ 

1.0U00L-U2 

I.OUOOE-02 

1.0000C.-02 

l.OOOOo-U^ 

1.00UOt-U2 

l.OOOOt-02 

i.0000t-U2 

l.OOOOc-02 

i.0000c:-02 

l.OOOOt-02 

i.UO0Ot-02 


l.OOOOE-02 

1«OOOOE-*02 

l.OOOOt-02 

l.OOUOb-02 

iiOOOOc-02 

l.OOOOE-02 

l.OOOOE-02 

l.OOOOE-02 

l.OOOOk-02 

i.OOOOc-02 

i.OOOOb-02 

1.0000t-U2 

1.0OU0c-U£ 

l.OUOUc-U^ 

i.0000c-U2 

x.uooot-u^ 

1.0000C-U2 
X.OUOUb-02 
1 .OOOOc-O^ 
i.ooo(;c-02 
l.OOUUc-02 
i.0000t-U2 
l.OOOOc-02 
l.OOOOh-02 


l*OO0OE*O2 
1*UOOOE-02 
l.OOOOb-02 
l.OUUOE-02 
l.OOOOt-02 
l.OOOOE-02 
l.OOOOfc-02 
l.OOOOE-02 
l.OOOOE-02 
l.OOOOc-02 
i.OOOUt-02 
i.OOOOc-02 
X.OOOoc-02 
l.OOOOc-02 
1.0000C'0<^ 
i.000oc”02 
l.U0Uoc-0<! 
l.UOOot-n2 
i.OUOOc-02 
l.OOOOc-02 
1 •OOOOc-02 
i.OOOOL-02 
1 . 0000c “02 

1 • 000uw“02 


OUTPUTS  MISSIUC  AT  tUCATlUN  1.  TlMt«  0. 


y 

w 

V 

X 

0«1U30 

-0,133bVAAAt-0i 

-^•1^04*96  ie6b*02 

Z 

0*1062 

-O.iSSbKAAAt-O^ 

-0*O9t>4^.96a>C'~02 

’i 

0*1126 

-0.iAau7b^bL-U< 

-W*b893436lt-U2 

4 

0*1187 

-O.ibObbHt/b-O^ 

-0*083036y6b«02 

S 

0*1260 

-Oi ibuAbbObt-O^ 

-U*O/bOi<*i6t-02 

6 

0*1312 

-0.ibU^V^4bb-U<: 

-O«b/06/162t-‘O4 

7 

0*1376 

-O.lSU^bb^bc-O^ 

-U«o643i:04yt«02 

b 
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FIG.  2  -  SPRING  RESTRAINT  SCHEME  FOR  BOUNDARY  CONDITIONS 
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Figure  3.  Spacing  of  Mesh  Points  Along  the  Generator 
of  the  Cone  for  n  *  24. 
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1.  INTRODUCTION  AND  PRESENTATION  OF  SAUDWICM  SHELL 
ASSUMPTIONS 

Following  i«  an  extenaion  of  the  analysis  of  homogeneous  shells  described 
in  Part  A  of  Chapter  IV  to  shells  of  sandwich  type  construction.  The  geometry 
considered  is  that  of  Figure  1;  i.  e. .  a  shell  of  revolution  with  core  material 
of  thickness  h^  and  with  an  outer  face-shell  of  thickness  h_  and  an  inner 
face-shell  of  thickness  h^.  The  coordinate  system  y^  y^  y^  is  the  same  as 
that  used  in  the  homogeneous  shell  analysis. 

In  accordance  with  the  accepted  practice  for  shell  theory  and  for 
sandwich  construction,  the  following  assumptions  are  made.  First,  in  the 


face  plates 

T*’  s  small  of  high  order  (1) 

s  small  of  high  order  .  (2) 

In  the  core,  infinite  rigidity  is  assumed  in  the  normal  direction 

i.e.,  T,,  =  0  (3) 

and  t’’  =  small  of  high  order  .  (4) 

Zero  rigidity  is  assiuned  in  the  tangential  directions; 

i.e.,  T**  =T**  =  0  (5) 

and  the  only  core  stresses  of  importance  are  ,  a  -  I,  2. 
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In  accordance  with  the  usual  assumptions  of  sandwich  theory,  the  face 

shells  are  assumed  to  be  thin  in  comparison  with  the  core  and,  therefore. 

3  3 

all  terms  depending  on  the  cube  of  the  face  thicknesses,  h^  or  h_  ,  are 
neglected. 
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2.  FORMULATION  OF  THERMOELASTIC  PROBLEMS  BY  MEANS 
OF  FREE  ENERGY 

Assxunptions  (1)  and  (2)  for  the  face  ihellf  imply  that  the  same  shell 
stress -strain  law  applies  to  the  face- shells  as  to  homogeneous  shells.  For 
the  tangential  stresses,  this  law  is  as  follows: 


l+v'  l-v  '  yb  l-v 

=  I,  2 

This  is  simply  Duhamel's  thermal  stress-strain  law  generalized  to  curvi¬ 
linear  coordinates  and  then  specialized  to  incorporate  the  vanishing  of  . 
In  the  core,  we  must  revert  to  the  more  general  relation  giving  and  then 
specialize  to  since  these  shear  stresses  are  the  only  core  stresses  of 
importance  in  the  analysis. 

1+  Vc  (7) 

where  use  has  been  made  of  the  geometric  relations  derived  in  the  homo¬ 
geneous  shell  theory  section  of  Part  A  of  this  chapter: 


*33  ~  I 


=  *03  =  0 


u  —  1»  2  . 


Now  the  thermoelastic  problem  for  sandwich  shells,  as  for  the  case  of 
homogeneous  shells,  is  formulated  through  the  use  of  the  so-called  free 
energy,  F. 


IrfAPD  TH  59-/2 


852 


This  free  energy  may  be  defined  in  general  as 

F  =  (1  J  dy*  dy*  dy*  t]. .  -  -  a^J  . 

2  JJj  V  iJ  I  -  2v  'ij/ 


In  the  present  context,  no  consideration  need  be  given  to  the  physical  nature 
of  this  quantity.  Rather,  F  may  be  thought  of  as  the  integral,  the  vanishing 
of  whose  first  variation  with  respect  to  displacements  leads  to  the  proper 
equilibrium  equations  and  natural  boundary  conditions  for  the  thermoelastic 
problem.  Now  from  the  assumed  vanishing  of  T**  it  follows  from  the 
general  thermal  stress- strain  law  that 


l+v  ( 


l-2v  [_ 


c  a  AT 
I  -  2v 


This  may  be  solved  for  t|jj  to  give: 


va  ^  _  (l+v  ) 

l-2v  '  l-2v 


-  V  ad 

n  —  ■  A  ^  n  + 

'jj  i-v  od 


f  H-v\ 

(^l-v  J 


a  AT 


853 


Substituting  this  into  the  integrand  of  (10)  and  eliminating  the  contribution 
of  T*’  ,  the  integrand  of  the  free  energy.  F.  becomes: 


-  r 

'ij  ’  ' 


<aAT 

l-2v 


'a^  l-', 


aP  1-1 


The  term  involving  (ATf  makes  no  contribution  to  the  variation  operation 
and.  hence,  is  dropped.  The  remainder  of  the  right-hand  side  of  (13)  reduces 


T'‘^  tT", 
aP  oJ 


c uAT  a 3 

TT-  *  • 


Thus,  the  free  energy  expression,  subject  only  to  the  restriction  that  =  0. 


is  given  by: 


=  TjjjJdyUy*  d/ 


«aAT  .  ^ 


Now  the  contribution  of  the  two  face- shells  to  F  is  clear.  Recognizing 
the  insignificance  of  in  the  face' shells,  this  contribution  becomes: 

F,.c.  .h.u.=n*.-(..4f j  nj 

'^’^•4ace  shells'  volume  _ ' 
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<10 

For  the  core,  the  term  T  tj^^rope  out  of  the  integrand  of  (15)  due  to 

(&  0 

the  aisvimed  inaignificance  of  core  in-plane  atreaaes  T  .  However,  it  ia 
neceaaary  to  evaluate  the  core  in-plane  atraina  which  enter  the  thermal 
term  in  (15).  Thia  ia  done  aa  followa:  Firat,  aince  ia  aaaumed 
inaignificant  in  the  core,  the  thermal  atreaa-atrain  law  for  the  core  in-plane 
atreaaea  ia  the  aame  aa  for  the  face  ahella  (6)  but  with  core  propertiea  <c  • 

d  0 

replacing  c,  v  .  Thua,  the  vaniahing  of  T  in  the  core  impliea: 


*c  ®c  * 


aj3 


AT 


I- 


(17) 


We  require,  for  F,  the  expreaaion  for  .  Thia  ia  found  by  multiplying 

( 17)  through  by  a^^  ,  contracting  on  a  and  recalling  that  the  metric  tenaora 
are  aymmetric  and  that  the  contra  variant  and  covariant  forma  of  the  metric 
tenaora  are  inveraea  of  each  other:  i.  e. , 


=  0 


=  I 


Y  =P 


(18) 


Performing  thia  multiplication  and  contraction  we  find  (noting  that  6^  =  2 
aince  we  are  effectively  in  2-apace  when  amnming  over  Greek  indicea) 


.-L.  .op, 

l-v^.  l-Vc  /  t6  l-Vc  * 


(19) 
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Thuf,  we  have  the  required  expreseion  core. 

Now  this  together  with  (7)  is  introduced  into  expression  (IS)  to  get  the 
core  contribution  to  F. 


core 


.1  ffr 
03 


r 


J  dy^  dy*  dy* 


l+v. 


ore 


^3^X3 


ZcCagAT)^ 
I-  Vc 


^core  "  \  1  \ 

(effective)  JJjcore 


J  dy*dy2  dy* 


(20) 


vdiere  the  (n^T)^  term  is  dropped  since  it  makes  no  contribution  in  the 
process  of  forming  a  first  variation  with  respect  to  the  displacements. 
The  final  effective  expression  for  free  energy  F  is  thus  given  by: 


^  “^face  shells  ^  ^core 


(21) 


or 


.  TVl  Jdy>dy*dy>/'ll-  + 

volume 


aj3 

l-v 


+  .  I"  J  dy‘d/ 


(22) 


afi 

*  n  n 

2(1+ V.)  ®>  PJ 
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Affuming  face  ehell  straina,  temperatures,  and  metric  elements  to  be 
replaceable  by  their  mean  values  (at  the  face  shell  mid-surfaces)  and  using 
(t)  notation  to  denote  inner  and  outer  face  shells,  respectively,  this  becomes 


If  the  face  shell  thicknesses  h^.^  vary  along  the  meridian,  it  is  merely 

necessary  to  keep  these  quantities  inside  their  respective  integrals. 

Thus,  in  the  core,  integration  through  the  core  thickness  (with  respect 
to  y^)  is  required  while  in  the  face  shells,  mean  values  (values  at  the  face- 
shell  mid- surfaces)  are  used.  More  will  be  said  subsequently  about  the 
meaning  of  the  face- shell  mean  temperatures  ^'^-) 

called  for  in  the  free  energy. 
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f 


) 


1  2 

The  pofsible  presence  of  external  pressure  p(y  .y  )  reqxiires  the 
formulation  of  a  work  term  W  as  follows: 

w  =  j  dy*  dy^  'JCy^.y*)  p(y*.y*)  u*  (y*.  y^)  •  (24) 

The  variational  formulation  of  the  problem  of  sandwich  shell  response 
to  both  mechanical  (pressure)  and  thermal  loads  is  then  simply  the  following: 


6  (F  -  W)  =  0  (25) 
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3. 


DISPLACEMENT  AND  STRAIN  EXPRESSIONS.  AND  REDUCTION  TO 


A  TWO-SPACE  PROBLEM 

In  accordance  with  usual  shell  theoretical  assumptions,  it  is  presumed 
that  the  tangential  physical  displacements  are  approximable  by  first-order 
polynomials  in  y^  and  the  normal  displacement  is  independent  of  y^. 


=  “(o)  • 


The  relation  between  tensor  and  physical  displacements  is: 


(26) 


u^  = 


11 


(27) 


no  sum  on  i 

wdiere  is  the  contravariant  form  of  the  displacement  vector  .  Thus, 


no  sum  on  a 


(28) 


u  =1 


The  strain-displacement  expression  within  the  assumption  of  linear 
elasticity  is 


I 


’^ij  2  j  i  ^ 


(29) 


commas  denoting  covariant  derivatives.  Following  from  (28)  and  (29), 
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together  with  the  orthogonality  of  the  curvilinear  coordinate  system  and  the 
independence  of  metric  elements  a^j^j  on  y^,  it  can  be  shown  that: 


,  ^  ±  ^  iui 

“  11  ft,,l  ?  *  tk\r^  ? 


2  By* 


-  .  r 

—  Bu^ 

u*  Be22 

^22 

-  Va22  - —  +  . 

By* 

2  By^ 

n  = 

U 

1 

2 

By* 

- 

4 

Ba,, 

Byi 

•  +  J- 
2 

■ 

T1  = 
13 

1 

2 

y  Bu* 

+  1 

2  By* 

_v^ 

4 

Ba^^ 

By* 

»l23 

.  1 

2 

y—  BU* 
'/*“  8/ 

1  Ini 

2  8y2 

ui 

4 

Ba2  2 
By* 

Bu> 


Substituting  into  (30)  from  (26)  the  strains  are  expressed  in  terms  of 


and  u^ 


u  'V  “  9yi  Qyi  J  2  By^  -^(o)  (ly  2  By^ 


*y’  -^V  ^  ^ 

.By*  177  2 


“  .,8y'  8y‘j'  “xSy'  Sy'y'  <  By*  l,'°>  V 


.  =T  i'*..  V,yT 

»3  2  \  (  i)  2  Qyi  ^  gy3 


U*  +  y» 

(o)  ^ 


u*  ) 

(,)/ 


_  I  -  2  I  Bu/  V  I  a‘ 


®®22  2  .3  2  ^ 

— T"  1  V  +  y  '3» »  > 

8y*  (o)  (ly 
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The«e  expressions  for  strain  now  permit  the  evaluation  of  ii . .  and 

'J(+) 

Ti-.  in  an  obvious  manner.  As  in  the  expressions  for  geometrical 

..  3 

quantities  aV  ,,  a^J  ,  etc. ,  it  is  merely  necessary  to  substitute  for  the 

(+)  (-) 


values 


-(^5  the  respective  integrands  of  and 


Now  the  metric  tensor  elements  a-^j  and  a^'^  are  known  functions  of 
yl,y^,  y^  as  is  the  Jacobian  J.  Thus,  integration  with  respect  to  y^  through 
the  core  thickness  is  possible  in  the  core  contribution  to  F.  The  core  con- 
tribution  may  be  integrated  with  respect  to  y’  either  analytically  or  by 
quadratures,  whichever  is  more  convenient  in  the  numerical  program.  The 
result,  in  either  case,  is  a  formulation  (25)  in  vdiich  now  both  F  and  W  are 
2-space  integrals  (over  y^  and  y^). 

It  is  important  to  note  that,  unlike  the  homogeneous  shell  analysis,  the 

sandwich  shell  theory  places  no  restriction  on  the  shear  strains  in  the 

a 

core  material.  Thus,  the  displacement  slopes  u^^^  cannot  be  evaluated 

a  priori  in  terms  of  the  u^^^^  but  must  be  found  as  a  result  of  the  overall 

calculation.  This  implies  that  five  functions  u^,  ,  uf’  ,  u*  are 

(o)  (o) 

being  sought  in  the  sandwich  problem,  whereas  only  three  were  required  in 
the  homogeneous  shell  problem.  Hence,  variations  must  be  performed  with 
respect  to  all  five  of  these  functions  in  the  present  case.  Furthermore,  it 
should  be  noted  that  AT  appears  only  in  the  face- shell  contributions  to  F, 


SO  that  it  is  unnecessary  to  introduce  the  7'  '  j'  '  7 


integral 


expressions  which  appeared  in  the  homogeneous  shell  analysis. 


jvadd  tr  s)-..: 


WAU)  TR  59-?2 


862 


CP 


11 


C08  m  8 

m  'll 


^122  =  >  xn^22  cos  m  6 


12 


sin  m  8 

III 


■n  sin  m  6 
m  o 


Tl2j  =  ^  ^Tijj  cos  m  8  .  (34) 

Analogous  expressions  apply  for  the  corresponding  stress  components,  i.  e.  . 

'rll 

I  IS  a  cosine  series,  etc. 

The  detailed  expressions  for  the  strain  components  are  as  follows: 


m 


m 


^  ^  8y 

.  y>  MlA  +  ^ 

\  8y^  0/  /  2  gy» 


m 


m 


1  •  ”  u’  -  fSL  /  1  .  .  1  \ 

11  2  r  n.  ^ 


2  \ 

V  0y* 

®*ii 

4 

0/ 

(35) 
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As  may  be  seen  by  inspection  of  the  free  energy  expression  (23)  and 
the  work  integral  (24),  the  harmonics  all  decouple  and  it  suffices  to  consider 
each  Fourier  component  separately.  The  m  =  0  component  corresponds  to 
axisymmetric  loading:  m  =  1,  to  the  simplest  case  of  asymmetric  loading, 
and  so  on.  The  linearity  of  the  analysis  then  permits  superposition  of  the 
Fourier  components  with  the  appropriate  cos  m  0  or  sin  m  0  multipliers  to 
give  total  displacements,  strains,  or  stresses  as  required. 
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5.  NUMERICAL  ANALYSIS 


The  eandwich  ehell  response  to  thermal  and/or  mechanical  loading  is 
to  be  solved  starting  from  the  formulation  (25)  together  vdth  expressions  (23) 
and  (24)  for  F  and  W>  respectively.  As  noted  earlier,  the  core  contribution 
to  F  is  integrated  numerically  or  analytically.  This  leaves  the  expression 
F-W  in  the  form  of  a  two -space  integral  with  integrand  involving  the  known 
geometry  (aj^j,  etc. )  and  the  five  unknown  functions  of  y  and  y^ 


^  %)*  '"o)  ’  "(o)  )  .  Resolution  into  Fourier  components  then  reduces  the 


(F-W)  integral  to  a  set  of  one-dimensional  integrals  over  y^,  the  meridional 
variable,  and  reduces  the  entire  analysis  to  a  set  of  uncoupled  problems. 


for  the 


6  (  F  -  W  )  =  0 
m  m  ' 

functions  (  „u® . ,  —U®”  ,  u*  ) 

'  m  (o)  m  (i)  ’  m  (o) 


(36) 


This  may  be  handled  by  any  one  of  a  number  of  techniques:  (a)  A  direct 
variational  method  may  be  used  in  which  each  of  the  unknown  displacement 
functions  is  expanded  in  a  series  of  boundary- condition  satisfying  functions 
and  the  variation  is  formed  with  respect  to  the  coefficients  of  these  functions, 
(b)  Alternatively,  the  direct  variational  method  may  be  applied  with  respect 
to  the  point  values  of  the  displacements  and  displacement  slopes  at  a  set  of 
mesh  points  along  the  shell  meridian.  (This  corresponds  to  the  method  udiich 
has  been  programmed  for  homogeneous  shells.)  (c)  Another  possibility  is 
to  form  the  required  variation  analytically,  thereby  deriving  the  differential 
equations  and  natural  boundary  conditions  of  the  problem.  These  could  then 
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be  folved  numerically.  Inasmuch  as  diis  corresponds  closely  to  the  direct 
variational  method  (b)  and  method  (b)  insures  the  positive  definite  nature  of 
the  matrix  for  the  difference  equations,  method  (b)  is  preferred.  Furthermore, 
method  (b)  eliminates  the  necessity  for  deriving  the  Euler  equations  analytical¬ 
ly  and,  at  the  sanie  time,  automatically  introduces  natural  boundary  conditions 
at  free  or  simply  supported  edges  in  proper  finite-difference  form.  Under 
this  method,  the  required  integrals  are  evaluated  by  quadratures  in  terms  of 
the  mesh-point  values  of  the  unknown  functions.  The  result  leads  directly  to 
the  matrix  and  the  forcing  vector  for  the  numerical  determination  of  these 
point  values.  The  matrix  algebra  problem  is  then  solved  for  these  point 
values,  and  strains  and  stresses  are  calculated  from  them  using  finite 
difference  approximations  to  the  required  derivatives. 
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where  =  T  (  x„j  ,  vAt)  . 

ni 

The  heat  fl\ix  Q  depends  on  time,  on  (y^  y^),  and  on  the  outer  wall 

V 

temperature  itself.  The  quasi-one-dimensional  analysis  must,  therefore, 
be  performed  at  a  number  of  space  points  (y^  y^)  depending  upon  the  nature  of 
the  physical  problem.  If  only  axisymmetric  conditions  are  of  interest  in  a 
particular  application,  then  no  y ^-dependence  of  Q  occurs.  The  reason  that 
Fourier  analysis  is  not  suggested  for  the  heat  problem  lies  in  the  essentially 
non-linear  character  of  the  outer  wall  boundary  conditions;  i.  e. ,  radiation 

V  ^ 

is  frequently  important,  and  this  involves  (  T^  )  ,  and  convective  heat 

V 

transfer  also  involves  T^  in  a  non-linear  manner.  Thus,  the  functions 
AT+  and  AT_  must  be  calculated  first  as  functions  of  (y^,  y^,  t  )  and  then 
the  quantities  aAT^  ,  aAT_  may  be  Fourier  analyzed  giving  Fourier 
components  which  depend  on  y^  and  time.  The  thermal  stress  analysis  is 
then  carried  out  at  a  number  of  fixed  values  of  time. 

A  mesh  is  chosen  for  the  core,  and  the  heat  equation  for  the  composite 
material  is  written  in  difference  form.  These  difference  equations  may  be 
put  in  explicit  difference  form  with  a  strict  upper  limit  on  permissible  time 
step  At  for  stability.  Alternatively,  they  may  be  written  in  implicit  form 
with  no  stability  limitation  on  At  .  The  implicit  form  requires  solution  of 
(N  -t  1)  simultaneous  equations,  but  the  tri- diagonal  nature  of  the  associated 
matrix  permits  rapid  and  efficient  solution.  A  more  serious  problem  which 
arises  in  connection  with  the  implicit  scheme  is  the  necessity  to  iterate  in 
order  to  express  heat  flux  Q  at  the  latest  time.  This  is  not  a  major  drawback. 
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however,  and  much  numerical  experience  has  been  accumulated  in  this  respect. 
The  large  gain  to  be  realized  in  time- step  size  At  recommends  the  use  of  the 
implicit  scheme.  The  difference  equations  for  the  heat  problem  are  well- 
known  and  need  not  be  reproduced  herein. 
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